物理学入門

　現在我が国の大学教育は一般教育と専門教育とに分かれているが、学問とはどのようなものであるかを学び取らせるということを第一の目的としている一般教育は専門教育の領域に一歩足を踏み入れなければその実行は不可能である。このことを考え、なかば専門教育の意味も含め、しかも対象を理工科方面に進む学生に絞って書いた一般物理学の教科書がすなわち本書である。

　本書はその性格からいって、物理学科もしくはこれに近い学科に進もうとする学生諸君にとっては物理学の入門編であり、その他の諸君にとっては完結編となるべきものであるから、執筆者はいずれも基礎的部分に重点を置くとともに、現代の情勢からみて、科学技術者としての必要な知識はできるたけ盛り込むように努めた。内容は多方面にわたり、多少専門領域にまで渡った部分もあるので、この全部を限られた授業時間内に教室で完結することは不可能である。教室で取り上げるのはこの中の基礎的部分だけというつもりである。全体を通じ、学生諸君の自習を期待したので、説明を平易にすると共に数式の問題などもなるべく途中を省かないように努めた。また、随所に例題を入れ、章末に演習問題を設け、さらにその解をも示しておいた。このように本書は教科書という従来の通念から離れた自習書である。自習の間に必ず湧いてくるであろう疑問はこれを教室において解決するという方法で学習の実を挙げていただきたい。実際このような方法をとらない限り、現在のように進歩した物理学を一般教育課程というわずかの時間内に修めることは不可能である。

　物理学の学習上いつも支障となるのは学術用語の不統一と、単位系の混用である。このような本質的でない問題に学力を割くのを防ぐため、本書では用語はなるべく文部省編集による学術用語集によることにし、単位記号や量をあらわす文字は国際純粋応用物理学連合（IUPAP）および日本工業標準（JIS）調査会の奨めているもの（SI単位系）に従うことにした。ただ、自習するにあたって過去の文献や論文を確認する必要が生じた際にMKS単位系、電磁気学ではMKSA単位系で表記されていることがある。そのような場合、用語や実験データを引用する際に現在利用されているものに変換が必要なことに注意が必要である。（参考文献：物理実験者のための13章）

物理学における研究法

　物理学では自然界におけるいろいろな現象を取り扱うが、その扱い方に特徴がある。すなわち、まず現象を分析してその中から大きさの考えられる要素を抜き出す。そして現象をこの要素の間の数量的関係として表す。例えば、空中に弧を描いて飛んでいるボールの運動という現象は、これをボールに与えられた”力”、これに働く”重力”、ボールの”速度”、”位置”、飛んでいる”時間”というように全て大きさの考えられる要素間の関係として捉える。このように大きさの考えられる要素を物理量あるいは単に量という。そして物理法則は一般に物理量の間の数量的関係として表される。物理量の中には例えば体積のように、その大きさだけしか考えられないもの、また力のように、その強弱（大きさ）の他に方向も合わせ考えなければなければならないものなどがあるが、いずれも物理量であることに変わりない。

　さて、一つの物理量を扱うには最初にその性格をはっきさせた上で、その大きさを決めるための物理的な方法を規定しなければならない。例えば、物体の重さという量を考えられるときに、その重さを単に手の上に乗せた時の我々の感じで決めるとすることはできない。なぜかといえば、感覚はそれほど精密でもなく、またそれほど信頼のおけるものでもないからである。どうしても何か基準となるはかりのようものを考え、はかりに現れる目盛りによって重さを表すというような方法によらなければならない。

　ところで、実際に自然界に起こっている現象というのは、これら多くの物理量が互いに関連しあい、複雑にからみあったものである。この中からある要素だけを単独に他と切り離して取り出すことは実際にはできない。例えば物体の運動に関して距離・時間・質量・力などの物理量が関係してくるが、この中から例えば時間という量を決める手段を考えてみよう。我々には時間の大きさというものが分からない。そのため、直感的に大小のわかるような長さあるいは角度のような1次元の空間的広がりを持つ量を用いて時間の測定を行う必要性が出てくる。加えて一定の速さで動くものを用意する必要もあるが、ここで一定の速度をもつことを確認するため事前に運動に関する法則を知る必要がある。しかし、まだ時間の決め方が確定していない以上その運動というものが決められるべきものでない。このように、ある一つの量を定義する際に他の量を定義しておく必要があった場合に、２つ目の量を定義するのに1つ目の量が関連してくるという循環的な関係が生ずる。

　この関係は、非常に数多くの未知数を持つ複雑な連立方程式の解を求めるときの事情に似ている。極めて多元で高次の方程式を解くとき、正攻法で各未知数を順次に求めることはできない。こういう場合、我々は視察によって解くという方法を用いる。すなわち、まず数個の未知数に予想される数値を仮に与え、他の未知数はこれから求める。そして次にこれらが全ての連立方程式を果たして満足するか否かを調べる。もし満足しなければ最初の予想値を修正する。こうして略値から漸近的に正しい答えに向かうのである。

　これを上の運動の例で言うと、まず地球が等速回転しているものと仮定してから時間の大きさの決め方を決める。次に適当な予想を許すことにより質量や力などの量の決め方を定め、運動の法則を導く。そして、この運動の法則を基礎として地球の自転運動を論じ、これが等速回転運動であることが証明されれば最初の予想に過ちが無かったということが言える。このような経緯で、物理量はまず最初にはっきりした定義を与え、明確な測定方法で一つ一つ順次に規定していくというわけにはいかない。どうしても、ある仮定を設けて進んだ後で訂正を行いつつ全体としてみた時にそこに矛盾がないように決めるという方法に頼らざるを得ない。ことに力学の最初の部分などにおいて、個々の量の決め方に論理的な必然性を欠くところがあるが、これらの量は全体の関連において定義され理解されるべきものである。

　以上のことは法則を導く場合でも同様である。我々が数個の要素の間の関係を知ろうとして実験しても、実際には必ずほかの要素が不可逆的に混入してくる。しかし我々はこの実験結果から混入要素を除いた時の結果を想像して、これから与えられる要素間の関係すなわり法則を帰納する。このように、物理学における研究方法は帰納的であって演繹的ではない。したがって仮に現在完全と思われる物理学の体系があったとしても、将来の新事実が果たしてこの中に矛盾なく収まるかどうかは保証できない。入らなければ体系を作り変える。こうして次第に完全な体系へ向かわせようというのが物理学の研究方法である。

単位と単位系

　一つの物理量を表すのには、まずその基準となる大きさの量すなわち単位を適当に選び、その物理学がこの単位の何倍（数値）にあたるか示せばよい。すなわち $物理量 = 数値 \times 単位$

の関係で表すことができる。

　さて、単位はそれぞれの物理量に対して決めるべきものであるから、それぞれ任意に選ぶことができる。例えば、長さと面積とは本来別種の量であるから別々に単位を決めて悪いことはないが、ある関係を持たせて決める方が現象を数量関係として表す場合に簡単になる。例えば長さの単位を任意に決めたら、面積の単位はこの長さの一辺とする正方形の面積と決める。また体積の単位としては単位長さの一辺とする立方体の体積をとり、さらに長さと時間の単位を適当に決めたら、速さの単位は単位時間に単位の長さ進む速さと決める。このように、最初まったく適宜に決めた単位を基本単位、これから導かれた単位を組立単位（誘導単位）という。問題は物理量全体を取り扱うのに何個の基本単位が必要かということであるが、力学的現象を扱う範囲内でいえば3個の基本単位が必要で、またそれだけあれば十分である。このことは実際にあたってみれば容易に理解されるであろう。基本単位はそれが絶対不変であると信じられるものでなければならない。そしてこの条件を満たすものとして長さ・質量・時間の3個が選ばれている。このようにして長さ・質量・時間の基本単位およびそれから導かれる組立単位の全体を絶対単位系という。

　単位自体、各人が自由に作り使用することができるが、いざ量の大小を比較するとなると単位間の換算まで都度行わなければならない。そのため、以前からそれぞれの集団や壮会において統一した単位を使用する動きがあった。世界的には、1875年に17カ国の批准によって成立したメートル条約から始まり、数年に一度ずつ国際度量衡総会（度は物差し、量は升、衡は天秤）が現在に至るまで開かれている。この中で決議され、各国で正式に使用が認められている単位系を国際単位系（Système International d'Unités 略して SI単位系）といい、7つの基本単位および2補助単位（平面角のrad(ラジアン)と立体角sr（ステラジアン））で構成される。下記は、2023年時点でのSI単位系における基本単位の定義を表している。（以下も古いため、修正が必要）

時間：秒（s）
セシウム周波数 $Δ ν_{Cs}$ 、即ちセシウム133原子の摂動を受けない基底状態の超微細構造遷移周波数を単位Hz（ s $^{- 1}$ に等しい）で表したときに、その数値を9192631770と定めることによって定義される。
長さ：メートル（m）
真空中の光の速さcを単位ms $^{- 1}$ で表したときに、その数値を299792458と定めることによって定義される。ここで、秒は $Δ ν_{Cs}$ で定義される。
質量：キログラム（kg）
プランク定数hを単位J・s（kg・m^ $^{2}$ ・s $^{- 1}$ に等しい）で表したときに、その数値を6.62607015×10 $^{- 34}$ と定めることによって定義される。ここで、メートル及び秒は、それぞれc及び $Δ ν_{Cs}$ を用いて定義される。
電流：アンペア（A）
電気素量 e を単位 C（ A・sに等しい）で表したときに、その数値を1.602176634×10 $^{- 19}$ と定めることによって定義される。ここで、秒は $Δ ν_{Cs}$ によって定義される。
熱力学的温度：ケルビン（K）
ボルツマン定数 k を単位 J・K $^{- 1}$ （ kg・m $^{2}$ ・s $^{- 2}$ ・K $^{- 1}$ に等しい）で表したときに、その数値を1.380 649×10 $^{- 23}$ と定めることによって定義される。ここで、キログラム、メートル、秒はそれぞれh、c、 $Δ ν_{Cs}$ を用いて定義される。
物質量：モル（mol）
1モルには、厳密に6.02214076×10 $^{23}$ の要素粒子が含まれる。この数は、アボガドロ定数 N $_{A}$ を単位mol $^{- 1}$ で表した時の数値であり、アボガドロ数と呼ばれる。系の物質量(記号は n )は、特定された要素粒子の数の尺度である。要素粒子は、原子、分子、イオン、電子、その他の粒子、あるいは粒子の集合体のいずれかであってもよい。
光度：カンデラ（cd）
周波数 540×1012 Hz の単色放射の視感効果度 Kcdを単位 lm・W $^{- 1}$ （ cd・sr・W $^{- 1}$ あるいは cd・sr・kg $^{- 1}$ ・m $^{- 2}$ ・s $^{3}$ に等しい）で表したときに、その数値を683と定めることによって定義される。ここで、キログラム、メートル、秒はそれぞれ、h、c、 $Δ ν_{Cs}$ によって定義される。

　物理学で取り扱う量は膨大なものから極微のものに至る。長さで言えば、天文学的の距離から原子的な長さまである。そこで一つの単位に対して副単位を考えておくと便利である。実際に用いられている副単位は基本単位名に以下のような接頭語をつけて呼ぶ。例えば nm は nanometerで $1 0^{- 9}$ m を意味し、Tg は teragram で $1 0^{12}$ g のことである。

基本単位に対する倍数	呼び方	記号
$1 0^{12}$	tera	T
$1 0^{9}$	giga	G
$1 0^{6}$	mega	M
$1 0^{3}$	kiro	K
$1 0^{- 1}$	deci	d
$1 0^{- 2}$	centi	c
$1 0^{- 3}$	milli	m
$1 0^{- 6}$	micro	μ
$1 0^{- 9}$	nano	n
$1 0^{- 12}$	pico	p

長さの単位

　長さはもともと点の相対的位置を表すための量で、二つの長さはこれらを重ね合わせることによって直感的にその大きさを判定することができる。長さがあらゆる物理量の中での最も基礎的なものになっている理由はここにある。そして他の物理量も、これを正確に知るためには長さの測定によることが非常に多い。例えば、時間は時計の針の動いた距離で、物体の重さははかりの目盛りで温度は水銀柱の長さで測るなどがこれである。

　長さの単位として物理学で用いられているものはメートルであるが、このメートルが最初測定されたのはフランスで1795年のことである。そして地球子午線の全長の4000万分の1の長さをもって単位とするというのが最初の考えであった。これにしたがって子午線の実測を行い、その後、幾多の経緯を経た後に1875年に至って初めて国際的にメートル条約ができ、メートル原器というものが制定された。このメートル原器は特殊の切り口（これは曲がりにくくするため）をもち、白金90％、イリジウム10％の合金で作られた棒で、その両端に近いところそれぞれ細い標線がある。0℃におけるこの標線間の長さが1m と決められた。この国際原器はパリ郊外の国際度量衡局に保管され、各国はこれとほとんど同じもの（副原器）を保管することとなった。それぞれ多少の誤差は免れ難く、我が国に配布されたものでは、標線間の長さは

$1 m - (1.3 - 8.667 t - 0.00100 t^{2}) μ$

であった。ただし $t$ は温度（℃）を、 $μ$ はミクロン（micron）で $1 0^{- 6}$ m である。

　このようにして決定されたメートル原器にもいろいろ難点がある。まず第１は永久不変性の問題で、人工的なものである以上、破損の危険がないとは言えないし、またこのような合金製のものが長年月の間にひとりでに変形しないとは保証できない。しかもそれ以上の欠点とされるのは精度である。原器の両端にある標線は細い線であるが、それでも幅が $68 μ$ 程度である。そこで両線の中心を採ったとしても全長に $0.2 μ$ 程度の誤差は免れ得ない。すなわち、相対誤差 $2 \times 1 0^{- 7}$ の精度しかなく、これでは進歩した現代の精密科学の要求を満足してくれない。そこで、永久不変性という意味で、ある決まった光の波長を長さの基準にするという考えが20世紀初頭からあった。そして、 Cd 蒸気中における放電の場合に発する赤色光をこのために用いるという提案がなされたこともある。

光の波長を基準にするというのは

"外から何らかの影響を受けていない静止した同一種の原子の出す光の波長は定まっている"¹

という考えに基づくものである。しかし現実には、元素が同位体の集合であって、完全に同一種の原子の集合でないこと、原子同士が互いに影響しあっていること、これを光らせるために電気的影響を受けていること、原子が前後左右に振動しているために、その発する光が例え一定でもDoppler効果によって異なった波長の光の混合として観測されること、などその他いろいろの理由で純粋に波長が一定した光というものが得られなかった。

　その後の研究によって $^{86}$ Kr² という同位体が分離されるようになり、理想に近い光が得られるようになった。そして1960年10月14日の国際度量衡総会で、今後長さの単位として $^{86}$ Kr の発するKr66橙色線の波長を用いることが決定され、わが国でも1961年7月1日からこれによるということが政令で決められた。

以下の図はこの $^{86}$ Kr ランプの要点を示したもので、Vは熱陰極を持った真空管で、この中に少量の $^{86}$ Kr が封入され、放電によって得る光の波長はBの方向から分光器で観測する。Vを液体窒素Nの中に入れ、この窒素はポンプで減圧し、Nの温度が（63±1）°Kの範囲になるようにして用いる。

Kr86ランプ

$^{86}$ Kr の低温低圧における発光を用いたのは、上に述べた種々の障害を少なくするためで、放電電流にも一定の限度が与えられている。

こうして出てくる光の中のいわゆる Kr86 橙色線の真空中における波長の1650763.73倍を 1 m と定義する。もちろん従来の 1 m とほとんど違いはない。

この方法は極めて精度の高いもので、メートル原器の精度が $2 \times 1 0^{- 7}$ であったのに、ここでは $1 0^{- 8}$ の精度が得られ、将来は $1 0^{- 9}$ の精度にまで上げられる見込みである。

数個の飛び飛びの波長の光の混合でるから、その中の一つを選ぶ。

もともとKrの中には質量数、78、80、82、83、84、86の同位体が含まれている。

時間の単位

　我々には時の経過したことは分かるが、経過した時間の量は精密には分からない。そこで、一定の周期をもつと想像されるもの、あるいは等速運動をすると想像されるものを仮定して時間を測る。そしてこの想像が果たして正しかったかどうかは、これから出発して作られる物理学の体系を全体的に見て判定する。

　それにしても、ある一定の時間、すなわち単位を決めなければならないが、最初の一定の時間として採用されたものは地球の自転から割り出したものであった。いま、地球の中心を中心としてきわめて大きな半径の球面を天空に想像し、全ての星をこの天球上に投影してみる。地球の自転を考慮しないことにすれば、全ての恒星はほぼ天球上の固定点となるが、太陽はこの天球上を運動する。それは地球が太陽の周りに公転しているからである。

この太陽の軌道 VDACV を黄道という。地球は太陽の周りに一定の角速度で運行しているわけではないから、黄道上の太陽の速さも一定ではない。また地軸はこの黄道面（すなわち地球の公転軌道の平面）に垂直になっていない。したがって、地球の自転角速度が一定としても、１日の長さ、すなわち太陽が南中してから翌日南中するまでの時間は四季によって変わってくる。

　そこで、地球の中心を通り、地軸に垂直な平面と天球との交線 EW を赤道と名づけ、この赤道上を一定の速さで動く仮想的な太陽を想像して、これを平均太陽と名付ける。赤道と黄道との交点 V, A の中で、太陽が南半球から北半球へ移る方の点 V を春分点、他の点 A を秋分点という。平均太陽の運動の速さは、真太陽と同時に春分点を出発して、翌年同時に春分点に戻るような速さと決められている。実は地軸が一定の方向を保たず、きわめて徐々にこまのようなすりこぎ運動をしているので、赤道面の方法は徐々に変わり、春分点も多少移動していく。春分点を出て再び春分点に戻るまでの時間は、厳密に太陽が天球上を１周する時間よりやや短く、これを 1太陽年（あるいは 1回帰年）という。

　それはとにかく、平均太陽は赤道上を一定の速さで動いているので、もし地球の自転の速さが一定ならば、平均太陽が今日南中してから、翌日南中するまでの時間は一定となるはずである。この時間を1日、その 1/24 を 1時間と定義する。¹

　以上が従来の時間の単位で、これは地球自転が一定不変の速さという仮定に基づいている。しかし。このようにして決めた今までの時間を用いると、地球から見た月の運動が計算から予想されるところと違い、永年の中には次第に早くなっていく傾向のあることが分かった。地球に万有引力をおよぼしているものは太陽の他に惑星もある。このために地球は太陽の周りに正確な長円運動をするのでなくて、きわめて僅かずつであるが、軌道の形が円形に近づいていく。月のみかけの運動が計算と合わない一つの原因がここにあるが、それだけでは計算の違いが十分説明できない。このようなことから、地球自転速度が一定という仮定にはかなり以前から疑いが持たれていた。近年、時計が非常に進歩し、ことに水晶時計が発達し、時間の測定がきわめて精密に行われるようになるにつれてこの疑いがますます濃くなってきた。

水晶時計

　水晶は、伸び縮みするとその両端に正負の電荷を発生する。そこで両端に交互に正負の電荷を与え、しかもその周期に合わせると共鳴してよく振動する。この振動の周期は非常に一定したもので、温度による影響などもはなはだ小さい。そこでこの水晶片を利用して一定周期の交流を作ることができる。ただこの場合の交流の周波数は極めて大きいが、周波数は真空管回路によって何分の一にでも減少することができる。こうして普通の交流が得られ、これに電気時計の原理を応用して一定の速さで動く時計が作られる。これが水晶時計の原理である。

　こうして現在では、もはや地球の自転速度は一定ではなく、100年後の1日の長さは現在よりも0.0016秒長くなると考えられるに至った。原因は主として、ベーリング海を通る潮汐が地球と摩擦してその回転運動を止めることなどのためであるといわれている。

　このように自転が一定でないことが明らかになり、時間の単位の基準は地球の自転から公転へ移った方が妥当であるということになった。そして1956年の国際度量衡委員会で、時間の単位としての秒（s）を1900年1月0日12時（1899年12月31日正午、日本時間同日午後9時）の回転速度をもって太陽が春分点を通るまでの時間（太陽年）の1/3156925.9747と定義するに至った。我が国ではあ昭和33年4月15日以後この基準によっている。

　さらに時間の単位として原子的基準を使うことも考えられる。たとえばアンモニア分子は、 $^{86}$ Krが極めて純粋な一定波長の光を発すると同様に、一定波長の純粋なマイクロ波を発生する。このアンモニア蒸気に外から種々の波長のマイクロ波をあてると、この一定波長に等しいものだけは強い吸収を受ける。この現象を利用すれば、きわめて純粋な一定の波長のマイクロ波が得られる。もっとも純粋なのは22834.185メガサイクル・毎秒（Mc・S $^{- 1}$ ）のものである。水晶時計と同様な方法で周波数を下げ、普通の交流にしてから電気時計の原理で時計を作る。原子時計とよばれるのがこれである。

日本では統計 135° の地点すなわち明石市を通る子午線を基準に選んでいる。そこで、平均太陽は正午に明石市で南中するが、真太陽は必ずしもそうでない。

重力質量とその単位

　水は熱すれば水蒸気となり、冷やせば水となる。またロウソクは燃えて空気中の酸素と化合し、水蒸気や二酸化炭素に変わる。しかし姿は変えても物質は物質である。フラスコの中に閉じ込めた物質が、中でどんな変化をしようと、物質としての量は変わりがない。これが経験から得られた我々の観念で、こうして物質の量という概念が導かれる。

　しかし、この物質の量というのは観念的なもので、その内容がきわめて漠然としている。このままでは到底物理量として取り入れるわけにはいかない。これを明確な物理量として取り入れるためには、それが何で測られるかを明瞭にしなければならない。仮に、そのものの体積で表すという方法を採ったとしてみよう。シリンダーの中に閉じ込めた気体はピストンで圧縮することによってその体積は何ほどでも小さくなる。物質が減ったわけでもないのに体積は減る。このことは、明らかに体積が実質の量を示すのに適当でない物理量であることを示している。物質の量を表す物理量は、その物体に物質が出入りしない限り一定のものでなくてはならない。

　物体に付随する量はいろいろあるが、その中で重さという物理量はこの役目を果たすのにちょうど適当であることが分かる。すなわち、閉じたフラスコをばねはかりで測ってみると、中で何が起ころうと重さに変わりがない。一方、例えば純粋のように全く一様なものについていえば、二倍の体積のものの重さは二倍であって、このことも物質の量という概念と一致する。

　ただ、物体の重さは測定する場所によって変わる。したがって重さそのものは物質の量を表すのに適切でない。しかしこの場合でも重さを用いて物質の量を比較することはできる。なぜかといえば、経験の示すところによると、場所による物体の重さの増減はすべての物体に同じ割合で起こるからである。

　そこで物質の量を比較するには次のようにすればよい。

まず厳密に一様と考えられる物質（例えば白金塊）を選び、これを等しい体積に細分したとき、各片のもつ物質量は等しいと考え、これを基準物体に選ぶ。
他の一般の物質の場合では、その重さが同一の地点で測って上の基準物体の $n$ 個の重さに等しい時には、その物質量は基準物体の $n$ 倍と考える。

　このようにして決めた量が重力質量という物理量である。重力質量はこのように、物質の量を示すのに最も適すると考えられて作られたる量ではあるが、物質の量そのものではない。物質の量が漠然とした観念的なものであるのに対し、重力質量は明確な物理量である。物理量としての重力質量は、物質の量などという漠然とした観念から独立して考えるべき量で、むしろ物体に付随し、その属性を示す１種の量と考えるべきである。

　重力質量の基準となるものはキログラム原器である。1875年のメートル条約に基づいてメートル原器を作ったときと同時に、同様な合金（Pt 90％、Ir 10％）で作られた円筒形の物体で、その重力質量を1kgとしてある。この円筒はその直径と高さが等しい（共に39mm）ように作られているが、これは表面積を最小にして、外からの影響を最小にするための配慮に基づいている。

　この原器で最も問題になるのはその永久性で、ことにその表面に表ずる化学変化である。しかしこの原器に関する限り、化学反応その他で質量的に変化すると認められるような確かなものは現在のところ存在しない。鋳造当時この合金の内部に吸蔵された気体が徐々に遊離して出ていくために、原器が物質的に減少していくと考える人もあるが、これも疑いの程度を出ていない。こうして、キログラム原器は以前まで重力質量の標準となっていた。

単位の変換と次元

　長さ、質量、時間の単位には、m,kg,sの三つが用いられているが、m,kg,sもまた同様に用いられる。この他、米英などでは長さにfoot、質量にpoundなどが現在もなお用いられている。組立単位の大きさは基本単位の大きさによって変わる。例えば面積の単位の大きさは、単位長さの一辺をもつ正方形の面積と規定されているから、1mを基本としてできる面積単位は、1cmを基本とするものに対して10000倍の広がりをもっている。そこで長さの単位 $L_{1}, L_{2}$ に対する面積の単位を $S_{1}, S_{2}$ で表すと、 $\frac{S _{2}}{S _{1}} = (\frac{L _{2}}{L _{1}})^{2}$ cmとmの場合でいえば $\frac{S _{2}}{S _{1}} = (\frac{m}{cm})^{2} = 10 0^{2}$ である。今、長さ・質量・時間に対する２組の基本単位をそれぞれ $L_{1}, M_{1}, T_{1}, L_{2}, M_{2}, T_{2}$ としてみると、各種の組立単位の大きさの比はそれぞれ次のように与えられる：

$\frac{[ 体積 ] _{2}}{[ 体積 ] _{1}} = (\frac{L _{2}}{L _{1}}) 、 \frac{[ 密度 ] _{2}}{[ 密度 ] _{1}} = (\frac{M _{2}}{M _{1}}) \times (\frac{L _{2}}{L _{1}})^{- 3} 、 \frac{[ 速さ ] _{2}}{[ 速さ ] _{1}} = (\frac{L _{2}}{L _{1}}) \times (\frac{T _{2}}{T _{1}})^{- 1} 、$ $\frac{[ 加速度 ] _{2}}{[ 加速度 ] _{1}} = (\frac{L _{2}}{L _{1}}) \times (\frac{T _{2}}{T _{1}})^{- 2} 、 \frac{[ 運動量 ] _{2}}{[ 運動量 ] _{1}} = (\frac{M _{2}}{M _{1}}) \times (\frac{L _{2}}{L _{1}}) \times (\frac{T _{2}}{T _{1}})^{- 1}$

これらの式の右辺の指数を組立単位の次元（dimention）とよび、また速さの次元は $L T^{- 1}$ であるともいう。一般に、一つの物理量を単位 $U_{1}$ で表した時の数値を $n_{1}$ とし、単位を $U_{2}$ で表した時の数値を $n_{2}$ とすれば

$物理量 = n_{1} \times U_{1} = n_{2} \times U_{2} ∴ \frac{n _{1}}{n _{2}} = \frac{U _{2}}{U _{1}}$

　そこで、単位 $U$ の次元を長さ、質量、時間についてそれぞれ $x, y, z$ とすれば

$\frac{n _{1}}{n _{2}} = (\frac{L _{2}}{L _{1}})^{x} \times (\frac{M _{2}}{M _{1}})^{y} \times (\frac{T _{2}}{T _{1}})^{z}$

となる。いまこの式を仮に次のように書き直してみる：

$n_{1} \times L_{1}^{x} \cdot M_{1}^{y} \cdot T_{1}^{z} = n_{2} \times L_{2}^{x} \cdot M_{2}^{y} \cdot T_{2}^{z}$ 　この式は先ほどと違って、文字通り解釈したのでは意味がない。それは長さの単位を $x$ 乗するということが無意味であるし、さらに長さの単位と質量の単位の相乗積などということも意味をなさないからである。しかし形式的に $L, M, T$ を普通の数と同様に取り扱えばここまでに出てきた式と同じことにある。このことはまた $U_{1} = L_{1}^{x} \cdot M_{1}^{y} \cdot T_{1}^{z} 、 U_{2} = L_{2}^{x} \cdot M_{2}^{y} \cdot T_{2}^{z}$ と書くこともできる。これが単位記号として広く用いられ、次元式とよばれているもので、例えばMKS制における密度、運動量の単位はそれぞれ、kg・m $^{- 3}$ 、m・kg・s $^{- 1}$ と表される。いうまでもなくm・kg・s $^{- 1}$ などは単なる記号に過ぎないが、このように表しておく単位の換算にあたってkg、m、sを数と同様に考え、数の乗除法を適用できるという便宜が得られる。

[例題] 4℃における水の密度は1立方cmあたり1gである。これは1立方mあたり何kgか。

[解]まず密度の次元は $M \cdot L^{- 3}$ で表される。そこで、求める値を $ρ$ とすれば $密度 = 1 g \cdot c m^{- 3} = ρ kg \cdot m^{- 3}$ kg=1000g、m=100cmを代入すると $1 = ρ \times 1000 \times (\frac{1}{100})^{3} 、 ∴ ρ = 1000 、 ∴ 密度 = 1000 kg \cdot m^{- 3}$

量の間の関係式

　すでに知っているように、数量 $m$ kgの物体に $f$ Nの力が働いた時 $a$ m・s $^{- 2}$ の加速度が得られたとすると、 $m \times a = f$ の関係がある。この式に現れる $m, a, f$ はいずれも数値であって、この式は数式である。

　さて、この場合の力の単位は、単位質量に働いて単位加速度を生ずる力として導かれるものであるから、単位の換算にあたって前式からも分かるように $N = kg \cdot m \cdot s^{- 2}$ と書くことができる。そこで各数値の代わりにそれらが表す量の単位の次元式を代入すると、 $kg \times m \cdot s^{- 2} = kg \cdot m \cdot s^{- 2}$ となり、両辺の単位を次元式で計算したとき等しいことを示す。以上のことから、形式的に $m kg \times a m \cdot s^{- 2} = f kg \cdot m \cdot s^{- 2}$ $∴ （質量） \times （加速度） = 力$ という式が導かれる。この式の（質量）、（加速度）、（力）は物理量である。したがって、この乗法記号の×は単なる数の相乗記号でない。数値は数値としての計算で両辺の等しいことを示している。

　このことから、 $m, a, f$ を単なる数値と考えずに超そのものを表していると考えることもできる。そしてその方が便利なことが多い。なぜかといえば、数値は単位の選び方によって変わるが、量そのものは単位のいかんに無関係だからである。本書でもこの意味において $m, a, f$ などのような記号を量記号として用いることが多い。

古典力学

　昔々、まだ夜道を照らす灯もなく、人々の住まいは一日の半ば近くも真っ暗な闇の中に閉ざされ、その闇の奥から時折狼の遠吠えが、野原を渡る風の音に混じって気味悪く聞こえて来ていた頃、恐らく星は人々にとって今よりも遥かに身近で重要な存在だった。一日の仕事を終えて粗末な褥（しとね）に身を横たえる前のひとときや、東の空の白むのももどかしく眠りから目覚めた時、彼らが決まって振り仰ぐ星空には幾つもの親しみやすい図形が待ち受けていた。その或るものは美しい乙女が天に舞う姿を思わせたし、また或るものは大きな獅子となって頭上から彼らを威圧するかのように見えた。

　そのような図形や相互の配置のおおよそを、恐らく多くの人々はそらんじていたのであろう。実際、満足なあかりも、そしてまた、その下で読むべき何物をも持たなかった人々にとって、星空は彼らの涯しない夢を育むただ一つの、そしてそれ故にこそ何回となく読み古された物語のようなものだった。

　月日の流れを知る目安として、子は親から、親はまたその親から教えられたことであっただろうが、夜明け前の東の空や日没間もない西の空に輝く星の一団は、何故か季節によっていつも決まっているように見えた。天球上の太陽の通り道としての黄道（地球の公転軌道面と天球との公転）は、このようにしていつとはなしに人々の知るところとなった。

　上図のような黄道十二星座が、その通り道をたどる道しるべとして人々の間に根をおろしたのは、今からもう五千年以上も昔、バビロニア地方でのことといわれている。この黄道近く、これらの十二の星座を横切って、時には後戻りを交えながら、夜毎に少しづつその位置を変える幾つかの明るい星があった。こうして先ず、水星、金星、火星、木星、土星と呼ばれる比較的明るい5つの惑星が人々の前にその姿をあらわした。これらに大きな関心を示すようになったのは、占星術上の必要もあったであろうが、多分に宗教的意味も込められていたと考えられている。

　今日、我々は例えば上図に示された一連の写真のように、いくつかの惑星がたまたま天球上の一か所に集まるように見えたり、あるいは惑星の一つが黄道十二星座の中の一等星のそばを通過するように見えるときなどに、このような惑星の動きを実感をもって捉えることができる。だが極く限られた天文愛好家たちを除けば、今日われわれのほとんどは、たとえ夜空を眺めたとしても、金星以外の惑星の存在には気づかずに見落としてしまうのが普通である。だから、古代の人々が来る夜も星との対話の時間を十分に持ち合わせたということは、物理学の萌芽にとってこの上もなく大事なことだった。彼らが天球上の幾つかの星の不思議な動きに目をとめた時、人は後年古典力学の成立のきっかけとなった。この天球上の動きというのは、実は身近な物体の動きと同じ法則に基づいているため、以降では最初に身近な物体の運動について論じていくことにしよう。

運動の表現方法

　我々は間口と奥行きと高さを持った3次元空間の中に住んでおり、物体の運動もその中で起こっている。運動とは物体の位置が時間とともに変化することであるから、ある運動を記述するためには、例えば以下の図のように時々刻々における物体の位置を追いかけていけばよい。

物体の位置を記述するとき、必ずなんらかの基準から見たものとなる。例えばAさんとBさんの二人がいて、AさんがBさんに床に置かれた物体の位置を伝えようとすると「Bさんから見て右側」というようになることからも伺えるであろう。それを聞いたBさんも右側を見て、自分を基準として目分量で物体の位置がこれぐらいだと把握するであろう。このように、物体の位置を把握するためには基準点O（Bさん自身の地点）を中心として座標系K（Bさんの目分量）を用意する必要がある。さらに物体のある地点Pに質点を考えると、その位置は基準点O（または原点）から見た位置として記号 $OP$ で表記することができる。あるいは太文字で $r$ と表すと、 $r 、 ∣ r ∣$ は長さOPを示し、太文字で書くと長さに加えて向きまで含めたものとなる。このときの $r$ を位置ベクトルという（ベクトルの中でも長さしか持たない量はスカラー量と呼ばれる）。また、ここでは座標系KをBさんの目分量としていたが、それだと見る人によって大きさが異なってしまうため、一般的に下図のように直交直線座標（デカルト座標）が多く用いられる。

このときのP点の位置ベクトル $r$ は座標 $(x, y, z)$ を用いて

$r = (x, y, z)$

で与えられる。また、この図での $e_{x} 、 e_{y} 、 e_{z}$ は、それぞれ x、y、z 軸の正の方向に向かう単位の大きさのベクトル（基本ベクトル）

$e_{x} = (1, 0, 0) 、 e_{y} = (0, 1, 0) 、 e_{z} = (0, 0, 1)$

を表している。これを用いると、位置ベクトル $r$ は

$r = x e_{x} + y e_{y} + z e_{z}$

と書くことができ、後で出てくるが別の座標系に変更する際にどの座標系で表現されているか見やすくて便利なものとなっている。一方で、このときの原点から質点までの距離 $r$ は

$r = ∣ r ∣ = x^{2} + y^{2} + z^{2}$

で与えられる。

　ここまでに位置ベクトルについて述べたが、これは質点の運動において時刻 $t$ ごと変化するため $r (t)$ というように $t$ の関数として表すことができる。そのため、仮に時刻 $t$ で位置ベクトル $r (t)$ なるA点にあった質点が、時間 $Δ t$ の後に $r (t + Δ t)$ なるB点まで動いたとする。このときの質点の変位（displacement）Δr は、以下の図に図示されたようなベクトル算法を用いて

$Δ r = r (t + Δ t) - r (t)$

と書ける。また、変位 $Δ r$ の長さも

$Δ r = ∣ r (t + Δ t) - r (t) ∣$

となるが、同じ $Δ t$ という時間でも変位が大きければ大きいほど速いということになり、同じ $Δ r$ という変位でも時間が短ければ短いほど速くなるので、速さを表す度合い（速度）は $v$ として

$v \equiv \frac{r ( t + Δ t ) - r ( t )}{Δ t} = \frac{Δ r}{Δ t}$

と定義すればよいと考えられる（ $\equiv$ は定義することを指しており、右側にその内容を記載する）。この方向は変位 $Δ r$ と同じであり、細かく言うと平均速度と呼ばれる量となる。一方で、速度というのは一般的に各時刻ごとに異なる（極端に短い時間で変化する）量であるため

$v \equiv Δ t \to 0 lim \frac{r ( t + Δ t ) - r ( t )}{Δ t} = \frac{d r}{d t}$

というように記載すると厳密なものとなる（こちらは瞬間速度と呼ばれる）。この方向は先ほどと異なり、下図のようにA点において引かれた質点の軌道に対する接線の方向に等しく、その大きさは軌道に沿っての質点の速さ（speed）に等しくなる。

ちなみに、先ほどのデカルト座標を用いると速度は以下の通りに書ける。

$v \equiv (v_{x}, v_{y}, v_{z}) = (\frac{d x}{d t}, \frac{d y}{d t}, \frac{d z}{d t})$

　ここまで運動をどのように表現するかを述べたが、次からは実際に運動がどのような法則に従ってなされるのかを見ていくことにする。

運動の法則

　物体の運動を論じる古典力学は、Newtonが著した「プリンキピア」の導入部分に掲げた3つの法則に基づいている。ここでは、その3法則を当時の表現（特に第2法則の表現）に若干手を加えた形で紹介することにしよう。

慣性の法則
　すべての質点はそれに加えられた力によってその状態が変化させられない限り、静止あるいは1直線上の等速運動の状態をつづける。
運動方程式
　質点の運動量（=質量×速度）の変化は、加えられた力の方向に沿って起こり、微小時間内における運動量の単位時間あたりの変化の大きさは加えられた力の大きさに等しい。
作用・反作用の法則
　全ての作用に対して、等しくかつ反対向きの反作用が常に存在する。すなわち、互いに働きあう2つ質点の相互作用は常に等しく、かつ反対方向へと向かう。

　全ての物体は、他の物体の影響が完全にない状況において運動状態を変えることなく、その状態をいつまでも持続するであろう。多くの経験から予想されるこのような物体の性質を慣性（inertia）という。この性質があるということを第1法則で掲示した上で、外部からの影響が働いた時の運動の変化を第2、3の法則により述べている。これらを質点で言い換えると、質点の運動状態はその速さと運動方向（速度ベクトル）を指定することにより指定されるわけであるから、すべての質点は孤立状態において速度ベクトルを変えずに保とうとする性質を持つといえる。すなわち、運動する質点の速度ベクトルの変化を追うことで第2法則により質点に働く刻々の力を求めることはできる。ただ、それだけだと質点の運動状態の変化を力という言葉に置き換えただけに過ぎないため、力学現象をより深く理解することには何等寄与もしていない。そのため、色々な運動の度重なる観測に基づいて、どのような種類の力がどのような状況のもとで働くのかを予め調べる必要が出てくる。Newtonが発見した万有引力の法則は、その一つの例であるといえる。そして、現在ではそれらの力（既知の力）のもとに質点の運動を第2法則から予測することが可能となっている。

　では、いよいよ第２法則の解析的表現に取り組むことにしよう。まず、速度 $v$ で動いている質量 $m$ なる質点の運動量を $p \equiv m v$ によって定義する。この量は質点に力 $F$ が働くとき変化するが、両者の間には $\frac{d p}{d t} = F$

力を用いた運動の記述

水平な地面に置かれた物体

何もしていない状態だと静止し続ける。そこに下図のように力をかけ続けると動き始め、この様子を運動方程式を用いてみる。 $m a = F \leftrightarrow {m a_{x} = F_{x} = F m a_{y} = F_{y} = 0$ このことから、 $x$ 方向に対しては $F / m$ の加速度で進むが、 $y$ 方向へは一定の速度（今は静止したまま）で動くことが分かる。だが、すぐ動き出すわけではなく最初はある程度静止し続けてから動き出す。これは、摩擦力 $W$ が働いているものと考える。 $m a = F \leftrightarrow {m a_{x} = F_{x} = F - W m a_{y} = F_{y} = 0$
上記と同様な状況で角度 $θ$ だけ傾けて力をかけた場合

先ほどと同様な考え方で $m a = F \leftrightarrow {m a_{x} = F_{x} = F cos θ - W m a_{y} = F_{y} = F sin θ$ となるわけだが、もしこの力だけ存在するならば $x, y$ 方向へ加速していくものと考えられる。しかし、実際はそうはならず $x$ 方向にのみ動く。そのため、下向きにも力が働くことが予想される。ただ、この場合だと地面が邪魔しているため、次に持ち上げた状態で落とした場合にどのような力がはたくか見てみる。
落下運動と2と同様の運動

物体を持ち上げて離すと、落下するわけだが、このときに働く力は重力 $G$ と呼ばれる。そのため、このときの運動方程式は以下の形になることが分かる。 $m a = F \leftrightarrow {m a_{x} = F_{x} = 0 m a_{y} = F_{y} = - G$ この重力が働くということを踏まえて、改めて先ほどの状況を考えてみる。この場合だと、重力と上にひく力が釣り合っていると思えるが、すると何も力が働いていないときに静止していることの説明がつかない。そのため、重力により下へ引き付けることで地面へ力がかかり地面から反発力（反作用）が働くものと考えられる。

$m a = F \leftrightarrow {m a_{x} = F_{x} = 0 m a_{y} = F_{y} = N - G = 0$ そのため、先ほど角度をつけた場合にも以下のようになる。 $m a = F \leftrightarrow {m a_{x} = F_{x} = F cos θ - W m a_{y} = F_{y} = F sin θ + N - G = 0$

以上のように、運動を観測することでどのような力が働いているのか予測することが可能ではあるが、実際にどのような値をとるかは調べる必要が出てくる。例えば、重力 $G$ の場合だと $m g$ という形になり、摩擦力 $W$ の場合は静止している場合には $μ N$ 、動き出した後には $μ^{'} N$ になることが分かっている（ただ、この形というのは地面の表面の状態などで変化するもので実際にこのようになる理由はいまだにわかっていない）。

万有引力の法則

　質量を持つ二体間には万有引力という力が働くことが知られており、例として物体1と物体2の質量、位置ベクトルを $m_{1}, r_{1}, m_{2}, r_{2}$ とすると、物体1にかかる力は $F = G \frac{m _{1} m _{2}}{∣ r _{2} - r _{1} ∣ ^{2}} (\frac{r _{2} - r _{1}}{∣ r _{2} - r _{1} ∣})$ というように書ける（ $G$ ：重力定数）。そして、この２体の運動方程式を記載すると以下の通りとなる。 $m_{1} \frac{d ^{2} r _{1}}{d t ^{2}} = F 、 m_{2} \frac{d ^{2} r _{2}}{d t ^{2}} = - F$ ここで二つの式の和をとった場合には $(m_{1} + m_{2}) \frac{d ^{2}}{d t ^{2}} (\frac{m _{1} r _{1} + m _{2} r _{2}}{m _{1} + m _{2}}) = 0$ というように全体の運動方程式となり、このときの質量 $m_{G}$ と位置 $r_{G}$ （重心）は $m_{G} = m_{1} + m_{2} 、 r_{G} = \frac{m _{1} r _{1} + m _{2} r _{2}}{m _{1} + m _{2}}$ とおけることから、全体的に見た場合には以下の式が成り立つ。 $m_{G} \frac{d ^{2} r _{G}}{d t ^{2}} = 0$

一方で、万有引力の形を踏まえて2つの方程式を $\frac{d ^{2} r _{1}}{d t ^{2}} = \frac{F}{m _{1}} 、 \frac{d ^{2} r _{2}}{d t ^{2}} = - \frac{F}{m _{2}}$ として2番目の式から1番目の式を差し引くと以下の式が得られる。 $(\frac{m _{1} m _{2}}{m _{1} + m _{2}}) \frac{d ^{2}}{d t ^{2}} (r_{2} - r_{1}) = - F$ この式は物体1から見た物体2の運動方程式となっており、質量 $μ$ （換算質量）と相対位置 $r$ は $μ = \frac{m _{1} m _{2}}{m _{1} + m _{2}} 、 r = r_{2} - r_{1}$ と書ける。そのため、上記の運動方程式は以下のようになる。 $μ \frac{d ^{2} r}{d t ^{2}} = - G \frac{m _{1} m _{2}}{r ^{3}} r$ ここで力が原点からの距離にしか依存していないことから、ベクトルを動径成分で記載したほうが便利なように思える。 $r = r e_{r} 、 \frac{d r}{d t} = \frac{d r}{d t} e_{r} + r \frac{d e _{r}}{d t} 、 \frac{d ^{2} r}{d t ^{2}} = \frac{d ^{2} r}{d t ^{2}} e_{r} + 2 \frac{d r}{d t} \frac{d e _{r}}{d t} + r \frac{d ^{2} e _{r}}{d t ^{2}}$ ここで、右辺に対して単位ベクトルを成分で書くと以下の通りにかける。 $e_{r} = (cos θ, sin θ) 、 \frac{d e _{r}}{d t} = \frac{d θ}{d t} (- sin θ, cos θ) = \frac{d θ}{d t} e_{θ}$ $\frac{d ^{2} e _{r}}{d t ^{2}} = \frac{d ^{2} θ}{d t ^{2}} e_{θ} + \frac{d θ}{d t} \frac{d e _{θ}}{d t} = \frac{d ^{2} θ}{d t ^{2}} e_{θ} - (\frac{d θ}{d t})^{2} (cos θ, sin θ) = \frac{d ^{2} θ}{d t ^{2}} e_{θ} - (\frac{d θ}{d t})^{2} e_{r}$ そのため、改めて置き換えて表すと以下の通りになる。 $r = r e_{r} 、 \frac{d r}{d t} = \frac{d r}{d t} e_{r} + r \frac{d θ}{d t} e_{θ} 、 \frac{d ^{2} r}{d t ^{2}} = [\frac{d ^{2} r}{d t ^{2}} - r (\frac{d θ}{d t})^{2}] e_{r} + [r \frac{d ^{2} θ}{d t ^{2}} + 2 \frac{d r}{d t} \frac{d θ}{d t}] e_{θ}$ これらの式と運動方程式を比較すると、 $μ \frac{d ^{2} r}{d t ^{2}} - μ r (\frac{d θ}{d t})^{2} = - G \frac{m _{1} m _{2}}{r ^{2}} 、 μ r \frac{d ^{2} θ}{d t ^{2}} + 2 μ \frac{d r}{d t} \frac{d θ}{d t} = 0$ ここで２番目の式に $r$ をかけることで $μ r^{2} \frac{d ^{2} θ}{d t ^{2}} + 2 μ r \frac{d r}{d t} \frac{d θ}{d t} = \frac{d L}{d t} = 0 、 L = μ r^{2} \frac{d θ}{d t}$ と変形できることから、この量（角運動量）が保存されている（ $L$ が定数になる）ことが分かる。そこで、これを変形して $\frac{d θ}{d t} = \frac{L}{μ r ^{2}}$ として、一つ目の式に代入してみる。 $μ \frac{d ^{2} r}{d t ^{2}} = \frac{L ^{2}}{μ r ^{3}} - G \frac{m _{1} m _{2}}{r ^{2}}$ 次に、一つ目の式の両辺に $r^{2}$ を乗じて $u = 1/ r$ と置き換えてみると $μ \frac{1}{u ^{2}} \frac{d ^{2} r}{d t ^{2}} = \frac{L ^{2}}{μ} u - G m_{1} m_{2} 、 \frac{d θ}{d t} = \frac{L}{μ} u^{2}$ となるが、ここで $r$ の時間微分がどうなるか見てみると $\frac{d r}{d t} = - \frac{1}{u ^{2}} \frac{d u}{d t} = - \frac{L}{μ} \frac{d t}{d θ} \frac{d u}{d t} = - \frac{L}{μ} \frac{d u}{d θ}$ であることから、二階微分は次の通りとなる。 $\frac{d ^{2} r}{d t ^{2}} = - \frac{L}{μ} \frac{d}{d t} (\frac{d u}{d θ}) = - \frac{L}{μ} \frac{d θ}{d t} \frac{d ^{2} u}{d θ ^{2}} = - \frac{L ^{2}}{μ ^{2}} u^{2} \frac{d ^{2} u}{d θ ^{2}}$ そのため、運動方程式は以下の形に整理される。 $\frac{d ^{2} u}{d θ ^{2}} = - u + \frac{G m _{1} m _{2} μ}{L ^{2}}$ この解は以下の通りとなる。 $u = A cos (θ - θ_{0}) + \frac{G m _{1} m _{2} μ}{L ^{2}}$ そのため、 $r$ は以下のように表される（ $ϵ$ ：離心率、 $l$ ：半直弦）。 $r = \frac{l}{1 + ϵ cos ( θ - θ _{0} )} 、 ϵ = \frac{A L ^{2}}{G m _{1} m _{2} μ} 、 l = \frac{L ^{2}}{G m _{1} m _{2} μ}$ この軌道は離心率 $ϵ$ の大きさによって決まる。 $（ 1 ） ϵ = 0 （円）、（ 2 ） 0 < ϵ < 1 （楕円）、（ 3 ） ϵ = 1 （放物線）、（ 4 ） ϵ > 1 （双曲線）$

　一方で、惑星の運動に関してKeplerの法則というものがあり、以下の3つの法則が存在している。

惑星の軌道は太陽を1焦点とする楕円である。
太陽と惑星とを結ぶ動径の描く面積速度は常に一定である。
惑星の公転周期の2乗は太陽からの平均距離（=軌道の長半径）の3乗に比例する。

先ほどの物体について、物体1を太陽、物体2を地球として考えると、離心率は0.0167になり範囲内にあることはわかる。次に、面積速度については以下のように動径によらずに一定になることが分かる。 $\frac{d S}{d t} = \frac{1}{2} r^{2} \frac{d θ}{d t} = \frac{L}{2 μ}$ そして、公転周期については楕円の面積が長半径が $a$ 、短半径 $b$ により $S = πab$ と書けるが、ここでの長半径、短半径が $r$ が最大と最小のときの長さの相加平均、相乗平均であることから $a = \frac{1}{2} (\frac{l}{1 - ϵ} + \frac{l}{1 + ϵ}) = \frac{l}{1 - ϵ ^{2}} 、 b = (\frac{l}{1 - ϵ}) (\frac{l}{1 + ϵ}) = \frac{l}{1 - ϵ ^{2}}$ というように書けるため、面積は以下の通りとなる。 $S = \frac{l ^{2} π}{( 1 - ϵ ^{2} ) ^{3/2}}$ 一方で、面積速度から面積を求めた際には $S = \int_{0}^{T} \frac{d S}{d t} d t = \frac{1}{2} \int_{0}^{T} r^{2} (\frac{d θ}{d t}) d t = \frac{L T}{2 μ}$ であるため、以下の式が成り立つ。 $\frac{l ^{2} π}{( 1 - ϵ ^{2} ) ^{3/2}} = \frac{L T}{2 μ} \to \frac{l ^{4} π ^{2}}{( 1 - ϵ ^{2} ) ^{3}} = \frac{L ^{2} T ^{2}}{4 μ ^{2}}$ ここで、 $L^{2} = G m_{1} m_{2} μ l 、 a = l / (1 - ϵ^{2})$ であるから、 $\frac{l ^{3} π ^{2}}{( 1 - ϵ ^{2} ) ^{3}} = \frac{G m _{1} m _{2} T ^{2}}{4 μ} \to \frac{a ^{3}}{T ^{2}} = \frac{G m _{1} m _{2}}{4 π ^{2} μ}$ となる。ただ、太陽に比べて地球の質量が軽いものとできるため近似的に以下の式で計算ができる。 $\frac{a ^{3}}{T ^{2}} ≃ \frac{G m _{1}}{4 π ^{2}} 、（ μ ≃ m_{2} ）$

質点系の力学

　ここまでは、単に物体を質点におきかえて正確に位置を指定してから、その時間変化を見ることで物体の運動を調べていた。このとき、質点を置く地点は物体内であれば何処でもよく特に指定があるわけではないが、場合によっては問題となってしまうことがある。

　例えば、上図のように質点が矢印に沿って運動を行ったものとすると、左側（並進運動）では他の位置に質点を指定しても矢印と同じ運動を行うので、全体的に物体は並進運動していることになる。しかし、右側（回転運動）にあるように途中で引っかかり回転した場合には質点の位置によって運動が異なるため、単純に質点として扱えないことが分かる、同様に、変形する場合でも場所ごとで運動の様子が違うことがわかる。このような運動はどのようにすればよいのか考えてみると、視点を変えてそもそも質点の指定する位置ごとで運動が違っているものと予想できる。言うなれば、それぞれの地点での質点は異なる運動方程式に従っており、全体（物体）の運動を表す際はそれら全て見ていく必要があることになる。そこで、それぞれの地点の質点が別の運動方程式に従うものとして、それらを踏まえて物体全体の運動を記述するという方法をとる。

　今、物体を質点の集合体としてその内 $i$ 番目の運動方程式を記述すると $m_{i} \frac{d ^{2} r _{i}}{d t ^{2}} = F_{i} + j = 1 (j \neq = i) \sum N F_{ij}$ と書ける。ここで、 $m_{i} 、 r_{i} 、 F_{i} 、 F_{ij} 、$ はそれぞれ質点の位置ベクトル、位置ベクトル、外部からの力（外力）、別の $j$ 番目の質点との間に働く力（内力）を表している。さて、内力に関して作用・反作用の法則 $F_{ij} = - F_{ji}$ が成り立つため、先ほどの式を $i$ に対して総和をとると以下のように内力の項が相殺されて0となる。 $i = 1 \sum N (m_{i} \frac{d ^{2} r _{i}}{d t ^{2}}) = i = 1 \sum N F_{i}$ ここで、右辺が質点全体にかかる力（合力）であることから全体的に見て $m \frac{d ^{2} r}{d t ^{2}} = F$ という形になっていると考えられる。そこで、 $m$ を質点全体（物体自体）の質量として $i = 1 \sum N (m_{i} \frac{d ^{2} r _{i}}{d t ^{2}}) = m \frac{d ^{2}}{d t ^{2}} (\frac{1}{m} i = 1 \sum N m_{i} r_{i})$ というように変形して比較してみると、 $r$ は以下の形で置くことができる。 $r = \frac{1}{m} i = 1 \sum N m_{i} r_{i}$ この位置ベクトルが指す地点は、例として質点が2個（ $N = 2$ ）の場合だと $r = \frac{m _{1} r _{1} + m _{2} r _{2}}{m _{1} + m _{2}}$ となっているため、内分点の公式と同じ形になっていることが見て取れる。

そのため、

r

は質点1と2の間を内分する地点の位置を指しており、この位置を物体の位置として扱うことになる。言い換えれば、物体を地点

r

にある質量

m (= m_{1} + m_{2})

の質点に置き換えるという操作を行っており、従来のときとほぼ同様なことをしていることが分かる。ただし、その位置はどこでもよいわけではなく、上式で制限されているところが従来と異なる。同じように、質点が3個の場合でも質量

m

の質点の位置は以下のようになる。

r = \frac{m _{1} r _{1} + m _{2} r _{2} + m _{3} r _{3}}{m _{1} + m _{2} + m _{3}} = \frac{( m _{1} + m _{2} ) r _{12} + m _{3} r _{3}}{( m _{1} + m _{2} ) + m _{3}} 、 r_{12} = \frac{m _{1} r _{1} + m _{2} r _{2}}{m _{1} + m _{2}}

　以上が物体を質点系とみなしたときの話であるが、これを見ても分かるように、以前と同じように物体を位置 $r$ にある質量 $m$ の質点として扱うことができ、並進運動にも通じる考え方であると言える。そのため、ここからは物体を質点系として議論を進めていくことにする。ただ、その際に質点の質量 $m_{i}$ を物体内での質量分布が一様であるものとして $m_{i} = ρ = \frac{m}{V}$ というように、密度 $ρ$ （ $m$ ：物体の質量、 $V$ ：物体の体積）でおいておく（このようにすることで、運動のほうに重点を置いて考察を進めることができる）。　

熱力学

　現在のように熱力学と呼ぶようになったのは19世紀に入ってからであり、それ以前は単に熱学であった。熱学の研究の起源については、他の科学の歴史と同様に、古代ギリシャ、ローマにまで遡らなければならない。当時、自然現象を研究する基本的方法は、見たり感じたりなどする簡単な観察が行われた程度であった。すでに当時でも、最も簡単な熱作用の現れである蒸発、沸騰、融解などが知られていたが、これらの現象を説明しようとする試みは成功していなかった。

　人間がこの分野における知識をはじめて充実させ、体系化して物理的解釈を与えたのはやっと17世紀になってからのことである。それは熱現象を単に観察するだけではなく、実験的に研究しようとするようになったからである。この頃になると温度計の発明がなされていた。たとえば、ガリレイは空気の熱膨張を利用した一種の気体を考案している。あとで述べるように、引き続いていろいろな温度計が作られた。こうして温度が定量的に測定できるようになると、それに応じて熱学も発展していった。特に18世紀に混乱をまねいた熱理論の2つの基本概念である温度と熱量を正確に区別することが可能となった。

　実験に基づいた熱現象の研究は、いわゆる熱と呼ばれるものの性質についての問題解決へと近づいて行った。しかしながら、その本質が何であるかについては、18世紀の終わりごろまでは熱を運動の一形態と考える学者と、熱を一種の物質とする熱素説（カロリック説）をとる学者とがあった。熱素説は18世紀半ばすぎに起こった。提唱者はフランスのラヴォアジエで、1787年の命名である。熱という一種の物質があるとする考え方はギリシャ時代からもあったが、それが18世紀に復活した背景には、光と同様にニュートン物理学が有力に働いていたことも事実であろう、18世紀のはじめにニュートンは光の粒子説を唱えていた。したがって、熱現象を熱素によって説明することは、当時多くの学者たちにとって全く当然に思われたのである。のみならず、カロリックによって純熱学的現象を説明することにある程度成功したのである。例えば、イギリスのブラックは熱素説の立場から熱容量や潜熱の概念を明らかにした。一方で同じ頃、独立に比熱、潜熱、熱膨張などについて研究していたイギリスのキャベンディッシュは、ブラックの潜熱という表現に反対して、水蒸気が凝結するときは”熱が発生する”という表現をとった。彼は熱が物体中の粒子の内部運動であるという考えを基にして、今日の気体分子運動論的立場をとっていた。そのほか、ロシアのロモノソフは熱の本性が物体を構成する微粒子の回転運動によるものと考えていた。このように、物体を構成する微粒子の運動が熱の本質であるという考えをもとに、逆に完全に運動が静止するような温度が存在しなければならないという結論も下している。これは後にトムソンによって発見された温度（絶対零度）であったのだが、熱素説が注目されていた当時としては影が薄いものとなっていた。

　しかし、熱学的現象に力学的現象が加わるような現象を説明するようになってから熱素説に陰りが出始めた。その一つとしてランフォードの実験がある。すなわち、彼がミュンヘンの造兵廠（ぞうへいしょう）で大砲作りの監督をしていたとき、砲身を穿孔（せんこう）する際に熱が無尽蔵に出てくることから考えて、1798年に熱が物質ではなく運動の中にあると結論を下したのである。1799年イギリスのデービィも空気ポンプを使って、真空中で2つの氷片を時計仕掛けで互いにすり合わせてできる水を観測し、運動と熱との関係に注目した。これらの実験は、これから先いろいろな形でほかの学者によって繰り返され、どれも熱素説に反することを証明したものとなっていた。このことから、熱というものに関する基本的な概念を根本的に再検討することの発端になった。そして、熱が運動の特殊な形態であることが熱の仕事当量の発見により確認されたのである。この仕事当量というのは、熱と仕事が互いに変換されるときの量で定数（仕事の状態によらない決まった値）をとっている。ドイツの医学者マイヤーが初めて自身の人体生理学の問題において取り上げ、その後イギリスのジュールとマイヤーが独立して正確な値を測定している。ここで、熱が仕事に置き換わることからエネルギーの一形態であることが言えるので、クラウジウス、トムソンによってエネルギー保存則の原理が打ち立てられ、1847年ドイツのへルムホルツが独立にNewtonの運動法則からエネルギー保存則を発表している。Newtonの運動法則とあるように、力学の知識が取り入られたことで熱学は熱力学となり、熱と運動との関係により熱現象を論じる学問として発展していった。

温度

　温度という語は、われわれの受ける温かい、冷たいという感覚から生じたもので、温冷の度合いを表す物理的な尺度と言える。これは力の概念が筋肉の緊張から生じたものに対比される。日常生活で火が熱く、氷が冷たいというのは誰にも共通の感じであろう。また、井戸の水が冬温かく、熱冷たいというのも同じであろう。しかし、火が氷より温度が高いというのは本当であっても、井戸の水の温度が冬高いというのは早計である。このように人間の感覚は温度の標準としては適当ではない。そこで温度が物理的な尺度となるためには、人間の感覚を離れて客観的に定義されることが必要である。

　実験の結果によると、2つの物体を外界と隔離して接触させておくと、いつも熱い方は冷え冷たい方は温まっていき、ついにはそれ以上変化が起こらなくなってしまう。このとき、この2物体は熱平衡に達したといい、このような状態に達した後では力学的つり合いを破らない（圧力や体積などを変化させない）ようにして２つの物体を離してもそれぞれの物体は熱平衡の状態にある。熱平衡の状態では温度は変化せず一定となることは経験的に理解しているところである。

　いま、物体AとBが接触して熱平衡にあり、AとCが接触して熱平衡にあれば、BとCを接触させてもやはり熱平衡である。この事実を熱力学第0法則ということがある。ここで、互いに熱平衡にある物体について等しい値をとるような量として温度を考えると、AとBが等しい温度でAとCが等しい温度であればBとCはまた等しい温度といえる。ゆえに、2つの物体BとCが等しい温度であることを確かめるのに、BとCを直接接触させなくても間接的にAという物体を別々にBとCを接触させて調べることができる。このAの役割を担うものが温度計である。こうして定められた温度は後述の理論温度に対して経験温度という。

　温度計としては何を選んでもよいが、容易に再現できるものを選ぶのが便利である。ただ、どの温度計も目盛りの示度が同じでなければ意味がない。この正しい温度目盛りを見出すためには、これまで述べた熱平衡の法則以外を調べなければならない。、歴史的には

エントロピーを次のように定義する。 $d S \equiv \frac{d ^{'} Q}{T}$

温度が異なるAとBの系（ $T_{A} > T_{B}$ ）を接触させた場合

$d S_{A} \equiv \frac{d ^{'} Q _{A}}{T _{A}}$

$d S_{B} \equiv \frac{d ^{'} Q _{B}}{T _{B}}$

となるが、ここで温度が熱の変化によって変化するため例えばAの場合だと

$\frac{d ^{'} Q _{A}}{T _{A}} = \frac{d ^{'} Q _{A}}{T _{A （一定）}}$ ではなく、 $d^{'} Q_{A} < 0 、 T_{A （一定）} \geq T_{A}$ であるため $\frac{d ^{'} Q _{A}}{T _{A （一定）}} \geq \frac{d ^{'} Q _{A}}{T _{A}}$ となり、以下の不等式が成り立つことが分かる。 $d S_{A} \geq \frac{d ^{'} Q _{A}}{T _{A}}$ Bについても同様な議論を行うことで、 $d^{'} Q_{B} > 0 、 T_{B} ≦ T_{A}$ であることをふまえて $d S_{B} \geq \frac{d ^{'} Q _{B}}{T _{B}}$ になる。エントロピー全体としても示量量であるため、 $d^{'} Q > 0 、 d^{'} Q_{B} = d^{'} Q 、 d^{'} Q_{A} = - d^{'} Q$ $d S = d S_{A} + d S_{B} = (\frac{1}{T _{B}} - \frac{1}{T _{A}}) d^{'} Q \geq 0$ これから孤立系では熱平衡出ない状態であればエントロピーが増大し続け、熱平衡になればエントロピーがある値をとるようになる。この値はBoltzmannの原理により $S = k_{B} ln W$ で求めることができる。

stefan-boltzmannの法則もここで説明して

電磁気学

　電磁気学が難しいとされる理由はいろいろあげられるが、その第一は力学と違って、そこで扱われる物理的対象がわれわれのに日常的な感覚を超えるものであることだろう。電荷や電流そのものも目に見えない。それでもそれらは、静電気や火花やコンセントのショートなどの体験によって、その存在を感知することができる。しかし、電磁気学の主役である電場や磁場となると、それらはまったく人間の五感とは無縁の存在であり、そのためにそれらは初学者にあたかも幻との格闘を強いている趣さえ感じられる。第２の理由は、力学と違って学ぶにつれて種々雑多の法則が次々と現れてきて、それらの法則の上下関係や相互関係が分からなくなり、全体の話の筋道が見失われがちな点にある。第3のつまずきの石は、そこに用いられている数学にあるようである。電磁気学で使用される数学は、主としてベクトル場の微分・積分という数学的にもかなり高度なものである。そのため、その数学に気をとられて物理どころではなくなってしまうのである。しかも、物理学と並行して学習しているはずの数学の内容が、物理学のそれに追いつかないことも原因の一つにあげられよう。

　そこで、上記のような困難を克服するために、まずはなぜ電場や磁場というような抽象的な外見を導入する

静電気

　コハクをこするとチリを引き付けることは、既に紀元前600年もの昔から知られていた。この種の現象は、冬の日にセーターを脱ぐときに起きる火花などとともに日常よく経験することである。これらの現象を詳しく調べるため、ガラス棒を絹布でこするか、あるいは樹脂棒を乾燥した毛皮でこすってみる。すると、これらの棒は紙片などをよく引き付ける。このときこれらの棒は、目に見えないが、何かを帯びたと考えられる。この何かを電気または電荷とよび、電気を帯びることを帯電するという。

電気の性質を調べるには、電気振り子を用いると便利である。それは右図に示すように、軽いコルク球を絹糸でガラスの支柱につるしたものである。絹でこすったガラス棒をコルク球に近づけると、はじめはコルク球は棒に引き付けられるが、それがガラス棒に触れた瞬間にコルク球は勢いよく反発して、ガラス棒から遠ざかろうとする。これと同じことをガラス棒の代わりに樹脂棒で試しても同様な現象を示す。

ところが、ガラス棒と樹脂棒に帯電している電荷が同じものかというとそうではない。実際、ガラス棒に触れたコルク球に帯電した樹脂棒を近づけると、ガラス棒とは反対にコルク球が引き付けられることが確認できる。さらに、そのコルク球に樹脂棒が触れると、それら2つの電気的性質は失われる。このことからも、2つの棒に帯電した電荷は異なることが分かるため、ガラス棒に帯びた電荷を正、樹脂棒の帯びた電荷を負と定義してみると、上記の現象を説明できる（同符号の電荷間には斥力、異符号の電荷間には斥力、また正負の電荷が代数的な数のように加え合わさった和が0となると電気的性質が失われると考える）。

　以上のように、正負の電荷を定義することで一種の電気的性質を説明することができるわけだが、ここで電荷間で働く力はどれくらいなのかというと1785年にColmonbが行った実験により確認されている。具体的には、以下の図に示すような装置を作り、点電荷間に作用する力の性質を調べた結果、次の法則が成立することを発見した。これをクーロンの法則と呼び、この力をクーロン力という。

2個の点電荷間の力は、両者を結ぶ直線上に働き、その大きさはそれぞれの電荷の量に比例し、またそれらの間の距離の2乗に反比例する。

Coulombが行った実験の装置
水晶の糸で吊るした絶縁体の棒の一端の小金属球Aに電荷を与え、また別に固定した絶縁体の下端の小金属球Bにも電荷を与えて、その間の力の大きさを糸のねじれとすり合わせることにより測った。このときの両球間の距離の決め方にあいまいさが生じ、それによる実験誤差が約10％にもなってしまう。しかし、現在ではほかの方法を用いたより精密な実験によって、上のCoulombの法則がきわめて正確に成立していることが確かめられている。

さて、二つの点電荷の電荷量をそれぞれ $q$ および $Q$ とし、それらの間の距離を $r$ とすると、その間に働く力の大きさ $F$ は $F = k \frac{qQ}{r ^{2}}$ と表わされる。ここで $k$ は正の比例定数であり、 $qQ > 0$ なら $F > 0$ で斥力を示し、 $qQ < 0$ なら $F < 0$ で引力を示す。比例定数 $k$ を決めるには、電荷量、距離および力の単位を決めておく必要がある。現在ではSI単位系の下で単位が決められており、1Cの電荷を持つ2個の点電荷を真空中で1mの距離を置いて静止させたとき、これらの間に働く力が $8.9876 \times 1 0^{9} N$ であるから $k = 8.9876 \times 1 0^{9} N ・ m^{2} / C^{2}$ となる（実際に測るのはクーロン力ではなく、真空中の光の速さであるが、この点については後述）。

　ここで、クーロン力の大きさがどれくらいのものかを知るために、同じ形で表される万有引力の法則 $F_{G} = G \frac{m M}{r ^{2}}$ で与えられる重力とクーロン力を比較してみよう。ここで $m$ と $M$ は2個の質点の重量で、 $r$ はその間の距離、また $G = 6.7 \times 1 0^{- 11} N \cdot m^{2} /k g^{2}$ は万有引力定数である。2個の質点として電子を考えると、電子の質量 $m = 9, 1 \times 1 0^{- 21} kg$ 、その電荷量 $e = - 1.6 \times 1 0^{- 19} C$ を用いて、二つの力の大きさの比をとると $\frac{F}{F _{G}} = \frac{k e ^{2}}{G m ^{2}} \approx 4.2 \times 1 0^{42}$ という、とてつもない巨大な値になる。このように、クーロン力は重力と比較して圧倒的に大きく、原子や分子の世界では重力の効果は無視できる。それにも関わらず地球上の物体あるいは天体間に重力が作用する理由は、原子や分子のほとんどすべてが電気的に中性であって、電気的な引力と斥力とが相殺してしまうからである。

　クーロン力というのは、その方向をも指定して決まるベクトル量である。そのため、以下の図において電荷 $Q$ が電荷 $q$ に作用する力をベクトルを用いて表すと次のようになる。

F (r_{q}) = k \frac{qQ}{∣ r _{q} - r _{Q} ∣ ^{2}} \frac{( r _{q} - r _{Q} )}{∣ r _{q} - r _{Q} ∣}

ここで、ベクトル $r_{q} = (x_{q}, y_{q}, z_{q})$ と $r_{Q} = (x_{Q}, y_{Q}, z_{Q})$ は、それぞれ図の原点Oから測った電荷 $q$ と $Q$ の位置ベクトルであり、 $F (r_{q})$ は電荷 $q$ に作用する力を表す。

遠隔作用と近接作用

　Newton力学において、万有引力というのは物体間にその力を伝える物体（糸など）があるわけでなく、空間を通して直接物体間に作用する力であると考えられる。このような力を遠隔作用の力という。電気現象でも、正負に帯電した粒子（帯電体）間にはCoulombの引力が働く。しかも、これらの力は全く同じ形をしているため、どちらも同じ遠隔作用の力であると考えることができる。しかし、帯電体の途中にある空気などのものに影響を与えず別の帯電体に力を及ぼすというのは如何にも不思議に思われるため、別の考えがFaradayによってなされた。彼は以下のような例を考えて、その力について考察を行った。

今、水を容器に入れて、その上に水をはじくようロウを塗ったコルク球を2個浮かべてみる。すると、コルク球の周りの水面はへこんでゆがんだ状態になる。この2つのコルク球を数mm程度まで近づけて手を離すと、コルク球は互いに引き寄せられ、そして接触する。これは、水がその表面張力によって表面積をできるだけ小さくしようとするためである。今、仮に人間にこの水が見えないものとすると、我々には2個のコルク球間に引力が作用しているように感じられるだろう。このようにコルク球がその周りの水面をゆがめ、その歪みが次々と水面上を伝わって、ほかのコルク球に力を与えるとする考え方を近接作用の考え方という。

　以上の考え方を踏まえると、帯電体があると周囲の空間はそれ自身の性質によりゆがむと考えられる（先ほどで言う水がコルク球でへこんだ状態になったことに相当する）。この歪みを電場（電気工学系では電界）という。結局、電磁気学とはこの空間内に実在する電場（と磁場）の性質を研究する学問である。そして、電場の中に別の帯電体を置くことでクーロン力が生じるわけだが、ここで $F (r_{q}) = q E_{Q} (r = r_{q})$ $E_{Q} (r) = \frac{k Q}{∣ r - r _{Q} ∣ ^{2}} \frac{( r - r _{Q} )}{∣ r - r _{Q} ∣}$ と表わすと次のように考えることができる。点電荷 $q$ のあるなしに関わらず、位置 $r_{Q}$ にある点電荷 $Q$ は空間のあらゆる場所 $r$ に電場 $E$ を作り、 $r = r_{q}$ の位置にもう一つの点電荷 $q$ があると、それに $F$ の力が作用する。

　電場がどのように発生するかを見るのに電気力線がよく用いられる。これは以下の図に示すように。その接線方向が電場の方向と一致するように線を描いたもので、一つの電荷の場合だと始点か終点のどちらかが無限遠に伸びた形になる。

これが2つになると、二つの電荷の電場を足したものになり、電荷の符号が同じか異なるかで電気力線の形が違ってくる（同じ電荷量であれば同じ本数だけ放出、吸収する、符号が反転すれば向きが逆になる）。

$E = E_{q} (r) + E_{Q} (r) = k [q \frac{( r - r _{q} )}{∣ r - r _{q} ∣ ^{3}} + Q \frac{( r - r _{Q} )}{∣ r - r _{Q} ∣ ^{3}}]$

以上のように、電荷が複数個あることで様々な形で電場が発生するわけだが、実際のところ電荷が物体内にあることを考えると物体表面から電場が生じるものと思われる。そこで、一つの電荷（電荷量 $Q$ ）の周囲に閉曲面があることを考えてみると、閉曲面が半径 $r$ の球面の場合は $E (r) 4 π r^{2} = 4 πk Q$ だけの電場が放出される。これが、任意の形状になると閉曲面 $S$ の微小面積が $d S^{'} = r^{2} sin ϕ d θ d ϕ$ となるのだが、実際は電場と平面に垂直な単位ベクトル $n$ のなす角度 $α$ により $d S = \frac{d S ^{'}}{cos α}$ だけの面積を貫く。

一方で、電場と

n

の内積が

E \cdot n = E cos α

であるが、面上の

E

の合計（積分）を求めるのに

cos α

が余計なので

\oint (E \cdot n) d S = \oint E d S^{'} = k Q \int_{0}^{2 π} d θ \int_{0}^{π} sin ϕ d ϕ = 4 πk Q

というようにすると、面上の電場の総和を求めることができる。ここで、左辺については一般的に

\oint (E \cdot n) d S = \oint E \cdot d S

というように面積分と呼ばれる形で書かれる。

　任意の閉曲面の結果から、閉曲面の形によらず電場の総量は電荷量に比例していることが予想される。そのため、２個の電荷がある場合にも $\oint E \cdot d S = \oint E_{q} \cdot d S + \oint E_{Q} \cdot d S = 4 πk (q + Q)$ というようになる。少しひねくれて、もし閉曲面の外部に電荷がある場合は電場が地点1から地点2に向けて出入りした場合は、 $\oint E \cdot d S = E_{1} cos (π - α_{1}) d S_{1} + E_{2} cos (α_{2}) d S_{2} = - E_{1} d S_{1}^{'} + E_{2} d S_{2}^{'} = 0$ である。したがって、閉曲面から出てくる電場というのは、その閉曲面内にある電荷量 $Q$ に比例することになる。これをGaussの法則と呼ぶ。 $\oint E \cdot d S = 4 πk Q$

光学

　昼は明るく夜は暗くなるという毎日の繰り返しや夏は暑く冬は寒くなるという四季の移り変わりも、太陽からくる日射量の変化によるものである。太陽の光により植物などが育ち、それを食料として数多くの生き物が生活していることは既にご存じのとおりである。その中でも、人類は長い歴史の中で自分の眼で周囲の状況を捉え身を守ったり食料を手に入れるということをしてきた。同時に、光は無意識のうちに人々にとって片時も欠かせないものであったと言える。実際、昔の人々がその存在を信じて疑わなかった神様や仏様の絵画や彫刻には必ずと言ってよいほど「後光」という光がその後ろに燦然と輝いている。ただ、その光景がそのまま眼に入っているわけでなく、外部からの光が眼の中にある水晶体という透明な物質でできたレンズを通り網膜の表面に像を作ることで物体が見えるようになっている。周囲の状況を正しくとらえるには、光とそれを集めて像を結ばせるレンズ、そしてその像を感知する網膜が欠かせないわけである。これと同じ原理で、Galileiは人工のレンズを巧みに組わせて望遠鏡を作り、肉眼でははっきり見えない遠くの小さな星を観察して、天文学に全く新しい発展をもたらした。一方でHookeもレンズの組み合わせを工夫して顕微鏡を作り、肉眼でははっきり見えない細かな微生物を観察して生物学を大きく進歩させた。一方で、晴れた夜空には星が輝いているが、特徴のある星が見える方向は時刻と方位を知ることのできる指標として、大海原を渡る船乗りや砂漠を横切る飛び人などに大昔から利用されてきた。このように自然の光の働きと輝きから、光は生きるのに必要な活力を与えてくれる生物に優しいエネルギーであると同時に、身の回りの状況を知らせてくれる実直な情報の運び手であることもよく分かる。

　以上のように、光に触れていくにつれて光を利用するための学問として光学が発展してきた。これを利用して、まず先ほど述べた望遠鏡や顕微鏡と同じ原理でカメラや夜を明るくしようと電気で灯すランプが登場した。19世紀の終わりにはカメラで撮影したフィルムで映像を映し出す映画が作られるようになり、それが普及するとともにランプが発端となって真空管を扱う電子工学が発展したことでテレビも生まれた。このあたりから、光学は「電気」と電気の運び手である「電子」の力を借りるようになった。電気工学、電子工学の分野で通信やコンピュータが発達してきたことを受けて、光学と電子工学の結びつきは決定的になってきました。このことは、光ファイバを使った光通信、CSなどの光ディスクメモリー、液晶などの薄型ディスプレイ、あるいはパソコンとも密接な関係にあるデジタルカメラや電子ゲームなどが急速に発展しつつあることから明らかである。このように、人々の日々の暮らしの中に徐々にではあるが光の能力を活かした製品は着実に入り込んできています。こうして、人工の光は社会の様々な所で使われるようになり、人々の生活環境が昔よりも格段に安全で快適なものになってきたと言えます。その光を取り扱うために、光学を学ぶ必要があるため見ていこう。

幾何光学

　光に関する研究は、既に遠くギリシャ時代においても行われ、例えば紀元前300年頃、Euclidは光学に関する本を著している。この本の中で、彼は光の直進性や反射の法則を彼の幾何学と結び付けて論じているが、このように光の本性は問題とせずに、光の進路を純粋に幾何学的に扱う分野を幾何光学と呼ぶ。幾何光学は中世の間にも大いに進歩し、球面鏡や放物面鏡、レンズなどを生み、老人用の眼鏡の発明なども行われた。さらにこれは、近代物理学の誕生期へと引き継がれ、Galileoによる屈折望遠鏡の製作（1609年）、Snellによる屈折の法則の発見（1621年）、Fermatによる光の進路を決める原理の確立（1657年）、Newtonによる反射望遠鏡の発明（1668年）などへとつながって行く。

　光の本性が本格的に問題とされるようになったのは、17世紀に入ってからであり、よく知られているように、光を微小な粒子の流れとみる粒子説と媒質中を伝搬する波動とみる波動説の間で、長い期間にわたって激しい論争が行われた。これは、例えばSnellの方に関しても、このどちらの考え方でも尤もらしい説明をすることが可能だからである。この法則は、二つの媒質の境界面に光が入射し、一部が屈折して別の媒質に侵入する場合、以下の関係式が $θ_{1}$ によらず媒質が固定されていれば決まった値をとるというものである。

\frac{sin θ _{1}}{sin θ _{2}} = n_{12}

これに対する粒子説の説明は次のように行われる。すなわち、二つの媒質ではポテンシャルエネルギーの差があり、このために粒子は境界面の所で面に垂直な方向に力を受ける。しかしこの場合、面に平行な方向の速度成分は変わらないから、それぞれの媒質中での光の速さを $v_{1}, v_{2}$ として、 $v_{1} sin θ_{1} = v_{2} sin θ_{2}$ が成立する。したがって、 $sin θ_{1}, sin θ_{2}$ の比は $θ_{1}$ に依存せず $\frac{sin θ _{1}}{sin θ _{2}} = \frac{v _{2}}{v _{1}}$ というようにSnellの法則での $n_{12}$ （相対屈折率）は速度の比 $v_{2} / v_{1}$ に等しいということになる。

　一方で、Huygensの原理を用いた波動説による説明は以下のとおりである。すなわち、図のように平面波が上方からやってきたとして、平面波が点Aから点Bに届くまでの時間 $Δ t$ を考えてみると、その間にできる波面というのは $Σ$ のようになり、関係としては $\overline{A B} sin θ_{1} = v_{1} Δ t 、 \overline{A B} sin θ_{2} = v_{2} Δ t$ が成り立つ。従って、それぞれの比をとることで $\frac{sin θ _{1}}{sin θ _{2}} = \frac{v _{1}}{v _{2}}$ となり、先ほどとは逆の関係（ $n_{12} = v_{1} / v_{2}$ ）になることが分かる。

どちらが正しいかどうかは空気中と水中での光の速さを比較すれば、粒子説と波動説のどちらが正しいか決めることができる筈である。実際、光が空気から水の中へと屈折する場合 $n_{12} > 1$ であるが、1850年、Foucaultは水中での光速が空気中でのそれよりも遅いこと（ $v_{2} < v_{1}$ ）を確かめ、これにより波動説が勝利をおさめる結果となった。

　上記のことをふまえて反射、屈折の法則をまとめると次のようになる。すなわち、入射光線と境界面との交点Oに立てた境界面に対する法線（接平面に垂直な直線）OQと、入射光線AOとで作られる平面を光の入射面と呼ぶが、反射光線も屈折光線もともにこの平面内にある。

そこで、この面内で図のように角度を決めると、入射角 $θ_{1}$ と反射角 $θ_{r}$ とは等しく、入射角 $θ_{1}$ と屈折角 $θ_{2}$ の間には入射側と屈折側それぞれの屈折率を $n_{1} 、 n_{2}$ として、次の関係（Snellの法則）が成り立つ。

$n_{1} sin θ_{1} = n_{2} sin θ_{2}$

ここで屈折率 $n$ は媒質中の光速 $v$ 、真空中の光速 $c$ と $n = c / v$ という関係があるため、光速に置き換えてみると先ほどの波動説の説明のところに出てきた式と同じ形になっていることが分かる。

$\frac{sin θ _{1}}{sin θ _{2}} = \frac{v _{1}}{v _{2}}$

また、屈折率は様々な色の光に対して異なった値を持つが、その違いというのは僅かであるから通常は光に対する媒質の屈折率は一定としても差支えがない。また、真空の屈折率が1であるが、空気の屈折率が1.00029であるから1として見なしてもよい。

　ここまでは $n_{2} > n_{1}$ に限った話であったが、一方で $n_{2} < n_{1}$ の場合だと入射角よりも屈折角の方が大きくなるから、入射角が大きくなると屈折角が90°をこえる場合も考えられる。しかし、そのように屈折することはなく、実際は全て反射（全反射）されることになる。また、全反射が始める入射角のことを臨界角（critical angle）と呼ばれる。この角度は屈折角が90°となったときの以下の式から求められる。 $θ_{c} = arcsin (\frac{n _{2}}{n _{1}}) 、 \frac{n _{2}}{n _{1}} < 1$ 例えば、媒質1を水、媒質2を空気とすると $n_{1} \approx 4/3, n_{2} \approx 1$ であるから、臨界角は約50°になる。 $θ_{c} = arcsin (\frac{3}{4}) = 48.59 \dots \approx 50°$ したがって、例えば以下の図のように水中に沈めた10円玉Pに光が反射し、その光が水中から空気中に出る際に屈折して目に届いたとする（このとき10円玉Pは点線で描いたように浮いて見えた状態になっている）。仮に $- 50° < r < 50°$ の範囲であれば屈折して10円玉が見えるが、それを超えると全反射をしてしまうため見えなくなることが分かる。

Fermatの原理

　屈折率が場所によって僅かずつ異なるような物質中を光がとおるときには、光線は曲がるがその場合の光の進路を決める法則はFermatの原理として知られている。これは、光の進む道はその近傍の経路と比べて光学距離が極値（正しくは停留値）になるようなものであるという表現される。光学距離とは、屈折率を経路に沿って積分したもので、 $r$ を空間座標、 $d s$ を経路に沿った不精な線要素とすると、A点からB点までの経路では、 $L (A, B) = \int_{A}^{B} n (r) d s$ と書くことができる、これは、AB間の経路の距離が光学的には真空中を伝わる平面波に換算して $L (A, B)$ であること、つまりB点の位相のA点のそれに比べた遅れが、真空中を平面波が距離 $L (A, B)$ だけ行ったときの億例に等しいことを意味している。この原理は、実際の経路についてはそれからの任意の無限小のずれに対して光学距離の変化分が零、つまり $δ L (A, B) = 0$ が成り立つと表現することもできる。さらに、光の速さは媒質の屈折率に反比例するから、もっと簡単には実際に光がとおる道はそれに少し変えた経路と比較して、同じ2点間を行くのに要する時間が極値（ないしは停留値）になるようなものである。といってもよい。

　Fermatの原理は、幾何光学の公理ともいうべきものであり、前節で述べた光線の可逆性や均一な媒質中での光の直進性はもちろん、光の反射や屈折の法則もこの原理の中に含まれる。すなわち、反射光線と屈折光線が入射面内にあることは、この原理から自明であるから、二つの媒質（屈折率 $n_{1}, n_{2}$ ）の境界面が $x z$ 面（ $x$ 軸を含む紙面に垂直な面）であり、光の入射面が $x y$ 面であるとして、以下の図のように座標を決める。

このとき反射、屈折のそれぞれの場合の経路AOA'およびAOBの光学距離は

L (A, A^{'}) = n_{1} (x^{2} + a^{2} + (b - x)^{2} + c^{2})

L (A, B) = n_{1} (x^{2} + a^{2} + (d - x)^{2} + e^{2})

となる。そこで、

x

を変えた時にこれらが極値となることから、それぞれ

δ L (A, A^{'}) = n_{1} (\frac{x}{x ^{2} + a ^{2}} - \frac{b - x}{( b - x ) ^{2} + c ^{2}}) δ x = 0

δ L (A, B) = (\frac{n _{1} x}{x ^{2} + a ^{2}} - \frac{n _{2} ( d - x )}{( d - x ) ^{2} + e ^{2}}) δ x = 0

となり、これらについて以下の関係が成り立つ。

sin θ_{1} = \frac{x}{x ^{2} + a ^{2}} 、 sin θ_{r} = \frac{b - x}{( b - x ) ^{2} + c ^{2}} 、 sin θ_{2} = \frac{( d - x )}{( d - x ) ^{2} + e ^{2}}

そのため、以下のように反射・屈折の法則が導かれる。

sin θ_{1} = sin θ_{r} （ θ_{1} = θ_{r} ）、 n_{1} sin θ_{1} = n_{2} sin θ_{2}

レンズの仕組み

　ところで、そもそも人の目というのは物体から反射された光をどのようにして認識しているであろうか。構造としては以下の左図のようになっており、角膜・水晶体のところで屈折をして、網膜上に光が集まることで、網膜上にはっきり画像が映し出される。その映し出された画像情報が脳に伝達されることで、物体がはっきり見えるという流れになる。この構造というのはレンズで再現することが可能であり、実際にカメラなどでは以下の右図のようになっている。

このことからも分かるように、レンズの仕組みを理解することが重要なため、順を追ってそのことについて述べていくことにする。まず曲面のガラスに光が当たるとどうなるかを考察してみる。曲面に光線があたると、光は以下の赤線のような軌道で反射をすると考えられる。

今の場合、点Aから点Mを介して点Bへ通る状況を考えているが、同じように点Pを介しても点Bに到達することができる（どちらも、Fermatの原理で言う $δ L (A, B) = 0$ を満たす）。すると、これら2つの光学距離は同じものとなるため、以下の式が成り立つ。 $(a - y)^{2} + x^{2} + (b - y)^{2} + x^{2} = a + b$ しかし、この式は大変複雑であるため、簡略化するため $a$ や $b$ に比べて $x$ や $y$ の絶対値が十分に小さい範囲の光線のみについて考えてみる。すると、左辺の各項目が $(a - y)^{2} + x^{2} = a 1 + \frac{x ^{2} + y ^{2} - 2 a y}{a ^{2}} \approx a + \frac{x ^{2} + y ^{2} - 2 a y}{2 a}$ $(b - y)^{2} + x^{2} = b 1 + \frac{x ^{2} + y ^{2} - 2 b y}{b ^{2}} \approx b + \frac{x ^{2} + y ^{2} - 2 a y}{2 b}$ となることから、以下の式に置き換えることができる。 $x^{2} + y^{2} - \frac{4 ab}{a + b} y = 0$ そして、さらにこの式を変形して $x^{2} + (y - R)^{2} = R^{2} 、 (R = \frac{2 ab}{a + b})$ とすることで円の方程式になっていることがわかるため、曲面鏡が半径 $R$ の球の一辺とみなせるものと考えられる。このときの $R$ というのは曲率半径と呼ばれ、この式をさらに変形することで $\frac{1}{a} + \frac{1}{b} = \frac{1}{f} 、 (f = \frac{R}{2})$ と書くことができる。これは、鏡の中心（極）Pと曲率半径Oを結ぶ線（主軸）に近い光線（近軸光線）に限った話であることに注意されたし。すなわち、主軸上にある点Aから出た光は球面により上記の式を満たす地点Bに反射して届くと結論づけられる。

解析力学

Lagrange形式

　以前まで、物体の運動は運動方程式を用いて見てきたわけだが、拘束力が働く場合は、その方向とは垂直に運動を行うため結果的に方程式から姿を消すことが多い。そのため、予め拘束力を考える必要がない方程式を考えると便利なように思う。そこで、拘束力 $R$ とは垂直の方向に $δ r$ だけの変位（同じ時間で二つの位置の差を測る）を考えると、以下のように仕事が0になることが分かる。 $R \cdot δ r = 0$ これは仮想仕事の原理というように呼ぶ。次に、運動方程式では $m \frac{d v}{d t} = R + F$ というように拘束力と外力とで力を分けて記載できるため、 $(m \frac{d v}{d t} - F) \cdot δ r = 0$ となる。次に、このベクトルが $f$ 個の変数で書けるものとして $r (q_{1}, q_{2}, \dots, q_{f})$ とすると

$δ r = i = 1 \sum f (\frac{\partial r}{\partial q _{i}}) δ q_{i}$

となるので、以下のようになる。 $i = 1 \sum f [m \frac{d v}{d t} \cdot (\frac{\partial r}{\partial q _{i}}) - F \cdot (\frac{\partial r}{\partial q _{i}})] δ q_{i} = 0$

ここで、第1辺を整理して

$m \frac{d v}{d t} \cdot (\frac{\partial r}{\partial q _{i}}) = \frac{d}{d t} [m v \cdot (\frac{\partial r}{\partial q _{i}})] - m v \cdot \frac{d}{d t} (\frac{\partial r}{\partial q _{i}})$

とすると、各項目について

$\frac{\partial v}{\partial q ˙ _{i}} = \frac{\partial}{\partial q ˙ _{i}} [j = 1 \sum f (\frac{\partial r}{\partial q _{i}}) \overset{q}{˙}_{i} + \frac{\partial r}{\partial t}] = \frac{\partial r}{\partial q _{i}}$ $\frac{d}{d t} (\frac{\partial r}{\partial q _{i}}) = j = 1 \sum f \frac{\partial}{\partial q _{j}} (\frac{\partial r}{\partial q _{i}}) \overset{q}{˙}_{j} + \frac{\partial}{\partial t} (\frac{\partial r}{\partial q _{i}}) = \frac{\partial}{\partial q _{i}} [j = 1 \sum f (\frac{\partial r}{\partial q _{j}}) \overset{q}{˙}_{j} + \frac{\partial r}{\partial t}] = \frac{\partial v}{\partial q _{i}}$ となるため、これらを代入すると以下の通りになる。 $i = 1 \sum f [\frac{d}{d t} \frac{\partial}{\partial q ˙ _{i}} (\frac{1}{2} m v^{2}) - \frac{\partial}{\partial q _{i}} (\frac{1}{2} m v^{2}) - F \cdot (\frac{\partial r}{\partial q _{i}})] δ q_{i} = 0$ 任意の仮想変位に対して成立するためには、中の項目が0になる必要があるため

$\frac{d}{d t} (\frac{\partial K}{\partial q ˙ _{i}}) - \frac{\partial K}{\partial q _{i}} = Q_{i} (i = 1, \dots, f)$ $K = \frac{1}{2} m v^{2} 、 Q_{i} = F \cdot (\frac{\partial r}{\partial q _{i}})$ が成り立つ。この運動方程式は一般的な座標に関して成立するものであり、Euler=lagrange方程式と呼ぶ。ただ、通常は力を保存力として $F = - \frac{\partial U}{\partial r}$ として、一般化力を以下のようにする。 $Q_{i} = - \frac{\partial U}{\partial r} \cdot (\frac{\partial r}{\partial q _{i}}) = - \frac{\partial U}{\partial q _{i}}$ すると、以下のように書くこともできる。 $\frac{d}{d t} (\frac{\partial L}{\partial q ˙ _{i}}) - \frac{\partial L}{\partial q _{i}} = 0 、 \frac{d}{d t} (\frac{\partial U}{\partial q ˙ _{i}}) = 0$ $L = K - U （ラグランジアン）$

Hamilton形式

Lagrange形式より対称性を高めたもの

一般化座標への変換を含む、ある種の $q_{i}, Q_{j}$ についての変換に対して運動方程式を不変とする正準方程式がある。 $q_{i} = q_{i} (Q_{j}, \dot{Q}_{j}) 、 \overset{q}{˙}_{i} = \overset{q}{˙}_{i} (Q_{j}, \dot{Q}_{j})$ ただし、 $(q_{i}, \overset{q}{˙}_{i})$ に対しての変換に関する理論はない。しかし、Lagrange形式でも出た $p_{i} = \frac{\partial L}{\partial q ˙ _{i}} \to \overset{q}{˙}_{i} = \frac{\partial G}{\partial p ˙ _{i}}$ ということを行うと、Lagrange形式のものより対称性が高い正準方程式が導かれる（Legendre変換）。 $G (q_{i}, p_{i}) = j \sum \overset{q}{˙}_{j} p_{j} - L$ この $G$ は一般的に $H$ と表わしHamiltonianと呼ばれている。

方程式の導出方法

ハミルトニアンの定義とオイラーラグランジュ方程式から

統計力学

　これから学ぼうとする統計力学では、物質が非常に多くの小さい粒子（原子・分子）から成り立つものと仮定して、これらの分子が古典力学に従い運動することでどのような性質（運動ではない）を示すかを研究するものとなっている。ここで言う原子（分子）は、実際のところ複雑な内部構造や相互作用を有しているが、そこまで考慮すると理論の構築が困難となるため次のような理想化を行う。

　分子は半径 $r_{1}$ の滑らかな表面を持った剛体弾性球であると考え、二つの分子は距離が半径の和 $r_{12} = r_{1} + r_{2}$ に接近したときにのみ斥力が働くものとする。二分子間のポテンシャルを相互距離 $r$ の関数とすれば、以下のグラフのようになる。しかしこの剛体平滑分子の模型はあまりにも素朴すぎるかもしれないので、次のような模型を考えることもある。分子を質点とみなし、中心力を及ぼしあうことを考える。二つの分子間の距離 $r$ のみの関数であるポテンシャルエネルギーを $U (r)$ とすれば以下のグラフのようになる。

　次に、これらの分子は温度の高低によって、運動が速やかになったり緩やかになったりする。そして互いに衝突しあい作用を及ぼしあう（運動量およびエネルギーを交換し合っている）。このとき、これらの気体分子の全体としての平均の並進エネルギーを除いた分子の運動エネルギーが気体の熱エネルギーであるものとする。なお、気体分子が回転、振動などの内部自由度をもつならば、その運動エネルギー、ポテンシャルエネルギーもまた熱エネルギーに含まれるものとする。以上の二つの仮定をもとに、通常の力学のように与えられた初期条件の下で運動方程式を積分すればよいのであるが、今考えている体系では非常に多くの分子から成り立っており、個々の分子に対して一々初期条件を知ることも、運動方程式を解くことも不可能である。それゆえに我々の問題は、適当な微小時間において平均化された、ある容積の中に存在する分子の平均数、平均エネルギー、あるいはこれらの量の時間的変化というような平均的な量のみ考えることにする。これは単に数学的な面ばかりでなく、物理的にも我々の測定しうる量は、このような平均的な量のみであることからもうなずける。このような問題を解くためには、後で分かるように統計的な考察を必要とする。

　おおよそ、上のような考え方で物質、気体の物理的性質を説明しようとする原子論の立場は、古くはギリシャ時代にまでさかのぼるのであるが、定量的な理論の形を取り始めたのは、もちろんNewton力学のできた後である。特に19世紀半ば頃MaxwellやBoltzmannによって一応完結された。その基礎に立って詳しい計算を行ったのは、EnskogおよびChapmanである。これらの理論の基底となっているのは、もちろん古典力学であるが、原子級の問題では量子力学を使わなければいけない。それで以下に述べる理論もその意味で改良しなければならいないのであるが、改良の結果は大したことでないようであるから説明を省く。この量子力学を取り入れ方は統計力学の場合と同じように二重である。一つは統計法の変化にあり、もう一つは衝突の確率に関することである。後者に関する限りは、古典理論に取り入れることはやさしい。理論の大部分が衝突の確率の詳しい知識には関係なく組み立てられるからである。気体運動論の方法を金属の中の自由電子の気体に応用しようとする場合には、上の二つの量子論的効果を十分に考慮しなければならない。

気体分子運動論

　気体が分子で構成されているものとみなした場合、容器内ではまず壁への衝突による力（圧力）が働くことが考えられる。そのため容器の壁に及ぼす圧力を考えると、圧力 $P$ は壁の単位面積が受ける力（撃力）の垂直成分の時間平均として定義される。そのため力の時間平均は分子の受ける力積の単位時間内総量に等しく、更にこれは運動量の変化の単位時間内総量に等しい。速さ $v$ の分子が壁に衝突して、弾性的に反発されるとすれば分子の運動量変化は壁に垂直に $2 m v_{⊥}$ となる。そのため、単位時間当たりの力は以下のようになる。 $\overset{ˉ}{F} = \frac{2 m v _{⊥}}{Δ t} 、 (v = v_{⊥} + v_{//})$ ただ、実際は一つの分子による力のみがかかるわけではなく、同じ速度の分子の数（ $N (v_{⊥})$ ）だけ力が働く。この個数というのは時間 $Δ t$ の間に範囲 $Δ S$ に力が生じるものとして、数密度 $n (v_{⊥})$ により $N (v_{⊥}) = n (v_{⊥}) (\frac{1}{2} v_{⊥} Δ t Δ S)$ となるため、全体の力は以下の通りとなる。 $F (v_{⊥}) = \overline{F} N (v_{⊥}) = mn (v_{⊥}) v_{⊥}^{2} Δ S$ 最終的な圧力は、全ての速度に対する力の合計を面積で割ったものであるため次のように書ける。 $P = \frac{1}{Δ S} v_{⊥} > 0 \sum F (v_{⊥}) = m v_{⊥} > 0 \sum n (v_{⊥}) v_{⊥}^{2}$ ここで、単位体積内の分子数を $n$ としたとき $P = mn ⟨ v_{⊥}^{2} ⟩ 、 ⟨ v_{⊥}^{2} ⟩ = v_{⊥} > 0 \sum v_{⊥}^{2} f (v_{⊥}) 、 (f (v_{⊥}) = \frac{n ( v _{⊥} )}{n})$ となる。ここで $⟨ v_{⊥}^{2} ⟩$ がどのくらいの値になるかは式を変形することで求めることができる（数密度： $n (= N / V)$ 、全質量： $M$ 、体積： $V$ ）。 $P = \frac{m N ⟨ v _{⊥}^{2} ⟩}{V} = \frac{M ⟨ v _{⊥}^{2} ⟩}{V} \to ⟨ v_{⊥}^{2} ⟩ = \frac{P V}{M}$ 例として、標準状態（0℃、1気圧）における1モルの気体を考えると、 $P = 1.013 \times 1 0^{5} N/ m^{2} 、 V = 0.0224 m^{3}$ であるため、各分子ごとに以下のようになる。

分子名	$M$ （kg）	$⟨ v_{⊥}^{2} ⟩$
H $_{2}$	0.0020	1134560
He	0.0040	567280
N $_{2}$	0.0280	81040
O $_{2}$	0.0320	70910

そのため、分子の平均速度はおおよそ $100 ～ 300 m/s$ となるわけだが、実際の気体の拡散速度（1 m/s程度以下）と比べるとずっと遅いため矛盾があるように感じる。これは衝突を繰り返すことで平均自由行路が短くなり、全体的な速度が遅くなるということによるものと考えられる。そのため、衝突に関して考察していくことにする。

　気体の分子には大きさがあるために必ず衝突する。分子の直径を $σ$ とし、一つの分子 M が M $_{1}$ 、M $_{2}$ 、M $_{3}$ ,・・・と次々に衝突して運動を変えていくありさまを想像したものが図である。

　さて、この場合に衝突の数を知るためにはMが半径 $σ$ の球、その他が半径0になったとしても、全く同じである。そこで、Mが全行程 $L$ を通過し終える間に何回他の分子と衝突するかということは、切り口の半径が $σ$ で中心軸の全長 $L$ の円柱の中に何個の分子中心がああるかということに一致する。そこで、いま数密度が $n$ であったとすると、円柱の体積が $π σ^{2} L$ であるから、 $衝突回数 = nπ σ^{2} L$ となる。ここで、衝突回数が1に相当する距離 $L$ は一つの衝突から次の衝突までの平均距離と考えられる。この距離を平均自由工程という。 $l = \frac{1}{nπ σ ^{2}}$ Maxwellは自身で導いた粘性係数 $η$ と平均自由工程の関係式 $η = \frac{1}{3} ρl \overset{v}{ˉ} 、 (\overset{v}{ˉ} = ⟨ v_{⊥}^{2} ⟩)$ を用いて、常温での平均自由工程の値を以下の値で推定している。 $l ≃ 5.7 \times 1 0^{- 8} m$ そのため、1sの間に分子が100m進んだ場合には $衝突回数 = \frac{L}{l} = \frac{100}{1 0 ^{- 8}} = 1 0^{10} 回$ というように多くの衝突が起きていることがうかがえる。さらに、Loschmidtは1m $^{3}$ の気体を圧縮して液体にしたときの体積 $V_{0}$ は、気体の分子の数密度（1m $^{3}$ に含まれる分子の個数）を用いて $V_{0} ≃ n \cdot \frac{4}{3} π (\frac{σ}{2})^{2}$ となると考え、空気分子1個の直径は $σ ≃ 1 0^{- 9} m （ 10 \overset{˚}{A} ）$ となることを推定している。これを踏まえると、数密度も $n ≃ 2 \times 1 0^{24} m^{- 3}$ となるわけだが、0℃、1気圧での1m $^{3}$ 中での分子数はロシュミット数と呼ばれている。この値を見ても分かるように、かなり多くの分子が存在していることが伺える（この値だと不正確であるため、現在ではより正確な値であるアボガドロ数が用いられる）。

では、これを進めて議論していきたいところであるが、まともに議論すると複数個の質点全ての運動を見ることに相当するので、現実的に解けるような問題でないことは予想できるであろう。そのため、これからは統計集団という考えを用いて複数個の分子を取り扱っていくことにする。

統計集団（小正準集団）

　気体の性質を調べる際には、よくピストンなどの容器に気体を閉じ込めておくことが多いが、このときに閉じ込められた気体（分子の集まり）は一種の集団を形成している。この集団を統計集団と呼ぶ。例えば、この気体が外部との熱のやり取り、漏れ出しなどもなく完全に孤立している集団（小正準集団）だとすると、熱力学第1法則 $d U = d^{'} Q + d W$ は右辺が0となるため、内部エネルギーが保存されることが分かる。さらに内部で熱平衡になっている場合にはエントロピーの増大もしなくなるため $d^{'} Q = T d S = 0$ であり、なおかつ容器に閉じ込められている関係で仕事も $d^{'} W = - P d V = 0$ と体積が変化しない関係で一定になる。体積については気体の性質とは直接関係がない（人的に変更されるもの）が、エントロピーについてはBoltzmannによると気体のとりうる状態数 $W$ により $S = k_{B} ln W$ というような関係があるとされている（ボルツマンの原理）。この状態数というのは全体の状態がエネルギーを用いて書けることから計算することができる。具体的には、解析力学においてエネルギー $U$ はハミルトニアン $H$ で記述される。例として、 $N$ 個の粒子がある場合は $i$ 番目の粒子の一成分の位置と運動量を $q_{i}, p_{i} (i = 1, \dots, 3 N)$ とすると $U = H (q_{1}, \dots, q_{3 N}, p_{1}, \dots, p_{3 N})$ というように $6 N$ の変数を持つ関数になる。これら $q . p$ の組（共役変数）が存在する空間のことを位相空間と呼ぶ。例えば、自由粒子の場合だとハミルトニアンが $H = i = 1 \sum 3 N \frac{p _{i}^{2}}{2 m}$ と $p_{1}, \dots, p_{3 N}$ のみ依存した形になる。今回はエネルギーが一定値 $U$ ということなので、この式を変形すると $i = 1 \sum 3 N p_{i}^{2} = 2 m U$ であることから、 $3 N$ 次元の球の方程式になっていることが分かる。簡単のため、粒子が一つで1次元の場合の状態数を考えると $q_{1}, p_{1}$ しかないため $p_{1}^{2} = 2 m U$ として、これを満たす $q_{1}, p_{1}$ の状態数（組み合わせの数）を数えていく。しかし、小数点以下も含めると無限個存在することになってしまうため、 $q_{1}, p_{1}$ がそれぞれ $L 、 P$ の幅しかとらず十分小さい間隔 $δ q 、 δ p$ ごとに一つの状態があるものとすると $W_{1} = \frac{L}{δ q} \times \frac{P}{δ p} = \frac{L P}{δ q δ p}$ として状態数を記述することができる。同様に2次元になった場合でも $W_{2} = \frac{L ^{2}}{δ q ^{2}} \times \frac{P ^{2}}{δ p ^{2}} = \frac{L ^{2} P ^{2}}{( δ q δ p ) ^{2}}$ となることから、 $3 N$ 次元になった場合は以下の形になる。 $W_{3 N} = \frac{V ^{N} P ^{3 N}}{h ^{3 N}} 、 V = L^{3} 、 h = δ q δ p$ しかし、自由粒子の場合は運動量が $3 N$ 次元の球の表面上に制限されていることから $W_{1} = \frac{L}{δ q} \times 2 、 W_{2} = \frac{L ^{2}}{δ q ^{2}} \times 2 π \frac{2 m U}{δ p}$

$3 N$ 次元になった場合は半径 $r$ の $f$ 次元の球の表面積が $S_{f} (r) = \frac{f π ^{f /2}}{Γ ( f /2 + 1 )} r^{f - 1}$ であることから、状態数は以下の通りとなる。 $W = \frac{L ^{3 N}}{δ q ^{3 N}} S_{3 N} (2 m U / δ p) = \frac{V ^{N}}{h ^{3 N}} \frac{2 π ^{3 N /2} δ p}{Γ ( 3 N /2 )} (2 m U)^{3 N /2 - 1}$

　一方で、熱力学第1法則をエントロピーとエネルギーを含む式に展開すると $d S = \frac{1}{T} d U + \frac{P}{T} d V$ であるため、以下の関係式が成り立つことが分かる。 $(\frac{\partial S}{\partial U})_{V} = \frac{1}{T} 、 (\frac{\partial S}{\partial V})_{U} = \frac{P}{T}$ そのため、まずエントロピーを求めてみると $S = N k_{B} ln V + \frac{1}{2} (3 N - 1) k_{B} ln (2 m U) + k_{B} ln \frac{2 π ^{3 N /2} δ p}{h ^{3 N} Γ ( 3 N /2 )}$ であることから、それぞれの式に代入すると $(\frac{\partial S}{\partial U})_{V} = \frac{1}{T} \to U = \frac{1}{2} (3 N - 1) k_{B} T ≃ \frac{3}{2} N k_{B} T$ $(\frac{\partial S}{\partial V})_{U} = \frac{P}{T} \to P V = N k_{B} T$ という二式が求められる。このうち、エネルギーをハミルトニアンと見比べると $i = 1 \sum 3 N \frac{p _{i}^{2}}{2 m} = i = 1 \sum 3 N \frac{1}{2} m v_{i}^{2} = \frac{3}{2} N k_{B} T \to \frac{1}{3 N} (i = 1 \sum 3 N \frac{1}{2} m v_{i}^{2}) = \frac{1}{2} k_{B} T$ というように、一成分当たりの平均エネルギーが $1/2 k_{B} T$ になっていることが分かる（エネルギー等分配則）。そこで、これを期待値として $⟨ v^{2} ⟩ = \frac{1}{3 N} i = 1 \sum 3 N v_{i}^{2} = \frac{k _{B} T}{m}$ というようにすると、圧力に関する式を得ることができる。 $P V = M ⟨ v_{⊥}^{2} ⟩ \to P V = N k_{B} T$ さらに、 $n$ モルでの理想気体の状態方程式は $P V = n RT$ であることを踏まえると、次のような関係があることが導ける。 $R = N_{A} k_{B} 、 N_{A} = \frac{N}{n}$ ここでの $N_{A}$ は1molでの分子数を表しておりアボガドロ数と呼ばれている。そして、この関係式から定積モル比熱（1molあたりの比熱）を求めてみると次のようになる。 $C_{V} = (\frac{\partial U}{\partial T})_{V} = \frac{3}{2} R$

この値と実験値を比較してみると以下の表の通りになる。

分子名	$C_{V}$ （実験値）	$3/2 R$
He（気体）	$12.60 J/K \cdot mol (291℃) 12.33 J/K \cdot mol (93℃)$	$12.47 J/K \cdot mol$
A（気体）	$12.86 J/K \cdot mol (288℃) 12.00 J/K \cdot mol (93℃)$	$12.47 J/K \cdot mol$

単原子分子の場合だと近い値をとっていることが見て取れるが、二原子分子などになると当てはまらなくなる。これは、単原子分子のときと自由度が異なるためで、例えば二原子分子の場合だと $5 N$ だけある（チアバトンのように棒の両端に原子があるイメージ）。

分子名	$C_{V}$ （実験値）	$5/2 R$
N $_{2}$ （気体）	$20.83 J/K \cdot mol (293℃) 19.78 J/K \cdot mol (92℃)$	$20.78 J/K \cdot mol$
O $_{2}$ （気体）	$20.84 J/K \cdot mol (293℃) 18.43 J/K \cdot mol (92℃)$	$20.78 J/K \cdot mol$

同様に、固体についても調和振動子と考えると自由度が $6 N$ だけあり、なおかつDulong-Putiの法則においても $3 R$ になることは確かめられている。

分子名	$C_{V}$ （実験値）	$3 R$
Al（固体）	$24.37 J/K \cdot mol (15 ～ 100℃ 間の平均)$	$24.94 J/K \cdot mol$
Pb（固体）	$22.70 J/K \cdot mol (15 ～ 100℃ 間の平均)$	$24.94 J/K \cdot mol$

正準集団（カノニカルアンサンブル）

　前回までは分子の集団が孤立している状態を考えていたが、実際のところ孤立系というのは現実にはなく最低でも外部と熱のやり取りをしていることが多い。そのため、以下の図ような系を考えてみる。

ここでの容器系の状態を見てみると、熱力学第一法則と第二法則により $d U = d^{'} Q + d^{'} W = d^{'} Q 、 \frac{d ^{'} Q}{T} \leq d S$ となるため、先ほどと異なりエネルギーが一定とならないため状態数を求めることができない。一方で、熱浴系においても $d U_{b} = d^{'} Q_{b} + d^{'} W_{b} = d^{'} Q_{b} 、 \frac{d ^{'} Q _{b}}{T _{b}} \leq d S_{b}$ となるわけだが、全体としては孤立系であることを踏まえると、全体の内部エネルギーが保存しており尚且つ熱浴系と容器系の間でしか熱のやりとりしないことになる。 $d U_{t} = d U_{b} + d U = 0 、 d^{'} Q_{b} = - d^{'} Q$ ここで熱力学第二法則の二式の和をとるとエントロピーに関して $\frac{d ^{'} Q}{T} + \frac{d ^{'} Q _{b}}{T _{b}} = (\frac{1}{T} - \frac{1}{T _{b}}) d^{'} Q \leq d S + d S_{b}$ となるわけだが、ここで熱と温度の関係を見てみると ${・ T < T_{b} (1/ T - 1/ T_{b} > 0) の場合・ T > T_{b} (1/ T - 1/ T_{b} < 0) の場合 \to d^{'} Q > 0 \to d^{'} Q < 0$ というように互いに同符号になっていることから、以下の関係が成り立つことが分かる。 $0 \leq (\frac{1}{T} - \frac{1}{T _{b}}) d^{'} Q \leq d S + d S_{b}$ ここでBoltzmannの原理よりエントロピーを状態数で置き換えると $d S + d S_{b} = d (k_{B} ln W) + d (k_{B} ln W_{b}) = d (k_{B} ln W W_{b})$ となり、一つのエントロピー $S_{x}$ としてこれが増大し続けることが分かる。 $S_{x} = k_{B} ln W W_{b} \to 0 \leq d S_{x}$ そして、そのまま放置するとやがて熱平衡になり熱の移動が緩やかになることから、準静的変化をするようになり近似的に以下の式が成り立つ。 $d S ≃ \frac{d ^{'} Q}{T} 、 d S_{b} ≃ \frac{d ^{'} Q _{b}}{T _{b}}$ そのため、内部エネルギーに関する式もエントロピーを含めた形で表現することができる。 $d S = \frac{1}{T} d U + \frac{P}{T} d V \to (\frac{\partial S}{\partial U})_{V} = \frac{1}{T} 、 (\frac{\partial S}{\partial V})_{U} = \frac{P}{T}$ $d S_{b} = \frac{1}{T _{b}} d U_{b} + \frac{P _{b}}{T _{b}} d V_{b} \to (\frac{\partial S _{b}}{\partial U _{b}})_{V_{b}} = \frac{1}{T _{b}} 、 (\frac{\partial S _{b}}{\partial V _{b}})_{U_{b}} = \frac{P _{b}}{T _{b}}$ また、二つのエントロピーの和（ $S_{x}$ ）も熱平衡（ $T ≃ T_{b}$ ）の状態になると $d S + d S_{b} = (\frac{1}{T} - \frac{1}{T _{b}}) d^{'} Q ≃ 0 \to d S_{x} = 0$ と増大しなくなり、極大値をとると考えられる。ここで $W_{x}$ は熱浴系も容器系も孤立系の場合だと、互いに独立していることから全体の状態数になるものと思われる。 $W_{t} = W_{x} = W W_{b}$ しかし、今の場合だと各々でエネルギーが様々な値をとることから、全体のエネルギーが一定の条件の下で以下のような総和をとった形になると考えられる。 $W_{t} (U_{t}) = U t = U + U_{b} \sum W (U) W_{b} (U_{b}) = U \sum W (U) W_{b} (U_{t} - U)$ ここで、左辺で右辺を割るとあたかも確率のようになることが分かる。 $1 = U \sum f (U) 、 f (U) = \frac{W ( U ) W _{b} ( U t - U )}{W _{t} ( U _{t} )}$ そして、右辺の $W_{b}$ のエントロピーに関して $S_{b} (U_{t} - U) ≃ S_{b} (U_{t}) - (\frac{\partial S _{b}}{\partial U _{b}})_{V} U = S_{b} (U_{t}) - \frac{U}{T _{b}}$ であることから、状態数に関しては以下のように書ける。 $W_{b} (U_{t} - U) = W_{b} (U_{t}) e^{- U / k_{B} T_{b}}$ そして、熱平衡（ $T_{b} ≃ T$ ）であることを踏まえると確率は次の通りになる。 $f (U) = \frac{W _{b} ( U _{t} )}{W _{t} ( U _{t} )} W (U) e^{- β U} 、 β = \frac{1}{k _{B} T}$ ここで、この総和をとったときに $\frac{W _{b} ( U _{t} )}{W _{t} ( U _{t} )} = \frac{1}{Z} 、 Z = U \sum W (U) e^{- β U}$ とすることで、以下のように置くこともできる。 $f (U) = \frac{1}{Z} W (U) e^{- β U}$ あるいは、状態数をエントロピーの形に置き換えることで $f (U) = \frac{1}{Z} e^{- βF} 、 F = U - TS$ とHelmholtzの自由エネルギー $F$ で表現することもできる。この関数は微小変化をとると $d F = d U - S d T - T d S = - S d T - P d V$ であるから、以下の関係式が成り立つ。 $S = - (\frac{\partial F}{\partial T})_{V} 、 P = - (\frac{\partial F}{\partial V})_{T} 、 U = F + TS = - T^{2} [\frac{\partial}{\partial T} (\frac{F}{T})]_{V}$ このことから、自由エネルギーの形で表示できることが分かる。そのため、期待値の方も $⟨ S ⟩ = - \frac{1}{T} (\frac{\partial}{\partial β} ln Z)_{V} - \frac{⟨ F ⟩}{T} = \frac{1}{T} (⟨ U ⟩ - ⟨ F ⟩)$ $⟨ P ⟩ = - \frac{1}{β} (\frac{\partial}{\partial V} ln Z)_{T} 、 ⟨ U ⟩ = - (\frac{\partial}{\partial β} ln Z)_{V}$ となり、分配関数から圧力とエネルギーが求まることが分かる。

　一方で、分配関数の方についてはエネルギーが連続していることから $\int_{0}^{\infty} f (U) d U = \frac{1}{Z} \int_{0}^{\infty} W (U) e^{- β U} d U = 1$ と積分の形で書けるものと思われるが、これだと分配関数にエネルギーの次元がついてしまう。そのため、連続している場合 $W = \frac{1}{β} (\frac{\partial W}{\partial U})_{V}$ というような関係が成り立つことを利用して、 $\frac{1}{Z} \int_{0}^{\infty} D (U) e^{- β U} d U = 1 、 D (U) = (\frac{\partial W}{\partial U})_{V}$ というように書ける。ただ、一般的にはハミルトニアンにより $Z = \frac{1}{h ^{3 N}} \int \dots \int e^{- β H} d q_{1} \dots d q_{3 N} d p_{1} \dots d p_{3 N}$ というように書くことが多い（このようにすると分配関数を状態数なしで求めることができる）。このようにすると分布関数も $ρ = \frac{e ^{- β H}}{Z}$ に置き換えた方がよいわけだが、これからエネルギーの期待値を求めてみると $⟨ H ⟩ = - \frac{1}{β} ⟨ ln ρ ⟩ - \frac{1}{β} ln Z \leftrightarrow ⟨ U ⟩ = T ⟨ S ⟩ + ⟨ F ⟩$ であるから、以下のような関係があることが予想される。 $⟨ S ⟩ = - k_{B} ⟨ ln ρ ⟩ 、 ⟨ F ⟩ = - \frac{1}{β} ln Z$ 確認のため、仮に小正準集団の場合を考えるとエネルギーが一定値 $U$ をとることから $ρ$ が $Z = W (U) e^{- β U} = e^{- βF (U)} 、 ρ = \frac{e ^{- β U}}{Z} = \frac{1}{W ( U )}$ であるため、エントロピーとヘルムホルツの自由エネルギーに関して等式が成り立つことが分かる。 $⟨ S ⟩ = - k_{B} ⟨ ln ρ ⟩ = k_{B} ln W (U) 、 ⟨ F ⟩ = - \frac{1}{β} ln Z = F (U)$ 　以上のことを踏まえて、例として自由粒子の場合を考えると分配間関数が $Z = \frac{V ^{N}}{h ^{3 N}} (\frac{2 mπ}{β})^{3 N /2}$ であることから、各期待値は以下の通りになる。 $⟨ P ⟩ = \frac{N k _{B} T}{V} 、 ⟨ U ⟩ = \frac{3}{2} N k_{B} T 、 ⟨ F ⟩ = - \frac{N}{β} ln [\frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2}]$ $⟨ S ⟩ = \frac{1}{T} (⟨ U ⟩ - ⟨ F ⟩) = N k_{B} {\frac{3}{2} + ln [\frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2}]}$

等温定圧集団

　ここまでは体積が一定の場合を考えていたが、実際の物体は圧力が外部と同じになることが多い。そのため、ここでは温度と圧力が変化する状況でどうなるかを見てみることにする。熱平衡になった場合 $S (U_{t} - U, V_{t} - V) = S (U_{t}, V_{t}) - (\frac{\partial S}{\partial U}) U - (\frac{\partial S}{\partial V}) V = S (U_{t}, V_{t}) - (\frac{U + P V}{T})$ であるから、状態数としては以下の通りになる。 $W_{t} (U_{t}, V_{t}) = U, V \sum W (U, V) W_{b} (U_{t} - U, V_{t} - V)$ $1 = U, V \sum f (U, V) 、 f (U, V) = \frac{W ( U , V ) W _{b} ( U t - U , V _{t} - V )}{W _{t} ( U _{t} , V _{t} )}$ $f (U, V) = \frac{W _{b} ( U _{t} , V _{t} )}{W _{t} ( U _{t} , V _{t} )} W (U, V) e^{- β (U + P V)} 、 β = \frac{1}{k _{B} T}$ ここで分配関数のようにエンタルピー $H$ を用いて $Y = \frac{W _{t} ( U _{t} , V _{t} )}{W _{b} ( U _{t} , V _{t} )} = U, V \sum W (U, V) e^{- β H} = 、 H = U + P V$ とすることで、以下の通りにかける。 $f (U, V) = \frac{1}{Y} W (U, V) e^{- β H} = \frac{1}{Y} e^{- βG} 、 G = H - TS$ また、ギブスの自由エネルギー $G$ の微小変化をとると $d G = d H - d TS + T d S = - S d T + V d P$ ここで、エンタルピーの微小変化の式を利用している。 $d H = d U + d P V + P d V = T d S + V d P$

上記の関係式から、以下の関係が成り立つ。 $S = - (\frac{\partial G}{\partial T})_{P} 、 V = (\frac{\partial G}{\partial P})_{T} 、 H = - T^{2} [\frac{\partial}{\partial T} (\frac{G}{T})]_{P}$ そのため、期待値は以下の通りにかける。 $⟨ S ⟩ = - \frac{1}{T} (\frac{\partial}{\partial β} ln Y)_{P} - \frac{⟨ G ⟩}{T} = \frac{1}{T} (⟨ H ⟩ - ⟨ G ⟩)$ $⟨ V ⟩ = - \frac{1}{β} (\frac{\partial}{\partial P} ln Y)_{T} 、 ⟨ H ⟩ = - (\frac{\partial}{\partial β} ln Y)_{P} = ⟨ U ⟩ + P ⟨ V ⟩$ そして、これが連続している場合には $W = \frac{1}{βP} (\frac{\partial W}{\partial V})_{U, N} \leftrightarrow Z \to (\frac{\partial Z}{\partial V})_{U, N}$ であることをふまえると、以下のように書くこともできる。 $ρ_{Y} = \frac{1}{Y} e^{- β (H + P V)} 、 Y = \int_{0}^{V_{t}} (\frac{\partial Z}{\partial V})_{U, N} e^{- βP V} d V 、$ このうち、 $ρ_{Y}$ からエンタルピーの期待値を求めると $⟨ H + P V ⟩ = - \frac{1}{β} ⟨ ln ρ_{Y} ⟩ - \frac{1}{β} ln Y \leftrightarrow ⟨ H ⟩ = T ⟨ S ⟩ + ⟨ G ⟩$ であるから、以下の関係になっていることが予想される（小正準集団においたときも成り立つ）。 $⟨ S ⟩ = - k_{B} ⟨ ln ρ_{Y} ⟩ 、 ⟨ G ⟩ = - \frac{1}{β} ln Y$ 理想気体での分配関数を求めてみると $Y = \frac{N}{h ^{3 N}} (\frac{2 mπ}{β})^{3 N /2} \int_{0}^{V_{t}} V^{N - 1} e^{- βP V} d V$ であるが、 $V_{t}$ が十分大きいとして $V_{t} \to \infty$ とすると $V_{t} \to \infty lim \int_{0}^{V_{t}} V^{N - 1} e^{- βP V} d V = \frac{( N - 1 )!}{( βP ) ^{N}}$ であるため、以下の通りになる。 $Y = \frac{( 2 mπ ) ^{3 N /2}}{h ^{3 N}} \frac{N !}{β ^{5 N /2} P ^{N}}$ となるので、各期待値は次の形で得られる。 $⟨ V ⟩ = \frac{N k _{B} T}{P} 、 ⟨ H ⟩ = \frac{5}{2} N k_{B} T 、 ⟨ U ⟩ = ⟨ H ⟩ - P ⟨ V ⟩ = \frac{3}{2} N k_{B} T$ $⟨ G ⟩ = - \frac{N}{β} ln \frac{( 2 πm ) ^{3/2}}{h ^{3} β ^{5/2} P} - \frac{1}{β} ln N!$ $⟨ S ⟩ = N k_{B} [\frac{3}{2} + ln \frac{⟨ V ⟩}{h ^{3}} (\frac{2 πm}{β})^{3/2}] + k_{B} (N - N ln N + ln N!)$ さらに、これからヘルムホルツの自由エネルギーも求めてみると $⟨ F ⟩ = ⟨ U ⟩ - T ⟨ S ⟩ = - \frac{N}{β} ln \frac{⟨ V ⟩}{h ^{3}} (\frac{2 πm}{β})^{3/2} - k_{B} T (N - N ln N + ln N!)$ というように、前回の形と同様な形が出てくるが後半に余計な項がついてくることが分かる。これについては、分配関数を $N!$ で割る（組み合わせの重複分を消す）ことで解消できる。実際、 $Z \to \frac{Z}{N !}$ とおくと、まず正準集団における分配関数は $Z = \frac{V ^{N}}{N ! h ^{3 N}} (\frac{2 mπ}{β})^{3 N /2}$ となるため、ヘルムホルツの自由エネルギーとエントロピーは以下の通りになる。 $⟨ F ⟩ = - \frac{N}{β} ln [\frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2}] + \frac{1}{β} ln N!$ $⟨ S ⟩ = \frac{1}{T} (⟨ U ⟩ - ⟨ F ⟩) = N k_{B} {\frac{3}{2} + ln [\frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2}]} - k_{B} ln N!$ 一方で、等温定圧集団においては分配関数が $Y = \frac{( 2 mπ ) ^{3 N /2}}{h ^{3 N}} \frac{1}{β ^{5 N /2} P ^{N}}$ であるから、ヘルムホルツの自由エネルギーとエントロピーは以下の通りになる。 $⟨ F ⟩ = ⟨ U ⟩ - T ⟨ S ⟩ = - \frac{N}{β} ln \frac{⟨ V ⟩}{h ^{3}} (\frac{2 πm}{β})^{3/2} + \frac{1}{β} (N ln N - N)$ $⟨ S ⟩ = N k_{B} {\frac{3}{2} + ln [\frac{⟨ V ⟩}{h ^{3}} (\frac{2 mπ}{β})^{3/2}]} - k_{B} (N ln N - N)$ $⟨ G ⟩ = ⟨ H ⟩ - T ⟨ S ⟩ = - \frac{N}{β} ln \frac{( 2 πm ) ^{3/2}}{h ^{3} β ^{5/2} P} = - N k_{B} T ln [\frac{⟨ V ⟩}{N h ^{3}} (2 πm k_{B} T)^{3/2}]$

これらを比較すると、以下の関係が成り立てば互いに等しい関係になっていることが分かる。 $ln N! = N ln N - N$ 実は、この関係はStirlingの公式と呼ばれており、近似的にこの関係が成り立つことが知られている。

化学ポテンシャル

　ここまで統計集団を考えることで、そのエネルギーや圧力などを期待値で求めることができていた。しかし、その中には化学的性質（気体などの構成物質）を考慮されていないため、熱力学第１法則に $d U = d^{'} Q + d^{'} W + μ d N$ というように化学ポテンシャル $μ$ を導入する（一つの物質のみ存在する場合）。このようにすると、粒子数の変化があった場合も取り扱うことが可能になる。各集団における化学ポテンシャルの形は $⎩ ⎨ ⎧ ・小正準集団： μ = - T (\frac{\partial S}{\partial N})_{U, V} ・正準集団： μ = - \frac{1}{β} (\frac{\partial}{\partial N} ln Z)_{T, V} ・等圧定圧集団： μ = - \frac{1}{β} (\frac{\partial}{\partial N} ln Y)_{T, P}$ として求まるので、自由粒子の場合を例にあげて各集団ごとにどうなるかを見ていくことにする。

　小正準集団の場合は、状態数ならびにエントロピーが前項の $N!$ を考慮して $W = \frac{V ^{N}}{N ! h ^{3 N}} \frac{2 π ^{3 N /2} δ p}{Γ ( 3 N /2 )} (2 m U)^{3 N /2 - 1}$ $S = N k_{B} ln \frac{V ( 2 πm U ) ^{3/2}}{h ^{3}} - k_{B} ln N! - k_{B} ln Γ (3 N /2) + ln \frac{2 δ p}{2 m U}$ であるが、Stringの公式を用いると $ln N! ≃ N ln N - N 、 ln Γ (3 N /2) = ln (3 N /2)! ≃ 3 N /2 ln (3 N /2) - 3 N /2$ と置けるため、以下の形になる。 $S = N k_{B} ln [\frac{8 V}{h ^{3}} (\frac{πm U}{3})^{3/2} (\frac{e}{N})^{5/2}] + ln \frac{2 δ p}{2 m U}$ これにより、化学ポテンシャルは以下の通りになる。 $μ = - T (\frac{\partial S}{\partial N})_{U, V} = - k_{B} T ln [\frac{V}{N h ^{3}} (2 πm k_{B} T)^{3/2}]$ 　同様に、正準集団においても分配関数が $Z = \frac{V ^{N}}{N ! h ^{3 N}} (\frac{2 mπ}{β})^{3 N /2} 、 ln Z ≃ N ln [\frac{V e}{N h ^{3}} (\frac{2 πm}{β})^{3/2}]$ であったため、化学ポテンシャルは以下の通りとなる。 $μ = - \frac{1}{β} (\frac{\partial}{\partial N} ln Z)_{T, V} = - k_{B} T ln [\frac{V}{N h ^{3}} (2 πm k_{B} T)^{3/2}]$ 　さらに、等圧定温集団においても分配関数が $Y = \frac{( 2 mπ ) ^{3 N /2}}{h ^{3 N}} \frac{1}{β ^{5 N /2} P ^{N}} 、 ln Y = N ln [\frac{( 2 πm ) ^{3/2}}{h ^{3} β ^{5/2} P}] = N ln [\frac{⟨ V ⟩}{N h ^{3}} (\frac{2 πm}{β})^{3/2}]$ であるから、上記と同じ形が得られることが確認できる。 $μ = - \frac{1}{β} (\frac{\partial}{\partial N} ln Y)_{T, P} = - k_{B} T ln [\frac{⟨ V ⟩}{N h ^{3}} (2 πm k_{B} T)^{3/2}]$ ここで、ギブスの自由エネルギと見比べると $N$ 以外の部分が一致しているため、 $μ = \frac{⟨ G ⟩}{N}$ となり、化学ポテンシャルは一粒子あたりの $G$ であることが予想される。これについては、ギブスの自由エネルギーが $G = U - TS + P V$ であったが、この微小変化をとると $d G = - S d T + V d P + μ d N$ となることから分かるように $G (T, P, N)$ の関数であることからきている。仮に、先ほどの式が成り立つとなると $μ (T, P)$ であることになる。 $G (T, P, N) = μ (T, P) N$ そして、この微小変化をとると $d G = d μ N + μ d N$ であるから、先ほどの式と比較すると $d μ = - \frac{S}{N} d T + \frac{V}{N} d P$ というように確かに $μ (T, P)$ となっていることが分かる。これはGibbs-Duhemの式と呼ばれており、化学ポテンシャルの変数に関することだけではなく、 $T, P, μ$ が互いに依存しており、独立して変化させることができないとを意味している。これが正しいものかは、 $- \frac{S}{N} = (\frac{\partial μ}{\partial T})_{P} 、 \frac{V}{N} = (\frac{\partial μ}{\partial P})_{T}$ であることから、化学ポテンシャルを圧力で表示した $μ = - k_{B} T ln [\frac{( 2 πm ) ^{3/2} ( k _{B} T ) ^{5/2}}{P h ^{3}}]$ を代入すると以前までに得たものと同じ形であることが分かる。 $S = N k_{B} ln [\frac{( 2 πm ) ^{3/2} ( k _{B} T ) ^{5/2}}{P h ^{3}} e^{5/2}] 、 P V = N k_{B} T$ 　　そして、ギブスの自由エネルギーが化学ポテンシャルで書けることから、それに関わる量も $U = TS - P V + μ N 、 F = - P V + μ N 、 H = TS - μ N$ というようになる。

大正準集団

　前回までは等温定圧集団まで扱ったが，それに加えて粒子数も変化することを考える。しかし、これを考えた時には、エントロピーの近似式が $S (U_{t} - U, V_{t} - V, N_{t} - N) = S (U_{t}, V_{t}, N_{t}) - (\frac{U + P V - μ N}{T})$ $U = TS - P V + μ N$ となるため、全体の状態数と熱浴の状態数が同じものとなってしまうことが分かる（完全に開放した状態になるため、完全に同化したともいえる）。ただ、実際のところ等温低圧集団よりも正準集団の方が計算量が少なくほぼ同じ結果が導けていたので、正準集団において粒子数が変化する集団を考えることにする。このときは、 $S (U_{t} - U, N_{t} - N) = S (U_{t}, N_{t}) - (\frac{U - μ N}{T})$ であることから、分布関数は以下の通りになる（ $J : グランドポテンシャル$ ）。 $f (U, N) = \frac{1}{Ξ} W (U, N) e^{- β (U - μ N)} = \frac{1}{Ξ} e^{- β J}$ $J = U - TS - μ N = - P V 、 Ξ = U, N \sum W (U, N) e^{- β (U - μ N)}$ グランドポテンシャルの微小変化は $d J = - P d V - S d T - N d μ$ であることから、 $(\frac{\partial J}{\partial V})_{T, μ} = - P 、 (\frac{\partial J}{\partial T})_{V, μ} = - S 、 (\frac{\partial J}{\partial μ})_{V, T} = - N$ が成り立つため、各気体値は以下のように書ける。 $⟨ P ⟩ = \frac{1}{β} (\frac{\partial}{\partial V} ln Ξ)_{T, μ} 、 ⟨ N ⟩ = \frac{1}{β} (\frac{\partial}{\partial μ} ln Ξ)_{V, T}$ $⟨ S ⟩ = - \frac{⟨ J ⟩}{T} + \frac{1}{β} (\frac{\partial}{\partial T} ln Ξ)_{V, μ} = k_{B} (\frac{\partial}{\partial T} T ln Ξ)_{V, μ}$ これを連続的な形で記載した場合には、 $ρ = \frac{1}{Ξ} e^{- β (H - μ N)} 、 Ξ = N = 0 \sum N_{t} Z e^{β μ N}$ というようになることに踏まえて、 $⟨ H ⟩ - μ ⟨ N ⟩ = - \frac{1}{β} ⟨ ln ρ ⟩ - \frac{1}{β} ln Ξ$ であるため、以下の関係が成り立つことが予想される。 $⟨ S ⟩ = - k_{B} ⟨ ln ρ ⟩ 、 ⟨ J ⟩ = - \frac{1}{β} ln Ξ$ 実際、小正準集団とみなした場合には、 $⟨ S ⟩ = k_{B} ln W (U, N) 、 ⟨ J ⟩ = J (U, N)$ となる。さらに、自由粒子の場合で解いてみると、 $Ξ = N = 0 \sum N_{t} \frac{V ^{N}}{N ! h ^{3 N}} (\frac{2 mπ}{β})^{3 N /2} e^{β μ N} = N = 0 \sum N_{t} \frac{1}{N !} [\frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ}]^{N}$ となるため、総和に関して $N = 0 \sum \infty \frac{x ^{N}}{N !} = e^{x}$ であることを利用して、 $N_{t} \to \infty$ とすると $Ξ = exp [\frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ}] 、 ln Ξ = \frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ}$ となるため、各期待値は以下の通りになる。 $⟨ P ⟩ = \frac{1}{β h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ} 、 ⟨ N ⟩ = \frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ} 、 ⟨ J ⟩ = - \frac{V}{β h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ}$ $⟨ S ⟩ = k_{B} \frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ} (\frac{5}{2} - β μ)$ $⟨ U ⟩ = T ⟨ S ⟩ - ⟨ P ⟩ V + μ ⟨ N ⟩ = \frac{3}{2} k_{B} T \frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ}$ あまりにも見にくいので、改めて整理すると以下のようになる。 $⟨ N ⟩ = \frac{V}{h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{β μ} 、 μ = - \frac{1}{β} ln [\frac{V}{⟨ N ⟩ h ^{3}} (\frac{2 mπ}{β})^{3/2}]$ $⟨ P ⟩ = \frac{⟨ N ⟩}{β V} 、 ⟨ J ⟩ = - ⟨ P ⟩ V 、 ⟨ U ⟩ = \frac{3}{2} ⟨ N ⟩ k_{B} T$ $⟨ S ⟩ = ⟨ N ⟩ k_{B} (\frac{5}{2} - β μ) = ⟨ N ⟩ k_{B} ln [\frac{V}{⟨ N ⟩ h ^{3}} (\frac{2 mπ}{β})^{3/2} e^{5/2}]$ このように、従来の結果と同じ結果が得られることが分かる。

相対性理論

　晴れた日に空を見上げると、朝から夕方にかけて太陽が東側から西側へと移動し、夜になると幾つもの星が流れていく様子を見ることができます。一見すると、全てのものが自分（地球）を中心に回っているように思われます。実際、2世紀ごろアレクサンドリアにいたプトレマイオスが\tb{天動説}を唱え、15世紀になるまで信じられてきました。ところが、それ以降にNicolaus\ CopernicusやGalileo\ Galileiなどが独自の考察や天体の観測をもとに\tb{地動説}（地球の方が太陽を中心として周囲を回っているという説）を提唱したことで、その状況が変わり始めます。当時としては、宗教的にも天動説が信じられていたので、宗教裁判なるものまでかけられるまでの大事になりました。しかし、ここで互いに見たときの運動の様子は違うけれども、運動の法則自体は同じものであるということにもなります。

　しかし、慣性の法則は”絶対空間”に対して一様な速度で運動しているすべてのゆがむことのない準拠系に対しても成り立つ。それは自由な粒子は、このような座標系に対しても一様な並進運動をするからである。慣性の法則が成り立つような準拠系が全て慣性系という。慣性系全体は互いに他に対して等速直線運動をする歪まない準拠系の三重に無限の集合である。これらのうちの１つが恒星全体に対して静止している絶対系であるけれども、慣性の法則が成り立つという点では、どの慣性系も完全に同等である。

　さて、力学における相対性原理というのは、すべての慣性系は、慣性の法則に対してのみならず他のあらゆる力学法則に対しても完全に同等であることを述べたものである。もしこれが正しいならば、力学現象はどの観測系においても同じ推移の経過をたどることになるから、このような現象を観察することによって、この系全体が”絶対空間”に対して行っている一様な運動を検知することは不可能である。したがって力学現象を調べるだけで”絶対”系を決定することはできない。

　次に、Newton力学の基礎方程式が実際にこの相対性原理にかなっていることを示そう。

同時の定義

　物体の運動を扱うとき、その位置（座標）と発生した時刻を観測する必要がある。各物理量は観測者のいる地点を原点として、座標は原点からの距離、時間は原点にて観測者が時計を用いて測定を行う。ただ原点から離れた地点で発生した事象を観測した場合だと、観測者の目に光が届くまでに時間を要する場合も考えられるため、厳密には事象が発生した地点でなければ同時に測定できないものと思われる（座標についても同時であることを確認しなければならない）。さらにそれらの時計が原点の時計と常に同時である（同期している）かを確認しなければならない。同期した状態で時計を持っていくにしても、動かすことで観測者と異なる状態になるため時間がずれる可能性があるため、動かして置いた後で再度時計を合わせる必要がある。そこで、一定の速度で移動する信号を往復させることで確認を行う。

$t_{C} - t_{B} = t_{B} - t_{A}$

この信号としては光がある（真空中でも伝わる）。さらに過去の実験結果を踏まえて、観測者ごとに光の速度が同じように見える（光速度不変の原理）ものとして観測者間の関係を見てみよう。

　まず、下図のように観測者Kと観測者K'が同じ位置にいてKから見てK'が時刻0にて速さ $V$ で $x$ 軸に沿って動いたとする。同時に光を原点から放出して地点Bにて折り返す。このとき、各観測者から見てどのような式が成り立つか見てみると以下のようになる。

[c (t_{C} - t_{B})]^{2} = (c t_{B})^{2} + (V t_{C})^{2} - 2 c V t_{B} t_{C} (\frac{x _{B}}{c t _{B}})

x_{B}^{2} + y_{B}^{2} + z_{B}^{2} = (c t_{B})^{2} 、 x_{B}^{'2} + y_{B}^{'2} + z_{B}^{'2} = (c t_{B}^{'})^{2}

x_{C} - V t_{C} = y_{C} = z_{C} = 0 、 x_{C}^{'} = y_{C}^{'} = z_{C}^{'} = 0

t_{C} = 2 \frac{c ^{2}}{c ^{2} - V ^{2}} (t_{B} - \frac{V}{c ^{2}} x_{B}) 、 t_{C}^{'} = 2 t_{B}^{'}

ここでC点において各観測者の間で $x_{C}^{'} \Leftrightarrow x_{C} - V t_{C} 、 y_{C}^{'} \Leftrightarrow y_{C} 、 z_{C}^{'} \Leftrightarrow z_{C} 、 t_{C}^{'} \Leftrightarrow t_{C} 、 t_{B}^{'} \Leftrightarrow \frac{c ^{2}}{c ^{2} - V ^{2}} (t_{B} - \frac{V}{c ^{2}} x_{B})$ という対応をしているものとして以下のように置いてみる。 $x^{'} = x - V t 、 y^{'} = y 、 z^{'} = z 、 t^{'} = \frac{c ^{2}}{c ^{2} - V ^{2}} (t - \frac{V}{c ^{2}} x)$ しかし、B点での関係式は満たさないという問題が起きる。 $x^{'2} - (c t^{'})^{2} \neq = x^{2} - (c t)^{2}$ ここで、この式以外で成立していることと上記が二次の方程式であることから $x^{'} = a x + b t 、 t^{'} = d x + e t （ a 、 b 、 d 、 e ：定数）$ というようにおいてみて、その他に成立する形がないか確認してみることにする。単純に代入して計算してみると $a^{2} - c^{2} d^{2} = 1 、 ab - c^{2} d e = 0 、 c^{2} e^{2} - b^{2} = c^{2}$ というようになるが、ここで条件として各観測者の原点で $x^{'} (x = V t) = 0$ という関係になっていることから $b = - aV$ であることから以下の通りになる。 $a^{2} - c^{2} d^{2} = 1 、 - a^{2} V - c^{2} d e = 0 、 c^{2} e^{2} - a^{2} V^{2} = c^{2}$ そして、このうちの2式において $a^{2} = 1 + c^{2} d^{2} 、 e = - \frac{V a ^{2}}{c ^{2} d}$ というようにすると、これらを代入することで以下の式が得られる。 $d^{2} = \frac{V ^{2}}{c ^{2} ( c ^{2} - V ^{2} )}$ すると、これにつられて他の定数も以下の形にかける。 $a^{2} = \frac{c ^{2}}{c ^{2} - V ^{2}}$ ここで $V$ が0となったときに、 $x^{'}$ と $x$ が一致していないといけないことから、 $a = \frac{c}{c ^{2} - V ^{2}}$ となる。一方で、 $d$ と $e$ の方は $d = \pm \frac{V}{c c ^{2} - V ^{2}} 、 e = \mp \frac{c}{c ^{2} - V ^{2}}$ となるが、こちらも $V$ が0となったときに、 $t^{'}$ と $t$ が一致しないといけないことから、 $e = \frac{c}{c ^{2} - V ^{2}} 、 d = - \frac{V}{c c ^{2} - V ^{2}}$ となる。以上のことから変換式を改めて記述すると以下の通りになる。 $x^{'} = γ (x - V t) 、 y^{'} = y 、 z^{'} = z 、 t^{'} = γ (t - \frac{V}{c ^{2}} x) (γ = \frac{1}{1 - V ^{2} / c ^{2}})$ このようにすると全ての関係式を満たすことが分かるが、ここで一つ気になるものとして $t_{C}^{'} = 1 - \frac{V ^{2}}{c ^{2}} t_{C}$ が出てくるが、これは光速度不変の原理により時間と座標の捉え方が変わったことで出てきたものと考えられる。すなわち、観測者KがK'と同じ視点に立つとある一点に向けて進み反射して戻ってくるのが往復するだけになるため、光速度が不変となるには座標と時間が短くなる必要があるということになる。

Lorentz変換

先ほど導いた変換式を改めて記載してみると $x^{'} = γ (x - V t) 、 y^{'} = y 、 z^{'} = z 、 t^{'} = γ (t - \frac{V}{c ^{2}} x) (γ = \frac{1}{1 - V ^{2} / c ^{2}})$ となるわけだが、これらが何を意味しているのか考えてみることにする。従来から変換式としてはGalilei変換が利用されており $x^{'} = x - V t 、 y^{'} = y 、 z^{'} = z 、 t^{'} = t$ という形であるが、 $V / c \to 0$ と観測者が光速に比べて低速で動いている場合は近似的にこの形になることが分かる。そのため、変換式を以下のように変形してみる。 $x^{'} = γ (x - β c t) 、 y^{'} = y 、 z^{'} = z 、 c t^{'} = γ (c t - β x) (γ = \frac{1}{1 - β ^{2}} 、 β = \frac{V}{c})$ すると、 $c t$ が $x$ と同じような座標の一部としてみなせるので $w$ として $x^{'} = γ (x - βw) 、 w^{'} = γ (w - β x)$ とおき、試しに $x$ と $w$ との関係を座標で描いてみると以下の図の通りになる。

この形を見ると斜交座標の形をしていることから仮に $x$ 軸と $x^{'}$ 軸あるいは $w^{'}$ 軸とのなす角をそれぞれ $θ$ 、 $ϕ$ とすると $w = w^{'} sin ϕ + x^{'} sin θ 、 x = w^{'} cos ϕ + x^{'} cos θ$ となるため、Lorentz変換を逆変換したもの $w = γ (w^{'} + β x^{'}) 、 x = γ (x^{'} + β w^{'})$ と比較すると以下の関係が成り立つことが予想される。 $sin ϕ = γ 、 sin θ = γ β 、 cos ϕ = γ β 、 cos θ = γ$ しかし、これでは三角関数の公式を満たさない。 $sin^{2} θ + cos^{2} θ = \frac{β ^{2} + 1}{1 - β ^{2}} \neq = 1 、 sin^{2} ϕ + cos^{2} ϕ = \frac{1 + β ^{2}}{1 - β ^{2}} \neq = 1$

そこで分子と分母の符号が反対であることを踏まると、双曲線関数であれば等式が成り立つことが分かる。実際、先ほどの三角関数と同じように $sinh ϕ = γ 、 sinh θ = γ β 、 cosh ϕ = γ β 、 cosh θ = γ$ としてみると、以下の関係式が成り立つことが分かる。 $cosh^{2} θ - sinh^{2} θ = 1 、 cosh^{2} ϕ - sinh^{2} ϕ = - 1$ そのため、Lonrentz変換は以下のように書けることになる。 $w = w^{'} sinh ϕ + x^{'} sinh θ 、 (sinh ϕ = γ 、 sinh θ = γ β)$ $x = w^{'} cosh ϕ + x^{'} cosh θ 、 (cosh ϕ = γ β 、 cosh θ = γ)$ この変換自体は以下の図のように双曲線に沿って回転するものとなっており、通常の回転とは異なっていることが分かる。一例として $w$ 軸が回転することで点線（漸近線）に近づいていき、やがて $w^{'}$ 軸は $w = x$ の直線と一致する。このとき、角度 $ϕ$ に関しては $ϕ \to \infty$ であることから, $β \to 1 (V \to c)$ というように観測者の速度が光速を上限とした値になっていると考えられる（ $θ$ も $β$ に依存するため同じようになっているといえる）。

　以上が変換自体の話であるが、前回観測者ごとに光速度が不変となるには座標と時間が短くなる必要があるということを述べていたことや最後のところででてきた式 $t_{C}^{'} = 1 - \frac{V ^{2}}{c ^{2}} t_{C}$ も説明なしに出てきたため、まずはある地点での時間に関する変換を見ていくことにする。

時間の遅れ

　観測者K、K'がおり、Kから見てK'がx軸に速度Vで移動しているものとする。このとき、K'においてある一つの地点Pでの時間を計測したとき

$x_{1}^{'} = x_{2}^{'} 、 t_{2}^{'} - t_{1}^{'} = Δ t^{'}$ が成立する。一方で、Lotentz変換に基づいてKからK'の立場で考えてみると $x_{2}^{'} - x_{1}^{'} = γ [(x_{2} - x_{1}) - V (t_{2} - t_{1})] = 0$ $t_{2}^{'} - t_{1}^{'} = γ [(t_{2} - t_{1}) - \frac{V}{c ^{2}} (x_{2} - x_{1})] = Δ t^{'}$ となる。ただし、 $γ = (1 - V^{2} / c^{2})^{- 1/2}$ である。すると、以下の関係式が得られる。 $(x_{2} - x_{1}) = V (t_{2} - t_{1}) = V Δ t$ $Δ t^{'} = 1 - \frac{V ^{2}}{c ^{2}} Δ t < Δ t$ 一つ目の式は点Pが動いていることから分かるが、2つ目の式については違和感を感じであろう。 $Δ t^{'}$ の方はK'の立場から見て静止している地点（同じ地点で）の時間経過を表しているのに対し、 $Δ t$ は別の地点で測定した時間を表している。つまり、動いているものの立場K'で見るとその速度に応じて時間が遅れているように見えることになる。先ほど出てきた式も時刻0から時間を計測しており $Δ t = t_{C}$ となるため、この事象によるものとなっている。

Lorentz収縮

　次に2点ABで上記と同様に時間を計測すると、その2点間において同時であることから $x_{B}^{'} - x_{A}^{'} = L^{'} 、 t_{A}^{'} = t_{B}^{'}$ が成立する。一方で、Lotentz変換に基づいてKからK'の立場で考えてみると $x_{B}^{'} - x_{A}^{'} = γ [(x_{B} - x_{A}) - V (t_{B} - t_{A})] = L^{'}$ $t_{B}^{'} - t_{A}^{'} = γ [(t_{B} - t_{A}) - \frac{V}{c ^{2}} (x_{B} - x_{A})] = 0$ となり、以下の関係式が得られる。 $V (t_{B} - t_{A}) = \frac{V ^{2}}{c ^{2}} (x_{B} - x_{A}) = \frac{V ^{2}}{c ^{2}} L = β^{2} L$ $L^{'} = 1 - \frac{V ^{2}}{c ^{2}} L < L$ このことから、まず片方で同時に測ったとしてももう片方が同じ立場で見たときはそうではなく $β^{2} L$ だけ移動しており、その分を $L$ から引いて $γ$ をかけることでKはK'と同じ立場で見れるが長さが短く見えてしまうという結果が得られる。

相対論的力学

　Lotentz変換から速度の式を求めてみると $v_{x}^{'} = \frac{v _{x} - V}{1 - \frac{v _{x} V}{c ^{2}}} 、 v_{y}^{'} = \frac{v _{y}}{γ ( 1 - \frac{v _{x} V}{c ^{2}} )} 、 v_{z}^{'} = \frac{v _{z}}{γ ( 1 - \frac{v _{x} V}{c ^{2}} )} 、 v_{w}^{'} = c$ となるが、ここで $v_{x}^{'}$ に着目すると $v_{x}$ に対して以下の条件下でないと発散することが分かる。 $v_{x} < \frac{c ^{2}}{V}$ 一方で $V$ に関しても、こちらも $c$ をこえると $γ (1 - V^{2} / c^{2})$ が虚数となってしまうことから以下の条件になることも分かる。 $V < c 、 v_{x} < c$ このことから、物体の速度というのは光速を超えることがなく、力が働いて加速しても光速で頭打ちになることがうかがえる。これは別の見方だと、速度が増すことで加速しにくくなる（質量が増す）ことになるため、質量が速度に依存しているともいえる。古典力学においても質量が速度（というより時間）に依存している場合は以下の方程式を解くことで運動を記述することができた。ここで $p$ は運動量、 $F$ は外部からの力である。 $F = \frac{d p}{d t} 、 p = m v$ 今回もこれと同様の式が成立しているものとして、質量がどのように速度に依存しているか確認してみることにする。これについてはすでにN.LwisとC.Tolmanの論文（The Principle of Relativity and Non-Newtonian Mechanics）において示されており、これでは外部の力がない状態において運動量が保存されることを利用して導いている。これによると $m (v) = \frac{m _{0}}{1 - \frac{v ^{2}}{c ^{2}}}$ というようになるとされている。ここで $m_{0}$ は物体が静止しているときの質量としている。では、観測者ごとにどうなるかを見ていくと、まず質量に関して $m (v^{'}) = \frac{m _{0}}{1 - \frac{v ^{'} ^{2}}{c ^{2}}}$ であり、速度に関しては $v^{'2} = v_{x}^{'2} + v_{y}^{'2} + v_{z}^{'2} = \frac{1}{( 1 - \frac{v _{x} V}{c ^{2}} ) ^{2}} [(v_{x} - V)^{2} + (1 - \frac{V ^{2}}{c ^{2}}) (v_{y}^{2} + v_{z}^{2})]$ となることから、以下のようになる。 $m (v^{'}) = γ (1 - \frac{v _{x} V}{c ^{2}}) \frac{m _{0}}{1 - \frac{v ^{2}}{c ^{2}}} = γ (1 - \frac{v _{x} V}{c ^{2}}) m (v) = \frac{d t ^{'}}{d t} m (v)$ ただし、ここで以下の関係式を利用している。 $\frac{d t ^{'}}{d t} = γ (1 - \frac{v _{x} V}{c ^{2}})$ そのため、運動方程式の左辺については以下の通りとなる。 $\frac{d}{d t ^{'}} (m (v^{'}) v_{x}^{'}) = γ \frac{d t}{d t ^{'}} [\frac{d}{d t} (m (v) v_{x}) - \frac{d}{d t} (m (v) V)] = \frac{d t}{d t ^{'}} γ (F_{x} - \frac{V}{c} F_{w})$ $\frac{d}{d t ^{'}} (m (v^{'}) v_{y}^{'}) = \frac{d t}{d t ^{'}} \frac{d}{d t} (m (v) v_{y}) = \frac{d t}{d t ^{'}} F_{y}$ $\frac{d}{d t ^{'}} (m (v^{'}) v_{z}^{'}) = \frac{d t}{d t ^{'}} \frac{d}{d t} (m (v) v_{z}) = \frac{d t}{d t ^{'}} F_{z}$ $\frac{d}{d t ^{'}} (m (v^{'}) v_{w}^{'}) = \frac{d t}{d t ^{'}} [\frac{d}{d t} (m (v) c) - \frac{V}{c} \frac{d}{d t} (m (v) v_{x})] = \frac{d t}{d t ^{'}} γ (F_{w} - \frac{V}{c} F_{x})$ ここで、右辺は力に等しいとしているため $F_{x}^{'} d t^{'} = γ (F_{x} - \frac{V}{c} F_{w}) d t 、 F_{y}^{'} d t^{'} = F_{y} d t 、 F_{z}^{'} d t^{'} = F_{z} d t 、 F_{w}^{'} d t^{'} = γ (F_{w} - \frac{V}{c} F_{x}) d t$ となり、力積（運動量の変化）が等しい形になる。そのため $d p^{'} = d p$ であるが、このままだと時間の部分が抜けた形になってしまっていることが分かる。一方で、質量の関係式をもう少し展開していくと $\frac{m _{0}}{1 - \frac{v ^{2}}{c ^{2}}} = \frac{d t}{d t ^{'}} \frac{m _{0}}{1 - \frac{v ^{'} ^{2}}{c ^{2}}}$ であるから、時間に関して以下の等式が成り立つ。 $d t^{'} 1 - \frac{v ^{'} ^{2}}{c ^{2}} = d t 1 - \frac{v ^{2}}{c ^{2}}$ こちらの量も変換に対して不変な形をしており、時間の代わりに扱うと同じ形で置き換えられることが分かる。そこで、 $d τ \equiv d t 1 - \frac{v ^{2}}{c ^{2}}$ として、運動方程式を以下のようにしてみよう。 $f = \frac{d p}{d τ} 、 f \equiv γ (v) F$ こちらも同じように見ていくと、質量に関しては同じとなる一方で各成分についてはそれぞれ $\frac{d}{d τ ^{'}} (m (v^{'}) v_{x}^{'}) = γ [\frac{d}{d τ} (m (v) v_{x}) - \frac{d}{d τ} (m (v) V)] = γ (f_{x} - \frac{V}{c} f_{w})$ $\frac{d}{d τ ^{'}} (m (v^{'}) v_{y}^{'}) = \frac{d}{d τ} (m (v) v_{y}) = f_{y}$ $\frac{d}{d τ ^{'}} (m (v^{'}) v_{z}^{'}) = \frac{d}{d τ} (m (v) v_{z}) = f_{z}$ $\frac{d}{d τ ^{'}} (m (v^{'}) v_{w}^{'}) = γ [\frac{d}{d τ} (m (v) c) - \frac{V}{c} \frac{d}{d τ} (m (v) v_{x})] = γ (f_{w} - \frac{V}{c} f_{x})$ となる。このことから、下記のように力に関して対応していることが分かる。 $f_{x}^{'} = γ (f_{x} - \frac{V}{c} f_{w}) 、 f_{y}^{'} = f_{y} 、 f_{z}^{'} = f_{z} 、 f_{w}^{'} = γ (f_{w} - \frac{V}{c} f_{x})$ また、運動量に対してエネルギーも保存するため、まずエネルギー（運動エネルギー）がどのような形になるか確認するため仕事率（力と速度の内積）を求めてみると $f \cdot v = \frac{d}{d τ} \frac{m _{0} v}{1 - \frac{v ^{2}}{c ^{2}}} \cdot v = \frac{d}{d τ} \frac{m _{0} c ^{2}}{1 - \frac{v ^{2}}{c ^{2}}} = \frac{d}{d τ} (m (v) c^{2})$ となるため、エネルギーは以下の形になることが分かる。 $E_{K} = m (v) c^{2} = m_{0} c^{2} + \frac{1}{2} m_{0} v^{2} + \dots$ このように、 $v$ が $c$ に比べて低速の場合は静止エネルギーと古典力学での運動エネルギーの和になることがわかる。また、仕事率の形を見ると $f_{w}$ とほぼ同じような形をしており以下の関係がある。 $f_{w} = \frac{1}{c} (f \cdot v) = \frac{1}{c} \frac{d E _{K}}{d τ}$

以上を踏まえて、運動方程式と変換式を整理すると以下のようになる。 $f = \frac{d p}{d τ} 、 (f \cdot v) = \frac{d E _{K}}{d τ} 、 p = m (v) v 、 E_{K} = m (v) c^{2}$ $f_{x}^{'} = γ [f_{x} - \frac{V}{c ^{2}} (f \cdot v)] 、 f_{y}^{'} = f_{y} 、 f_{z}^{'} = f_{z} \to (f^{'} \cdot v^{'}) = γ [(f \cdot v) - V f_{x}]$ これらの式に基づいて、実際に重力がかかる場合の運動を考えてみる。重力には質量が関わってくるが、この質量というのは重力質量と呼ばれ、実質的には違うものとなっている。先ほどの質量は慣性（運動の状態を維持しようとする性質）がどのくらいかを表す慣性質量（加速のしにくさ）に相当している（詳細は、また一般相対性理論の項で述べる）。そのため、運動方程式は以下の通りになる。 $\frac{d}{d t} \frac{m _{0} v}{1 - \frac{v ^{2}}{c}} = m_{0} g 、 \frac{d}{d t} \frac{m _{0} c ^{2}}{1 - \frac{v ^{2}}{c}} = m_{0} (g \cdot v)$ ここで、静止した状態から落とした状況を考えると、z成分しかないため $\frac{m _{0} v _{z}}{1 - \frac{v _{z}^{2}}{c ^{2}}} = - m_{0} g t 、 E = m_{0} c^{2} - m_{0} g \cdot (z - z_{0})$ となる。ここで、一つ目の式を整理すると速度と位置が求められる。 $v_{z} = - \frac{g t}{1 + ( \frac{g t}{c} ) ^{2}} 、 z = z_{0} + \frac{c ^{2}}{g} 1 - 1 + (\frac{g t}{c})^{2}$ このとき、光速に比べて $g t$ が小さいものとすると $v_{z} = - g t 、 z = z_{0} - \frac{1}{2} g t^{2}$ というように従来と同じ形になることが分かる。次に電場や磁場が働くとどうなるかを見るため、電磁気学の方を見てみることにする。

相対論的電磁気学

電磁気学においては、以下のMaxwell方程式を用いて説明することができた。 $div D = ρ 、 div B = 0 、 rot E = - \frac{\partial B}{\partial t} 、 rot H = i + \frac{\partial D}{\partial t}$

このうち、２番目と３番目の式に関して $\frac{\partial B _{x}}{\partial x} + \frac{\partial B _{y}}{\partial y} + \frac{\partial B _{z}}{\partial z} = 0$ $\frac{\partial E _{z}}{\partial y} - \frac{\partial E _{y}}{\partial z} = - \frac{\partial B _{x}}{\partial t} 、 \frac{\partial E _{x}}{\partial z} - \frac{\partial E _{z}}{\partial x} = - \frac{\partial B _{y}}{\partial t} 、 \frac{\partial E _{y}}{\partial x} - \frac{\partial E _{x}}{\partial y} = - \frac{\partial B _{z}}{\partial t}$ と各成分ごとに書き下せるが、ここでLotrentz変換を導入すると $\frac{\partial}{\partial x} = γ \frac{\partial}{\partial x ^{'}} - γ \frac{V}{c ^{2}} \frac{\partial}{\partial t ^{'}} 、 \frac{\partial}{\partial y} = \frac{\partial}{\partial y ^{'}} 、 \frac{\partial}{\partial z} = \frac{\partial}{\partial z ^{'}} 、 \frac{\partial}{\partial t} = - γV \frac{\partial}{\partial x ^{'}} + γ \frac{\partial}{\partial t ^{'}}$ となるため、これらを適用すると以下の形が得られる。 $γ \frac{\partial B _{x}}{\partial x ^{'}} + \frac{\partial B _{y}}{\partial y ^{'}} + \frac{\partial B _{z}}{\partial z ^{'}} = γ \frac{V}{c ^{2}} \frac{\partial B _{x}}{\partial t ^{'}}$ $- γV \frac{\partial B _{x}}{\partial x ^{'}} + \frac{\partial E _{z}}{\partial y ^{'}} - \frac{\partial E _{y}}{\partial z ^{'}} = - γ \frac{\partial B _{x}}{\partial t ^{'}}$ $\frac{\partial E _{x}}{\partial z ^{'}} - \frac{\partial}{\partial x ^{'}} γ (E_{z} + V B_{y}) = - \frac{\partial}{\partial t ^{'}} γ (B_{y} + \frac{V}{c ^{2}} E_{z})$ $\frac{\partial}{\partial x ^{'}} γ (E_{y} - V B_{z}) - \frac{\partial E _{x}}{\partial y ^{'}} = - \frac{\partial}{\partial t ^{'}} γ (B_{z} - \frac{V}{c ^{2}} E_{y})$ であるが、ここで1番目の式と2番目の式を $- γV \frac{\partial B _{x}}{\partial x ^{'}} - V \frac{\partial B _{y}}{\partial y ^{'}} - V \frac{\partial B _{z}}{\partial z ^{'}} = - γ \frac{V ^{2}}{c ^{2}} \frac{\partial B _{x}}{\partial t ^{'}}$ $γV \frac{\partial B _{x}}{\partial x ^{'}} - \frac{\partial E _{z}}{\partial y ^{'}} + \frac{\partial E _{y}}{\partial z ^{'}} = γ \frac{\partial B _{x}}{\partial t ^{'}}$ という形にして一番目の式に2番目の式に代入することで $\frac{\partial B _{x}}{\partial x ^{'}} + \frac{\partial}{\partial y ^{'}} γ (B_{y} + \frac{V}{c ^{2}} E_{z}) + \frac{\partial}{\partial z ^{'}} γ (B_{z} - \frac{V}{c ^{2}} E_{y}) = 0$ となる。または足し合わせると以下のようになる。 $\frac{\partial}{\partial y ^{'}} γ (E_{z} + V B_{y}) - \frac{\partial}{\partial z ^{'}} γ (E_{y} - V B_{z}) = - \frac{\partial B _{x}}{\partial t ^{'}}$ 以上のことをひっくるめて、改めて整理した式を書くと以下の通りとなる。 $\frac{\partial B _{x}}{\partial x ^{'}} + \frac{\partial}{\partial y ^{'}} γ (B_{y} + \frac{V}{c ^{2}} E_{z}) + \frac{\partial}{\partial z ^{'}} γ (B_{z} - \frac{V}{c ^{2}} E_{y}) = 0$ $\frac{\partial}{\partial y ^{'}} γ (E_{z} + V B_{y}) - \frac{\partial}{\partial z ^{'}} γ (E_{y} - V B_{z}) = - \frac{\partial B _{x}}{\partial t ^{'}}$ $\frac{\partial E _{x}}{\partial z ^{'}} - \frac{\partial}{\partial x ^{'}} γ (E_{z} + V B_{y}) = - \frac{\partial}{\partial t ^{'}} γ (B_{y} + \frac{V}{c ^{2}} E_{z})$ $\frac{\partial}{\partial x ^{'}} γ (E_{y} - V B_{z}) - \frac{\partial E _{x}}{\partial y ^{'}} = - \frac{\partial}{\partial t ^{'}} γ (B_{z} - \frac{V}{c ^{2}} E_{y})$ これらの式から、以下の関係式が成り立つことがわかる。 $E_{x^{'}} = E_{x} 、 E_{y^{'}} = γ (E_{y} - V B_{z}) 、 E_{z^{'}} = γ (E_{z} + V B_{y})$ $B_{x^{'}} = B_{x} 、 B_{y^{'}} = γ (B_{y} + \frac{V}{c ^{2}} E_{z}) 、 B_{z^{'}} = γ (B_{z} - \frac{V}{c ^{2}} E_{y})$ 確認の例として、電荷 $q$ の点電荷に電場 $E$ ' と磁場 $B^{'}$ が働いているとすると $\frac{d p ^{'}}{d τ ^{'}} = f = \frac{1}{1 - \frac{v ^{'2}}{c ^{2}}} [q E^{'} + q (v^{'} \times B^{'})]$ $\frac{d E ^{'}}{d τ ^{'}} = \frac{1}{1 - \frac{v ^{'2}}{c ^{2}}} [q E^{'} + q (v^{'} \times B^{'})] \cdot v^{'} = \frac{1}{1 - \frac{v ^{'2}}{c ^{2}}} (v^{'} \cdot q E^{'})$ となるが、各成分に対して変換すると以下の通りになる。 $\frac{d p _{x}^{'}}{d τ ^{'}} = \frac{1}{1 - \frac{v ^{2}}{c ^{2}}} γ [q E_{x} + q (v_{y} B_{z} - v_{z} B_{y}) - \frac{V}{c ^{2}} (v \cdot q E)]$ $\frac{d p _{y}^{'}}{d τ ^{'}} = \frac{1}{1 - \frac{v ^{2}}{c ^{2}}} [q E_{y} + q (v_{z} B_{x} - v_{x} B_{z})]$ $\frac{d p _{z}^{'}}{d τ ^{'}} = \frac{1}{1 - \frac{v ^{2}}{c ^{2}}} [q E_{z} + q (v_{x} B_{y} - v_{y} B_{x})]$ $\frac{d E ^{'}}{d τ ^{'}} = \frac{1}{1 - \frac{v ^{2}}{c ^{2}}} [(v \cdot q E) - V {E_{x} + q (v_{y} B_{z} - v_{z} B_{y})}]$ これにより、力に関して以下の対応になっていることが再確認できる。 $f_{x}^{'} = \frac{1}{1 - \frac{v ^{'2}}{c ^{2}}} [q E_{x}^{'} + q (v_{y}^{'} B_{z}^{'} - v_{z}^{'} B_{y}^{'})] \leftrightarrow f_{x} = \frac{1}{1 - \frac{v ^{2}}{c ^{2}}} [q E_{x} + q (v_{y} B_{z} - v_{z} B_{y})]$

$f_{y}^{'} = \frac{1}{1 - \frac{v ^{'2}}{c ^{2}}} [q E_{y}^{'} + q (v_{z}^{'} B_{x}^{'} - v_{x}^{'} B_{z}^{'})] \leftrightarrow f_{y} = \frac{1}{1 - \frac{v ^{2}}{c ^{2}}} [q E_{y} + q (v_{z} B_{x} - v_{x} B_{z})]$ $f_{z}^{'} = \frac{1}{1 - \frac{v ^{'2}}{c ^{2}}} [q E_{z}^{'} + q (v_{x}^{'} B_{y}^{'} - v_{y}^{'} B_{x}^{'})] \leftrightarrow f_{z} = \frac{1}{1 - \frac{v ^{2}}{c ^{2}}} [q E_{z} + q (v_{x} B_{y} - v_{y} B_{x})]$ $f_{w}^{'} = \frac{1}{c 1 - \frac{v ^{'2}}{c ^{2}}} (v^{'} \cdot q E^{'}) \leftrightarrow f_{w} = \frac{1}{c 1 - \frac{v ^{2}}{c ^{2}}} (v \cdot q E)$

重力と慣性（等価原理）

　特殊相対性理論において、各地点の座標と時間について議論していたと思うが、一定速度で動く観測者というのは実際にないことが多い。というのも、地球上で観測を行うにあたって重力や地球の自転といった影響を受けていることからも伺える。一例として、等速円運動で自転している地球上で球の運動を見てみると $m \frac{d ^{2} z}{d s ^{2}} = - m g$ であることから、外部から見ると $m \frac{d ^{2} z ^{'}}{d s ^{'2}} = m r^{2} ω^{2} - m g$ というように慣性力が働く。ただ、ここでこの式を $m r^{2} ω^{2} - m g = - m (g - r^{2} ω^{2}) = - m g^{'}$ というようにすると、あたかも重力が変化したように見ることもできる。これ以外にも、加速する系では慣性力と重力が働いているが、互いに等価的なものとしてみなすことができる。実際、Newtonの運動方程式において $m_{I} \frac{d ^{2} r}{d s ^{2}} = F$ と書かれるが、ここに現れる質量は物体の慣性の度合いを表すものであるので慣性質量と呼ばれる。必要なときには添え字 $I$ をつけて表す。 Newton力学にはもう一つ別の種類の質量が登場する。それは万有引力の法則 $F = G \frac{m _{1 G} m _{2 G}}{r ^{2}}$ に使われる質量である。これは物体の慣性とは概念として全く別のものであり重力質量と呼ばれる。普通、添え字 $G$ をつけて慣性質量と区別する。上記の式によれば、地球表面で重力質量 $m_{G}$ の物体が受ける重力は $F = m_{G} g$ となる。ここで $g$ は、 $M_{G}, R$ を地球の重力質量、半径として $g = G \frac{M _{G}}{R ^{2}}$ で与えられる。さて、地球上での質点の自由落下の方程式というのは、簡単に鉛直方向 $z$ だけを考慮すると $m_{I} \frac{d ^{2} r}{d s ^{2}} = - m_{G} g$ となる。この式の両辺を $m_{I}$ で割ると $\frac{d ^{2} r}{d s ^{2}} = - \frac{m _{G}}{m _{I}} g$ を得る。すなわち自由落下の加速度は比 $m_{G} / m_{I}$ に比例する。ところがGalileiの有名なPisaの斜塔の実験以降、他の摩擦力などなければ自由落下の加速度は全ての物体で等しいことが知られている。これは $m_{G} / m_{I}$ が全ての物体について等しい値を持つことを意味する。

　斜塔の落下実験はあまり精度のよいものとはいえないであろうが、振り子の実験などはもっと精密な測定ということができよう。2種類の質量の比例関係についてはその後さらに正確な実験が行われた。特に19世紀末から20世紀末にかけてE $\overset{o}{¨}$ tv $\overset{o}{¨}$ sが優れた実験を行っている。それによれば $m_{G} / m_{I}$ を多くの物体について測定した結果、 $1 0^{- 9}$ という実験制度の範囲ですべて同じ値になった。さらにもっと最近のDicke、Braginskiらの測定によると $1 0^{- 12}$ くらいの精度で $m_{G}$ と $m_{I}$ の比例関係が成り立っていることが確認されている。2種類の質量が正確に比例しているのであれば、その比の値を1に選ぶのが便利である。もし比の値が1でなかったならば、どちらかの質量測定の目盛りを変更して比が1になるようにすればよい。要するにすべての物体について $m_{G} = m_{I}$ と考えればよいのである。

　以上の議論から、等加速度運動する系も含めた変換を考慮することで重力も含めた方程式を導出できるものと考えられるため、まずは等加速度運動する座標系でどうなるかを見ていくことにする。

特殊相対論から一般相対論へ

　ここからは、等加速度運動している系同士での変換がどのようになるか見ていこう。例えば、今 $x$ 方向に一定加速度 $a$ で移動しているものとすると、

$\frac{d ^{2} x}{d t ^{2}} = a 、 \frac{d ^{2} x ^{'}}{d t ^{'2}} = 0 、 t^{'} = t$

となるため、以下の関係が成り立つものと考えられる。

$x^{'} = x - \frac{1}{2} a t^{2} = x - \frac{a}{2 c ^{2}} w^{2} 、 w^{'} = w 、 (w = c t)$

これをLorentz変換のところでも見せた時空図で描写すると以下のようになる。

このことから、加速する座標への変換というのは曲がったものになることが予想される。この曲がった座標へ変換する理論としてRiemann幾何学というものがある。だが、この分野はかなり難解であるため、まずRiemann幾何学の記法について述べていくことにする。

　一般的にベクトルは $x$ というように太字で表記されるが、用いる座標系によってその中身の表記が異なってくる。簡単な例として、2次元直交座標系の場合だと $(X, Y)$ というようになるが、2次元斜交座標系を用いると $(x, y)$ というように表記される。

ここで注意することは、

x = (X, Y) = (x, y)

というような等式が成り立つわけではなく、比較する座標系の関係で変化する。このことは、上記の座標間にて以下の関係が成り立つことからもうかがえる。

$X = x + \frac{1}{2} y 、 Y = \frac{1}{2} y \leftrightarrow x = X - Y 、 y = 2 Y$

この関係により、各座標系でのベクトルは

$x = (X, Y) = (x + \frac{1}{2} y, \frac{1}{2} y) （直交座標系）$ $x = (x, y) = (X - Y, 2 Y) （斜交座標系）$

というようになるわけだが、この表記だと一目見てどの座標系を利用しているのかが分からず不便である。そこで、どの座標系を利用しているか分かりやすいよう各座標ごとに単位ベクトル（各成分の大きさ1のベクトル）というものを用意する。例として、直交座標系での単位ベクトルを $e_{X}, e_{Y}$ 、斜交座標系での単位ベクトルを $e_{x}, e_{y}$ として

$x = X e_{X} + Y e_{Y} = (x + \frac{1}{2} y) e_{X} + (\frac{1}{2} y) e_{Y}$ $x = x e_{x} + y e_{y} = (X - Y) e_{x} + 2 Y e_{y}$

というように表記すると、それぞれ直交座標系と斜交座標系を利用していることが分かる。

更に言うと、上記のベクトルの関係は

$x = X e_{X} + Y e_{Y} = x e_{x} + y e_{y}$

とも書けるわけだが、もし $e_{X}, e_{Y}$ がそれぞれ $e_{x}, e_{y}$ により置き換えることができれば直交座標系から斜交座標系へ変換ができることになる。そこで、単位ベクトル同士の関係を見るために、一般的に

$x (X, Y) = x (x, y)$

というように直交座標系と斜交座標系で表したものが等しいものとする。そして微小変化をとることで $d x = u_{X} d X + u_{Y} d Y = u_{x} d x + u_{y} d y$ $(u_{X} = \frac{\partial x}{\partial X} 、 u_{Y} = \frac{\partial x}{\partial Y} 、 u_{x} = \frac{\partial x}{\partial x} 、 u_{y} = \frac{\partial x}{\partial y})$

と展開できる。すると、前の式と同じように各成分のベクトルの和になっていることが分かる。そこで、 $u$ を基底ベクトル（大きさが1でない各成分のベクトル）と考え

$e_{X} = \frac{u _{X}}{∣ u _{X} ∣} 、 e_{Y} = \frac{u _{Y}}{∣ u _{Y} ∣} 、 e_{x} = \frac{u _{x}}{∣ u _{x} ∣} 、 e_{y} = \frac{u _{y}}{∣ u _{y} ∣}$

という関係を満たしているとすると、まず各基底ベクトルは $e_{X}, e_{Y}$ により

$u_{X} = e_{X} 、 u_{Y} = e_{Y} 、 u_{x} = e_{X} 、 u_{y} = \frac{1}{2} (e_{X} + e_{Y})$

となる。次に、基底ベクトルの大きさについては直交座標系において単位ベクトル同士の内積が

$e_{X} \cdot e_{X} = 1 、 e_{X} \cdot e_{Y} = 0 、 e_{Y} \cdot e_{X} = 0 、 e_{Y} \cdot e_{Y} = 1$

と設定されていることから、各成分の基底ベクトルの大きさも求められる。

$∣ u_{X} ∣ = u_{X} \cdot u_{X} = 1 、 ∣ u_{Y} ∣ = u_{Y} \cdot u_{Y} = 1 、$ $∣ u_{x} ∣ = u_{x} \cdot u_{x} = 1 、 ∣ u_{y} ∣ = u_{y} \cdot u_{y} = 1$

したがって、各成分の単位ベクトルは $e_{x} = e_{X} 、 e_{y} = \frac{1}{2} (e_{X} + e_{Y})$ と言う関係になっていることが分かる。そのため、この関係を先ほどの式に代入してみると

$X = x + \frac{1}{2} y 、 Y = \frac{1}{2} y$

というようになるが、先ほどの変換式と見比べてみることで全て0になることが分かる。 $r_{θ} = r_{ϕ} = 0$

各座標系の基底ベクトルの関係とそれらの内積が分かると座標系間の関係が導かれることになる。このことをふまえて、今度は前回までに登場した変換（Galielei変換、Lorentz変換）を見ていくと、まずGaliei変換は $w^{'} = w 、 x^{'} = - βw + x 、 y^{'} = y 、 z^{'} = z 、 (β = \frac{V}{c})$ という変換であったため、ベクトル $x$ は以下の通りに書ける。 $x = w e_{w} + x e_{x} + y e_{y} + z e_{z} = w^{'} e_{w} + (x^{'} + β w^{'}) e_{x} + y^{'} e_{y} + z^{'} e_{z}$ 一方で、変換先の座標を用いて $x = w e_{w} + x e_{x} + y e_{y} + z e_{z} = w^{'} e_{w}^{'} + x^{'} e_{x}^{'} + y^{'} e_{y}^{'} + z^{'} e_{z}^{'}$ となるが、まず基底ベクトルの方は以下の通りになる。 $u_{w}^{'} = \frac{\partial x}{\partial w ^{'}} = e_{w} + β e_{x} 、 u_{x}^{'} = \frac{\partial x}{\partial x ^{'}} = e_{x} 、 u_{y}^{'} = \frac{\partial x}{\partial y ^{'}} = e_{y} 、 u_{z}^{'} = \frac{\partial x}{\partial z ^{'}} = e_{z}$ ここで、基底ベクトルの内積が次のように設定されているものとする。 $e_{w} \cdot e_{w} = 1 、 e_{x} \cdot e_{w} = 0 、 e_{y} \cdot e_{w} = 0 、 e_{z} \cdot e_{w} = 0$ $e_{w} \cdot e_{x} = 0 、 e_{x} \cdot e_{x} = 1 、 e_{y} \cdot e_{x} = 0 、 e_{z} \cdot e_{x} = 0$ $e_{w} \cdot e_{y} = 0 、 e_{x} \cdot e_{y} = 0 、 e_{y} \cdot e_{y} = 1 、 e_{z} \cdot e_{y} = 0$ $e_{w} \cdot e_{z} = 0 、 e_{x} \cdot e_{z} = 0 、 e_{y} \cdot e_{z} = 0 、 e_{z} \cdot e_{z} = 1$

すると、基底ベクトルの大きさは以下の通りになる。 $∣ u_{w}^{'} ∣ = u_{w}^{'} \cdot u_{w}^{'} = 1 + β^{2} 、 ∣ u_{x}^{'} ∣ = u_{x}^{'} \cdot u_{x}^{'} = 1 、$ $∣ u_{y}^{'} ∣ = u_{y}^{'} \cdot u_{y}^{'} = 1 、 ∣ u_{z}^{'} ∣ = u_{z}^{'} \cdot u_{z}^{'} = 1$

そのため、Galilei変換後の単位ベクトルは $e_{w}^{'} = \frac{1}{1 + β ^{2}} (e_{w} + β e_{x}) 、 e_{x}^{'} = e_{x}^{'} = e_{x} 、 e_{y}^{'} = e_{y} 、 e_{z}^{''} = e_{z}^{'}$ となることから、関係式は以下の通りになる。

$w = \frac{w ^{'}}{1 + β ^{2}} 、 x = x^{'} + \frac{β w ^{'}}{1 + β ^{2}} 、 y = y^{'} 、 z = z^{'}$

このように成立していないことが分かるため、内積の値に問題があることが予想される。そのため、

$e_{w} \cdot e_{w} = 1 、 e_{x} \cdot e_{w} = - \frac{1}{2} β 、 e_{y} \cdot e_{w} = 0 、 e_{z} \cdot e_{w} = 0$ $e_{w} \cdot e_{x} = - \frac{1}{2} β 、 e_{x} \cdot e_{x} = 1 、 e_{y} \cdot e_{x} = 0 、 e_{z} \cdot e_{x} = 0$ $e_{w} \cdot e_{y} = 0 、 e_{x} \cdot e_{y} = 0 、 e_{y} \cdot e_{y} = 1 、 e_{z} \cdot e_{y} = 0$ $e_{w} \cdot e_{z} = 0 、 e_{x} \cdot e_{z} = 0 、 e_{y} \cdot e_{z} = 0 、 e_{z} \cdot e_{z} = 1$

というようにしてみると、一方の内積も

$e_{w}^{'} \cdot e_{w}^{'} = 1 、 e_{x}^{'} \cdot e_{w}^{'} = \frac{1}{2} β 、 e_{y}^{'} \cdot e_{w}^{'} = 0 、 e_{z}^{'} \cdot e_{w}^{'} = 0$ $e_{w}^{'} \cdot e_{x}^{'} = \frac{1}{2} β 、 e_{x}^{'} \cdot e_{x}^{'} = 1 、 e_{y}^{'} \cdot e_{x}^{'} = 0 、 e_{z}^{'} \cdot e_{x}^{'} = 0$ $e_{w}^{'} \cdot e_{y}^{'} = 0 、 e_{x}^{'} \cdot e_{y}^{'} = 0 、 e_{y}^{'} \cdot e_{y}^{'} = 1 、 e_{z}^{'} \cdot e_{y}^{'} = 0$ $e_{w}^{'} \cdot e_{z}^{'} = 0 、 e_{x}^{'} \cdot e_{z}^{'} = 0 、 e_{y}^{'} \cdot e_{z}^{'} = 0 、 e_{z}^{'} \cdot e_{z}^{'} = 1$

となり変換式を満たすことが分かる。さらには、このときの各座標は以下の図のような関係にある。

　次にLorentz変換の場合、変換式は以下の通りであった。 $w^{'} = γ (w - β x) 、 x^{'} = γ (- βw + x) 、 y^{'} = y 、 z^{'} = z 、 (γ = \frac{1}{1 - β ^{2}})$ このときの基底ベクトルを求めてみると $u_{w}^{'} = γ (e_{w} + β e_{x}) 、 u_{x}^{'} = γ (e_{x} + β e_{w})$ であり、このときも内積が不変な形になっていないため、こちらでも試しに $(e_{w}^{'} \cdot e_{w}^{'}) = (e_{w} \cdot e_{w}) 、 (e_{w}^{'} \cdot e_{x}^{'}) = (e_{w} \cdot e_{x})$ というように内積が不変になるところがあるか見てみると $∣ u_{w}^{'} ∣^{2} (e_{w} \cdot e_{w}) = γ^{2} [(e_{w} \cdot e_{w}) + 2 β (e_{w} \cdot e_{x}) + β^{2}]$ $∣ u_{w}^{'} ∣∣ u_{x}^{'} ∣ (e_{w} \cdot e_{x}) = γ^{2} [(e_{w} \cdot e_{x}) + (1 + β^{2}) (e_{w} \cdot e_{w}) + β]$ $∣ u_{x}^{'} ∣^{2} = γ^{2} [1 + 2 β (e_{x} \cdot e_{w}) + β^{2} (e_{w} \cdot e_{w})]$

、単位ベクトルの関係とそれらの内積を求めてみると $e_{w}^{'} = γ (e_{w} + β e_{x}) 、 e_{x}^{'} = γ (e_{x} + β e_{w})$ $e_{w} \cdot e_{w} = 1 、 e_{w} \cdot e_{x} = e_{x} \cdot e_{w} = - β = cos (\frac{π}{2} + α)$ $e_{w}^{'} \cdot e_{w}^{'} = 1 、 e_{w}^{'} \cdot e_{x}^{'} = e_{x}^{'} \cdot e_{w}^{'} = β = cos (\frac{π}{2} - α)$ となり、座標も記載すると以下の図（ローデル図）の通りになる。

　このように既存の座標系の単位ベクトル（

e_{x}, e_{y}, e_{z}

など）から基底ベクトル

u

とその内積

u \cdot u

から単位ベクトルを求めることで変換前と後の座標系がどのような関係にあるか分かることになる。ただ、そもそも基底ベクトルは微小ベクトル

d x

で関連つけられていたため、

w = x^{0} 、 x = x^{1} 、 y = x^{2} 、 z = x^{3} 、 u_{w} = u_{0} 、 u_{x} = u_{1} 、 u_{y} = u_{2} 、 u_{z} = u_{3}

とおき、

d x

を以下のように書く。

d x = μ \sum d x^{μ} u_{μ} = μ^{'} \sum d x^{μ^{'}} u_{μ^{'}}

(d s)^{2} = d x \cdot d x = μ, ν = 0 \sum 3 (u_{μ} \cdot u_{ν}) d x^{μ} d x^{ν} = μ^{'}, ν^{'} = 0 \sum 3 (u_{μ^{'}} \cdot u_{ν^{'}}) d x^{μ^{'}} d x^{ν^{'}}

ここで片方の内積が分かることで、もう片方の内積や返還前後の関係が分かるため、これはリーマン計量と呼ばれている。

　一方で、相対論的力学でもあったようにLorentz変換においては以下の形で不変となっていた。 $(d w)^{2} - (d x)^{2} - (d y)^{2} - (d z)^{2} = (d w^{'})^{2} - (d x^{'})^{2} - (d y^{'})^{2} - (d z^{'})^{2}$ そのため、これを内積として定義した空間（Minkofsky空間）を考えると、内積の形が不変となって便利であることがうかがえる。

これを見ても分かるように内積から求めると変換式がおかしくなってしまうため、逆に

あるいは、これをさらに扱い易くするよう添え字に行列の番号を振って $x^{'0} x^{'1} x^{'2} x^{'3} = α_{0}^{0} α_{0}^{1} α_{0}^{2} α_{0}^{3} α_{1}^{0} α_{1}^{1} α_{1}^{2} α_{1}^{3} α_{2}^{0} α_{2}^{1} α_{2}^{2} α_{2}^{3} α_{3}^{0} α_{3}^{1} α_{3}^{2} α_{3}^{3} x^{0} x^{1} x^{2} x^{3}$ とおくと、各成分ごと（ベクトルではない）に簡略化して書ける。 $x^{' μ} = ν = 0 \sum 3 α_{ν}^{μ} x^{ν} (μ = 0, 1, 2, 3)$ あるいは、今ここで $ν$ に対して総和をとっているが、行列においては下付き添え字（行）と上付き添え字（列）の掛け算は足し合わせることになるので、以下のように総和記号を省いた表記（Einsteinの縮約記法）がよく用いられている。 $x^{' μ} = α_{ν}^{μ} x^{ν} (μ = 0, 1, 2, 3)$ 一方で、Lorentz変換においても $w^{'} = γ (w - β x) 、 x^{'} = γ (- βw + x) 、 y^{'} = y 、 z^{'} = z 、 (γ = \frac{1}{1 - β ^{2}})$ となるため、同じように行列にしてみると $w^{'} x^{'} y^{'} z^{'} = γ - γ β 00 - γ β γ 00 00100001 w x y z$ であるため、先ほどと同様に $α$ を用いた表記をすることができる。ところが、つい先ほど扱った加速度の系の場合だと $w^{'} = w 、 x^{'} = - \frac{a}{2 c ^{2}} w^{2} + x 、 y^{'} = y, z^{'} = z$ というように、 $w^{2}$ が関わってきてしまうため、一概に同じような関係式で表せないことが分かる。そこで、一般的に以下のような関数で表されるものとする。 $w^{'} = w^{'} (w, x, y, z) 、 x^{'} = x^{'} (w, x, y, z) 、 y^{'} = y^{'} (w, x, y, z) 、 z^{'} = z^{'} (w, x, y, z)$ そして、これらの微小変化をとると $d w^{'} = (\frac{\partial w ^{'}}{\partial w}) d w + (\frac{\partial w ^{'}}{\partial x}) d x + (\frac{\partial w ^{'}}{\partial y}) d y + (\frac{\partial w ^{'}}{\partial z}) d z$ $d x^{'} = (\frac{\partial x ^{'}}{\partial w}) d w + (\frac{\partial x ^{'}}{\partial x}) d x + (\frac{\partial x ^{'}}{\partial y}) d y + (\frac{\partial x ^{'}}{\partial z}) d z$ $d y^{'} = (\frac{\partial y ^{'}}{\partial w}) d w + (\frac{\partial y ^{'}}{\partial x}) d x + (\frac{\partial y ^{'}}{\partial y}) d y + (\frac{\partial y ^{'}}{\partial z}) d z$ $d z^{'} = (\frac{\partial z ^{'}}{\partial w}) d w + (\frac{\partial z ^{'}}{\partial x}) d x + (\frac{\partial z ^{'}}{\partial y}) d y + (\frac{\partial z ^{'}}{\partial z}) d z$ であるのだが、これ書き直すと以下のような形で書けることが分かる。 $d w^{'} d x^{'} d y^{'} d z^{'} = \frac{\partial w ^{'}}{\partial w} \frac{\partial x ^{'}}{\partial w} \frac{\partial y ^{'}}{\partial w} \frac{\partial z ^{'}}{\partial w} \frac{\partial w ^{'}}{\partial x} \frac{\partial x ^{'}}{\partial x} \frac{\partial y ^{'}}{\partial x} \frac{\partial z ^{'}}{\partial x} \frac{\partial w ^{'}}{\partial y} \frac{\partial x ^{'}}{\partial y} \frac{\partial y ^{'}}{\partial y} \frac{\partial z ^{'}}{\partial y} \frac{\partial w ^{'}}{\partial z} \frac{\partial x ^{'}}{\partial z} \frac{\partial y ^{'}}{\partial z} \frac{\partial z ^{'}}{\partial z} d w d x d y d z$ これは先ほどの各成分ごとの表記と同様な形で $d x^{' μ} = (\frac{\partial x ^{' μ}}{\partial x ^{ν}}) d x^{ν} (μ = 0, 1, 2, 3)$ と記述できるため、先ほどの一定加速度での行列は次の通りになる。 $d x^{'0} d x^{'1} d x^{'2} d x^{'3} = 1 - a x^{0} / c^{2} 00 010000100001 d x^{0} d x^{1} d x^{2} d x^{3}$ さらに、ここでの偏微分の部分を分母の添え字を用いて簡潔に表現すると以下のようになる。 $e_{ν}^{μ} = \partial_{ν} x^{' μ} = (\frac{\partial x ^{' μ}}{\partial x ^{ν}}) \to d x^{' μ} = e_{ν}^{μ} d x^{ν}$ あるいは、行列の見方を変えると $d w^{'} d x^{'} d y^{'} d z^{'} = d w \frac{\partial w ^{'}}{\partial w} \frac{\partial x ^{'}}{\partial w} \frac{\partial y ^{'}}{\partial w} \frac{\partial z ^{'}}{\partial w} + d x \frac{\partial w ^{'}}{\partial x} \frac{\partial x ^{'}}{\partial x} \frac{\partial y ^{'}}{\partial x} \frac{\partial z ^{'}}{\partial x} + d y \frac{\partial w ^{'}}{\partial y} \frac{\partial x ^{'}}{\partial y} \frac{\partial y ^{'}}{\partial y} \frac{\partial z ^{'}}{\partial y} + d z \frac{\partial w ^{'}}{\partial z} \frac{\partial x ^{'}}{\partial z} \frac{\partial y ^{'}}{\partial z} \frac{\partial z ^{'}}{\partial z}$ というように、単位ベクトルの和の形にもなるため、次のように書くこともできる。 $d x^{'} = d x^{ν} e_{ν} (e_{ν} = \frac{\partial x ^{'}}{\partial x ^{ν}})$ そして、この大きさを求めるために内積をとると $d s^{2} = d x^{'} \cdot d x^{'} = (e_{μ} \cdot e_{ν}) d x^{μ} d x^{ν} = g_{μν} d x^{μ} d x^{ν} (g_{μν} = e_{μ} \cdot e_{ν})$ となるが、ここで現れる $g_{μν}$ が計量と呼ばれており、Riemann幾何学では重要な役目を果たすものとなっている。実際、この大きさを次の通りにすると互いに不変な形で記載することができる。 $d s^{2} = g_{μν}^{'} d x^{' μ} d x^{' ν} = g_{λ τ} d x^{λ} d x^{τ} 、（ g_{μν}^{'} = e_{μ}^{'} \cdot e_{ν}^{'} 、 g_{λ τ} = e_{λ} \cdot e_{τ} ）$ この関係が成り立つかどうかは、まず $d x^{' μ} d x^{' ν} = \partial_{λ} x^{' μ} d x^{λ} \partial_{τ} x^{' ν} d x^{τ} = \partial_{λ} x^{' μ} \partial_{τ} x^{' ν} d x^{λ} d x^{τ}$ と展開することができ、これに計量 $g_{μν}^{'}$ をかけると $g_{μν}^{'} d x^{' μ} d x^{' ν} = g_{μν}^{'} \partial_{λ} x^{' μ} \partial_{τ} x^{' ν} d x^{λ} d x^{τ}$ であり、同じように $g_{λ τ} d x^{λ} d x^{τ}$ も展開してみると $g_{λ τ} d x^{λ} d x^{τ} = g_{λ τ} \partial_{μ}^{'} x^{λ} \partial_{ν}^{'} x^{τ} d x^{' μ} d x^{' ν}$ というようになるため、以下の関係式が成り立ち互いに代入することで元の $d s^{2}$ の関係式を満たしていることが確認できる。 $g_{μν}^{'} = g_{λ τ} \partial_{μ}^{'} x^{λ} \partial_{ν}^{'} x^{τ} 、 g_{λ τ} = g_{μν}^{'} \partial_{λ} x^{' μ} \partial_{τ} x^{' ν}$ ここまでで、 $d s^{2}$ の式というは今までのベクトルの積の形をしていないように見えるが $d x_{μ} = g_{μν} d x^{ν} 、 d x_{ν} = g_{μν} d x^{μ}$ というように定義することで、次のように積の形で表記することもできる。 $d s^{2} = d x^{' μ} d x_{μ}^{'} = d x_{ν}^{'} d x^{' ν} = d x^{λ} d x_{λ} = d x_{τ} d x^{τ}$ また、 $g_{μν}$ の逆行列で $g^{μν}$ というように書くと $d x^{ν} = g^{ν λ} d x_{λ} 、 d x^{μ} = g^{μτ} d x_{τ}$ であるから、以下のように逆行列の性質を利用して等式が成り立つことが分かる。 $d x_{μ} = g_{μν} d x^{ν} = g_{μν} g^{ν λ} d x_{λ} = δ_{μ}^{λ} d x_{λ} = d x_{μ}$ $d x_{ν} = g_{μν} d x^{μ} = g_{μν} g^{μτ} d x_{τ} = δ_{ν}^{τ} d x_{τ} = d x_{ν}$

実は、この関係が成り立つことは相対論的力学でも出てきており、このときは $d t^{'} 1 - \frac{v ^{'} ^{2}}{c ^{2}} = d t 1 - \frac{v ^{2}}{c ^{2}}$ であったが、変位の形に整理して二乗にすると不変な形をしていることが分かる（Galilei変換、等加速度系の変換でも別の形で不変な形になる）。 $d w^{'2} - d x^{'2} - d y^{'2} - d z^{'2} = d w^{2} - d x^{2} - d y^{2} - d z^{2}$ そのため、このときの計量は以下の形をしているものと考えられる。 $g_{μν} = 1000 0 - 1 00 00 - 1 0 000 - 1$ しかし、このままだと各座標ごとの単位ベクトルが $e_{0} = γ - γ β 00 、 e_{1} = - γ β γ 00 、 e_{2} = 0010 、 e_{3} = 0001$ であることから、計量の値は $g_{00} = e_{0} \cdot e_{0} = \frac{1 + β ^{2}}{1 - β ^{2}} \neq = 1 、 g_{11} \dots$ というように別の値がでてきてしまう。そこで、 Lorentz変換の形を虚時間を用いて $x^{'0} = γ (x^{0} - i β x^{1}) 、 x^{'1} = γ (x^{1} + i β x^{0}) 、 x^{'2} = x^{2} 、 x^{'3} = x^{3} 、 (x^{0} = i w)$ というようにすると、単位ベクトルが $e_{0} = γ i γ β 00 、 e_{1} = - i γ β γ 00 、 e_{2} = 0010 、 e_{3} = 0001$ となり、このときの計量は以下の形になる。 $g_{μν} = 1000010000100001$ この場合だと、以下の形で不変な形になることが分かる。 $- d w^{'2} + d x^{'2} + d y^{'2} + d z^{'2} = - d w^{2} + d x^{2} + d y^{2} + d z^{2}$

以上のことから、座標に含まれた虚数によってあたかも計量が別の形になっているように見えていたことが考えられ、この見かけの計量はよくMinkofsky計量と呼ばれている。 $η_{μν} = 1000 0 - 1 00 00 - 1 0 000 - 1 または η_{μν} = - 1 000 010000100001$

このように、特殊相対性理論にもRieman幾何学との間に計量を通じて関係しているため、以降ではこれを用いて一般相対性理論について述べていくことにする。

測地線方程式

　特殊相対性理論において、運動方程式は $\frac{d p}{d τ} = f \leftrightarrow \frac{d p}{d t} = F 、 p = \frac{m _{0} v}{1 - \frac{v ^{2}}{c ^{2}}}$ という形で書けたわけなので、これを解析力学へ応用してみることにする。まず、解析力学において運動方程式はラグランジアン $L$ と一般化座標 $q$ を用いて以下のようになっていた。 $\frac{d}{d t} (\frac{\partial L}{\partial q ˙}) - (\frac{\partial L}{\partial q}) = 0$ これと先ほどの運動方程式を比較すると、力がポテンシャル $U$ で記述されるもの（保存力）としたとき $\frac{\partial L}{\partial v _{x}} = \frac{m _{0} v _{x}}{1 - \frac{v ^{2}}{c ^{2}}} 、 \frac{\partial L}{\partial v _{y}} = \frac{m _{0} v _{y}}{1 - \frac{v ^{2}}{c ^{2}}} 、 \frac{\partial L}{\partial v _{z}} = \frac{m _{0} v _{z}}{1 - \frac{v ^{2}}{c ^{2}}}$ $\frac{\partial L}{\partial x} = - (\frac{\partial U}{\partial x}) 、 \frac{\partial L}{\partial y} = - (\frac{\partial U}{\partial y}) 、 \frac{\partial L}{\partial z} = - (\frac{\partial U}{\partial z})$ であるから、ラグランジアンが以下の通りになることが分かる（静止エネルギーがポテンシャルエネルギーの一部のような振る舞いをする）。 $L = - m_{0} c^{2} 1 - \frac{v ^{2}}{c ^{2}} - U ≃ \frac{1}{2} m v^{2} - (m_{0} c^{2} + U)$ 次に、前回でも登場した計量 $g_{μν}$ により $1 - \frac{v ^{2}}{c ^{2}} = \frac{1}{c} (\frac{d s}{d t}) 、 \frac{d s}{d t} = η_{μν} \frac{d x ^{μ}}{d t} \frac{d x ^{ν}}{d t} \to g_{μν} \overset{x}{˙}^{μ} \overset{x}{˙}^{ν}$ と置き換えて、改めて運動方程式を記述すると $\overset{s}{˙} = g_{μν} \overset{x}{˙}^{μ} \overset{x}{˙}^{ν} 、 q = x^{λ}$ として $\frac{\partial L}{\partial x ˙ ^{λ}} = - m_{0} c \frac{\partial}{\partial x ˙ ^{λ}} g_{μν} \overset{x}{˙}^{μ} \overset{x}{˙}^{ν} = - \frac{m _{0} c}{2 s ˙} (g_{λ μ} \overset{x}{˙}^{ν} + g_{λ μ} \overset{x}{˙}^{μ}) = - m_{0} c (g_{λ μ} \frac{d x ^{μ}}{d s})$ $\frac{\partial L}{\partial x ^{λ}} = - m_{0} c \frac{\partial}{\partial x ^{λ}} g_{μν} \overset{x}{˙}^{μ} \overset{x}{˙}^{ν} - \frac{\partial U}{\partial x ^{λ}} = - \frac{m _{0} c s ˙}{2} (\frac{\partial g _{μν}}{\partial x ^{λ}} \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s}) - \frac{\partial U}{\partial x ^{λ}}$ となる（ここで $g_{μν} = g_{νμ}$ を利用している）。以上の結果を上記の運動方程式に代入すると以下のような形になる。 $m_{0} c \overset{s}{˙} [\frac{d g _{λ μ}}{d s} \frac{d x ^{μ}}{d s} + g_{λ μ} \frac{d ^{2} x ^{μ}}{d s ^{2}} - \frac{1}{2} (\frac{\partial g _{μν}}{\partial x ^{λ}} \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s})] - \frac{\partial U}{\partial x ^{λ}} = 0$ ここで、計量に関しての微分が $\frac{d g _{λ μ}}{d s} = \frac{\partial g _{λ μ}}{\partial x ^{ν}} \frac{d x ^{ν}}{d s}$ であるため、これを踏まえて整理すると次の通りになる。 $m_{0} c \overset{s}{˙} [g_{λ μ} \frac{d ^{2} x ^{μ}}{d s ^{2}} + (\frac{\partial g _{λ μ}}{\partial x ^{ν}} - \frac{1}{2} \frac{\partial g _{μν}}{\partial x ^{λ}}) \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s}] - \frac{\partial U}{\partial x ^{λ}} = 0$ 仮に外部からの影響がなく（ $U = 0$ ）、 $\overset{s}{˙} \neq = 0$ であることを踏まえると、 $g_{λ μ} \frac{d ^{2} x ^{μ}}{d s ^{2}} + (\frac{\partial g _{λ μ}}{\partial x ^{ν}} - \frac{1}{2} \frac{\partial g _{μν}}{\partial x ^{λ}}) \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s} = 0$ となるわけだが、ここで以下のように変形する。 $(\frac{\partial g _{λ μ}}{\partial x ^{ν}} - \frac{1}{2} \frac{\partial g _{μν}}{\partial x ^{λ}}) \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s} = \frac{1}{2} (\frac{\partial g _{λ μ}}{\partial x ^{ν}} + \frac{\partial g _{λ ν}}{\partial x ^{μ}} - \frac{\partial g _{μν}}{\partial x ^{λ}}) \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s}$ そして、計量の微分を展開してみると $\frac{\partial g _{λ μ}}{\partial x ^{ν}} = \frac{\partial}{\partial x ^{ν}} (e_{λ} \cdot e_{μ}) = \frac{\partial e _{λ}}{\partial x ^{ν}} \cdot e_{μ} + e_{λ} \cdot \frac{\partial e _{μ}}{\partial x ^{ν}}$ であり、ここで第１種クリストフェル記号を $Γ_{ν λ, μ} = \frac{\partial e _{λ}}{\partial x ^{ν}} \cdot e_{μ} = e_{μ} \cdot \frac{\partial e _{λ}}{\partial x ^{ν}}$ というように定義することで、 $\frac{\partial g _{λ μ}}{\partial x ^{ν}} = Γ_{ν λ, μ} + Γ_{νμ, λ}$ となる。同様に、その他の偏微分についても $\frac{\partial g _{λ ν}}{\partial x ^{μ}} = Γ_{μ λ, ν} + Γ_{μν, λ} 、 \frac{\partial g _{μν}}{\partial x ^{λ}} = Γ_{λ μ, ν} + Γ_{λ ν, μ}$ となるため、以下のような関係が成り立つことが分かる。 $Γ_{μν, λ} = \frac{1}{2} (\frac{\partial g _{λ μ}}{\partial x ^{ν}} + \frac{\partial g _{λ ν}}{\partial x ^{μ}} - \frac{\partial g _{μν}}{\partial x ^{λ}})$

そして、元の式に $g^{kλ}$ をかけることにより以下の形が得られる。 $\frac{d ^{2} x ^{k}}{d s ^{2}} + g^{kλ} Γ_{μν, λ} \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s} = 0$ これは第２種クリストッフェル記号を用いて $Γ_{μν}^{k} = g^{kλ} Γ_{μν, λ} = \frac{1}{2} g^{kλ} (\frac{\partial g _{λ μ}}{\partial x ^{ν}} + \frac{\partial g _{λ ν}}{\partial x ^{μ}} - \frac{\partial g _{μν}}{\partial x ^{λ}})$ により、次の測地線方程式の形で記載される。 $\frac{d ^{2} x ^{k}}{d s ^{2}} + Γ_{μν}^{k} \frac{d x ^{μ}}{d s} \frac{d x ^{ν}}{d s} = 0$ この式は物体がRiemann空間内を最短距離で通るということを前提にして求められるもので、例として点Aから点Bへ物体が移動するものとすると $s = \int_{A}^{B} d s = \int_{A}^{B} \frac{d s}{d t} d t = \int_{A}^{B} g_{μν} \overset{x}{˙}^{μ} \overset{x}{˙}^{ν} d t$ が最小となる（ $δs = 0$ ）ときの関係式からも確認できる。特殊相対性理論においては、 $g_{μν} = η_{μν}$ （Lorentz変換）とするとクリストッフェル記号は0となり $\frac{d ^{2} x ^{k}}{d s ^{2}} = 0 、 \frac{d s}{d t} = η_{μν} \overset{x}{˙}^{μ} \overset{x}{˙}^{ν} = c^{2} - v^{2}$ となるが、ここで $s$ の時間微分に関して固有時 $τ$ を $d τ = \frac{s ˙}{c} d t = \frac{d s}{c} = 1 - \frac{v ^{2}}{c ^{2}} d t$ とおくと、先程測地線方程式に $m_{0} c \overset{s}{˙}$ がかかっていたことを踏まえて $m_{0} c \overset{s}{˙} \frac{d ^{2} x ^{k}}{d s ^{2}} = 0 \to \frac{d p ^{k}}{d τ} = 0 、 (p^{k} = m_{0} \frac{d x ^{k}}{d τ})$ というように、相対論的力学で登場した運動方程式で力を0とした式が求められることが分かる。ここでの $p^{k}$ は4元運動量呼ばれており、右辺に力がある場合には $f^{k}$ （4元力）となる。このように、今までは速度に応じて質量が変化するものと思われてきたが、実は時間の進み方が変化していたことが分かる。

量子力学

　物理学の中に量子力学という分野ができたのは比較的新しいことである。それはPlanckが今世紀のはじめ（1900年）にいわゆるエネルギー量子というものを発見したときにはじまると考えてよい。前世紀の終わりごろまで物理学は主に目に見え耳で聞いて手に触れるような世界の自然現象をその対象にしていた。ところが世紀の終わりに近いころ、陰極線、X線或いは放射能などという珍しい現象が発見されるようになって、物理学の領域は直接我々の感覚に触れる世界から、分子とか原子とかという小さな世界に食い込んでいったのである。ところで、これらの小さな世界を支配している物理法則は当時のもので説明できるものであるか。実験物理学者はこの問いに答えるために、新しく発見された珍しい現象の研究に胸を躍らせながら従事した。一方で、18世紀後半にイギリスで産業革命が起こり鉄を精錬する溶鉱炉が飛躍的に進歩したが、よりよい鉄を作るためには温度を正確に測ることが重要であった。この温度を理論的に説明するため、一部の理論物理学者は「空洞輻射の問題」としてこの問いに答える手がかりを求めて、この問題を解くのに尽力することとなった。

　前世紀においても小さな（分子や原子の）世界の研究はMaxwellやBoltzmannによって大規模に行われていた。これらの人々は分子や原子の世界の法則も、ただ寸法が小さいだけで大きな世界のものそのままであると仮定して、遠い昔ギリシャの原子論者の想像していた原子論を数学的に築き上げようとしていた。そしてその努力は非常な成功をおさめたのである。この目に見える世界の法則というのはNewtonの力学とMaxwellの電磁気学のことである。そして分子や原子はこの力学に従って運動し、この電磁気学に従って光を射出したりするものと考える。目に見える物質は非常にたくさんのこれらの分子や原子からできていて、それらがこの法則に従って複雑に運動し、それらの統括的な効果が我々の目に見える物の性質として現れると考えるのである。

　この統括的な効果を取り扱うには統計的な方法を用いなければいけない。我々の目に見えるのは一つ一つの分子や原子の運動、それの出し入れする光ではなく、多くの分子や原子が行うそれらの現象のある種の統計的平均である。例えば、入れ物に入れられたたくさんのガス分子は、容器の中を飛び回るうちに絶えず壁にぶつかって、その度にこれに衝撃的な力を及ぼすものであるが、我々がガスの圧力として観測するのはこの力の平均値（この時には時間的平均値と考えてよい）である。また我々が温度として観測するのは、その分子の一自由度当たりの運動エネルギーの平均値なのである。

　ところで、力学と電磁気学の法則がそのまま小さな世界でも成り立つと仮定して統計的考察を行うこの原子論は、物質の色々な性質を説明するうえでかなり成功したけれども、空洞輻射の問題では全くうまくゆかない。その結果、小さい世界の物理法則としてNewton力学やMaxwellの電気学をそのまま持ってくることが疑われ始めた。すなわち、この空洞輻射の問題を解決するためには、どうしてもこれらの法則と異なったものを土台として新しい物理的世界像を築かねばならぬ結果となったのである。この新しい理論によりNewton力学やMaxwell電磁気学は古典理論という地位を有するようになり、後に量子力学として理論も発展していくのである。

空洞輻射の理論

　空洞輻射の問題とは一口に言えばこうである。

「熱せられた物体は温度ごとにどういう色の光を出すであろうか」

　物体を熱したとき出る光は、温度の低い時には赤い色をしているが、温度が上がるにつれて次第に白い色を呈するようになる。この現象を原子論的に説明するとどういうことになるかという問題である。この問題を物理学的に精密に言い表し、且つ取り扱いやすい形に理想化すると次のようにいってもよい。ある物体が壁で囲まれた室（空洞）の中にあったとする。この物体を熱すると光を発するのだが、その光により壁の温度も変化する。そして、空洞内が熱平衡になり温度 $T$ になり特定の光を生じるようになる。この光の強度（スペクトル）は温度 $T$ だけに関係し、壁の物質、空洞の形、あるいは大きさによらないことがまずKirchhoffによって証明された。

　次にスペクトルがどのように温度に依存するか考察するにあたって、統計力学において一自由度あたり $1/2 k T$ のエネルギーが分配されるという法則（エネルギー等分配則）が利用できると思われる。そこで、光が電磁波であり一方向にあたり電場と磁場によるエネルギーを1つずつ有すること、横波で2つの自由度が別に存在することから $2 k T (= 1/2 k T \times 4)$ 分のエネルギーを有すると考えてみる。また、スペクトルが空洞の形によらないことから一辺 $L$ の立方体の部屋を仮定して、この中で電磁波がどのようになるかというと固定端反射をすることから波長は以下のように制限される。 $λ = \frac{2 L}{n} (n = 1, 2, 3, \dots)$ 三次元の場合だと、振動数において $ν^{2} = ν_{x}^{2} + ν_{y}^{2} + ν_{z}^{2}$ というように各成分の二乗和で書けるため、先ほどの波長の式を $ν_{x} = \frac{c n _{x}}{2 L} 、 ν_{y} = \frac{c n _{y}}{2 L} 、 ν_{z} = \frac{c n _{z}}{2 L} (n_{x}, n_{y}, n_{z} = 1, 2, 3, \dots)$ としてこれらの式を代入すると以下のようになる。 $ν^{2} = (\frac{c}{2 L})^{2} (n_{x}^{2} + n_{y}^{2} + n_{z}^{2})$ 一見すると離散的な値をとっているため一概には計算ができないが、ここで $L \to \infty$ として振動数が連続していると仮定すると、振動数 $ν$ と $ν + d ν$ での自由度は $(\frac{c}{2 L})^{3}$ につき1自由度あることから以下のようになる。 $\frac{1}{8} \frac{4 π ν ^{2} d ν}{( \frac{c}{2 L} ) ^{3}} = \frac{4 π L ^{3}}{c ^{3}} ν^{2} d ν$ スペクトルが空洞の形によらないことを踏まえて、エネルギー密度を求めると以下のようになる。 $u d ν = \frac{1}{L ^{3}} (\frac{4 π L ^{3}}{c ^{3}} ν^{2} d ν \times 2 k T) = \frac{8 πk T}{c ^{3}} ν^{2} d ν$ これがRayleigh–Jeansの法則と呼ばれるものとなっている。実はこの式は低周波数でしか成立していないが、これは先ほど $L \to \infty$ としたことで波長、振動数が $λ \to \infty 、 ν \to 0$ となることからもうかがえる。一方で、高振動数の領域に対してはWienが提唱した以下の式（Wienの輻射式）が成り立つが知られている。 $u = \frac{8 πk α}{c ^{3}} e^{- αν / T} ν^{3} 、 (α ：定数)$ ここで式の導出過程を述べていないが、これは導出に用いられた仮定が不自然なものであり、なおかつMaxwellの電磁気学を全く無視していて説明するには不適切と判断したためである。

　以上のことを踏まえて、Planckは低振動数でRayleigh–Jeansの法則、高振動数でWienの輻射式が成り立つような内挿式（Planckの公式）を導出した。 $u = \frac{8 πk α}{c ^{3}} \frac{1}{e ^{αν / T} - 1} ν^{3}$ ここで、定数に関して $h = k α$ （Planck定数）に置き換えることで、以下のように $k T$ を含んだ形にまとめることができる。 $u = \frac{8 πh}{c ^{3}} \frac{1}{e ^{h ν / k T} - 1} ν^{3} 、（ h = 6.626 \times 1 0^{- 34} Js ）$ この式において、 $h ν << k T （ ν \to 0 ）$ とすると $u ≃ \frac{8 πh}{c ^{3}} \frac{1}{( 1 + h ν / k T ) - 1} ν^{3} = \frac{8 πk T}{c ^{3}} ν^{2}$ というようになり、反対に $h ν >> k T （ ν \to \infty ）$ では $e^{h ν / k T} >> 1$ であるため $u ≃ \frac{8 πh}{c ^{3}} \frac{1}{e ^{h ν / k T}} ν^{3} = \frac{8 πh}{c ^{3}} e^{- h ν / k T} ν^{3}$ と先ほどと同様な式が得られることが分かる。ここで、この式の意味を考えるため、 $u = \frac{8 π ν ^{2}}{c ^{3}} \frac{h ν}{e ^{h ν / k T} - 1}$ とすると、 $k T$ の代わりに別の形になっていることが分かる。この部分というのは、電磁波のハミルトニアンが振動体の形で書けることから自由度 $f$ に対して $H = i = 1 \sum f (a q_{i}^{2} + b p_{i}^{2})$ となり、この分の分配関数と期待値が $Z = \frac{1}{f ! h ^{f}} \int \dots \int e^{- β H} d q_{1} \dots d q_{f} d p_{1} \dots d p_{f} = \frac{1}{f !} (\frac{π}{h β ab})^{f}$ $⟨ H ⟩ = - \frac{\partial}{\partial β} ln Z = f k T 、 (β = \frac{1}{k T})$ であるため、 $k T$ が出てきているということになる（分配関数を除けば、1自由度の場合を $f$ 倍したものと同じ）。そのため、ここで1自由度あたりのエネルギーを $U_{1} = nh ν 、 (n : 自然数)$ という離散的な値をとる（エネルギーの最小単位いわゆる量子が $h ν$ である）と仮定すると $Z = n_{1}, \dots, n_{f} \sum W (n_{1}, \dots, n_{f}) e^{- (n_{1} + \dots + n_{f}) β h ν}$ となる。ここで各準位ごとに対象の数が $N_{0}, \dots, N_{n}, \dots$ だけあるとすると $n_{1} + \dots + n_{f} = N_{1} + 2 N_{2} + \dots + n N_{n} + \dots$ であり、状態数も区別がつかないことを考慮して $W (n_{1}, \dots, n_{f}) = \frac{f !}{N _{0} ! N _{1} ! \dots N _{n} ! \dots}$ となるため、多項定理より以下の通りになる。 $Z = N_{0} + N_{1} + \dots = f \sum \frac{f !}{N _{0} ! N _{1} ! \dots N _{n} ! \dots} e^{- (0 N_{0} + \dots + n N_{n} + \dots) β h ν} = (n = 0 \sum \infty e^{- n β h ν})^{f}$ このときは、分配関数も一体の分配関数を $N$ 乗したものとなっており、 $Z_{1} = n = 0 \sum \infty e^{- n β h ν} = \frac{1}{1 - e ^{- β h ν}}$ であるから、最終的なエネルギーの期待値は以下の通りになる。 $⟨ U ⟩ = - \frac{\partial}{\partial β} ln Z = - f \frac{\partial}{\partial β} ln Z_{1} = \frac{f h ν}{e ^{h ν / k T} - 1}$ そのため、エネルギーが1自由度ごとに等分配されるわけではなく、振動数に応じて分配されるエネルギーに制限があることになる。

光の粒子仮説

　前回は電磁波のエネルギーの最小単位が $h ν$ であると仮定することで、Planckの公式を導出していた。ただ、これはエネルギーが質量と等価であることを踏まえると、電磁波は $h ν$ のエネルギーを持つ粒子として考えることもできると予想される。それを用いてEinsteinは、光が $h ν$ のエネルギーを持つ粒子としてもふるまうもの（光量子仮説）として光電効果の説明を行った。　光電効果とは高い振動数の光線（紫外線域）いわゆるX線のようなものを金属などの物質の表面にあてると、そこから電子が飛び出してくる現象である。実験を行ったLenardにより、以下のことが確認されている。

飛び出る電子の一つ一つのエネルギーを測ってみると、それは照らす輻射線の強さに無関係である。
照らす輻射線の強さを大きくすると、飛び出る電子の個数（単位時間あたりの個数、出現頻度）が多くなるだけである。
飛び出る電子一つ一つのエネルギーは照らす輻射線の色に関係する。すなわちその輻射線の波長が短いほどエネルギーの大きな電子が飛び出してくる。

これに対して、先ほどの光量子仮説を適用すると電子を飛び出させるのに必要なエネルギーを $W$ （仕事関数）としたときに電子のエネルギーは $E = h ν - W$ と書ける。そのため、実験で得られた事実と矛盾しないことが分かる。（1・3：波長、振動数にしか依存しない、2：強さ（量子の個数）が大きくなれば電子の数も比例して多くなる）このことは、後にMillikanが行った実験においても確かめられている。

　また、光の粒子性が現れるもう一つの現象として1923年にComptonによって発見されたCompton効果がある（ここで発見されたという意味はその現象が光の粒子性によるもであるということが発見されたという意味で、現象自体は以前より知られていたものとなっている）。X線を物体にあてるとそれが四方に散乱されることはよく知られている。ところが、このとき単色のX線を用いて実験を行い、この散乱されたX線の波長を詳しく測ってみると、その一部は入射X線と同じ波長を持っているが、その外に入射するX線よりも長いものがあるという事実が知られていたのである。このX線の散乱について有名なThomsonの理論が存在していたが、散乱されたX線の波長は入射したものと同じである必要があり、事実とは異なるものとなっていた。そこで、Comptonは光の粒子説を用いてきわめて単純に説明できることを示したのである。具体的には、まず光は $h ν$ という量子としてみなせることを述べたが、同時に運動量も有している。これについては、空洞輻射において光が壁に反射する際に圧力を及ぼすことから伺える。実際、光が等方的に壁に当たる場合、圧力 $P$ はエネルギー密度 $u$ と $P = \frac{1}{3} u$ という関係にあり、Stefanの法則やWienのずれの法則を導くときの土台となっている。さらに、気体分子運動論にて圧力は運動量により $P = n ⟨ c_{⊥} p_{⊥} ⟩ = \frac{N}{V} ⟨ c_{⊥} p_{⊥} ⟩$ であり、 $h ν$ の粒子が $N$ 個ある場合は $u V = N h ν \to u = \frac{N}{V} h ν$ となることから、以下の式が得られる。 $⟨ c_{⊥} p_{⊥} ⟩ = \frac{1}{3} h ν$ さらに、運動量と速度の方向が同じであることと等方的にエネルギーが利用されることを踏まえて $c_{⊥} = c (\frac{p _{⊥}}{p}) 、 ⟨ p_{⊥}^{2} ⟩ = \frac{1}{3} ⟨ p^{2} ⟩$ となることや、今は振動数 $ν$ のみ考えていることから $\frac{1}{3} c p = \frac{1}{3} h ν \to p = \frac{h ν}{c}$ となることが分かる。このことについては、特殊相対性理論からも推測できる。このときは、物体が光速では質量を持てない（静止質量が $0$ ）ことからエネルギーと運動量の関係式より $(h ν)^{2} = (m_{0} c^{2})^{2} + (c p)^{2} \to p = \frac{h ν}{c}$ というようになる。

運動量保存の法則により、 $\frac{h ν _{0}}{c} = \frac{h ν _{ϕ}}{c} cos ϕ + γ m_{e} v cos α 、 (γ = \frac{1}{1 - ( v / c ) ^{2}})$ $0 = \frac{h ν _{ϕ}}{c} sin ϕ - γ m_{e} v sin α$ であることから、 $(\frac{h ν _{0}}{c} - \frac{h ν _{ϕ}}{c} cos ϕ)^{2} = (γ m_{e} v cos α)^{2}$ $(\frac{h ν _{ϕ}}{c} sin ϕ)^{2} = (γ m_{e} v sin α)^{2}$ となるため、これら二式を足し合わせる。 $(\frac{h ν _{0}}{c})^{2} - 2 \frac{h ^{2} ν _{0} ν _{ϕ}}{c ^{2}} cos ϕ + (\frac{h ν _{ϕ}}{c})^{2} = (γ m_{e} v)^{2}$ さらに、エネルギー保存則から $h ν_{0} + m_{e} c^{2} = h ν_{ϕ} + γ m_{e} c^{2}$ であり、両辺を二乗すると $(h ν_{0} - h ν_{ϕ} + m_{e} c^{2})^{2} = (γ m_{e} c^{2})^{2}$ $v^{2} = c^{2} [1 - (\frac{m _{e} c ^{2}}{h ν _{0} - h ν _{ϕ} + m _{e} c ^{2}})^{2}]$ であるため、右辺は以下の通りになる。 $(γ m_{e} c^{2})^{2} = \frac{1}{c ^{2}} [(h ν_{0} - h ν_{ϕ} + m_{e} c^{2})^{2} - (m_{e} c^{2})^{2}]$ そのため、式を整理すると最終的に以下の式が得られる。 $ν_{ϕ} = \frac{ν _{0}}{1 + \frac{2 h ν _{0}}{m _{e} c ^{2}} sin ^{2} \frac{ϕ}{2}}$ これが散乱されたX線の振動数であり、 $ϕ$ によらず $ν_{ϕ} < ν_{0}$ と小さくなること及び実験とも一致した結果が得られる。ただ、この形よりも波長の形により $λ_{ϕ} = λ_{0} + (\frac{h}{m _{e} c}) 2 sin^{2} \frac{ϕ}{2}$ と表示した方が式が容易になり、ここでの $h / m_{e} c$ はCompton波長と呼ばれている。 $\frac{h}{m _{e} c} = 2.42 \times 1 0^{- 13} m = 0.0242 \overset{˚}{A}$

電子の波動説

　色々のガスのガイスラー管（蛍光灯の一種）や、いろいろの物質を電極にするアーク燈（ランプ）から出てくる光を分光器でみると、それぞれの物質に特有なスペクトルがみられる。このスペクトルのうち線スペクトルといわれているものは、その物質の原子が射出するものである。この線スペクトルは夫々の物質原子に配列を持っているから、原子の構造をうかがう最も直接なデータである。だからこのスペクトル線の配列の仕方の法則性を求めることは非常に重要なことである。

　いろいろな物質のうち水素の発するスペクトルは最も簡単な配列を示している。以下の図は可視部から紫外線部にかけて現れる水素のスペクトルの写真であるが、これを見ただけでもいかにも意味ありげな線の並び方になっている。

　この水素の線の並び方を数式的に表すことに成功したのはスイスの中学校の先生Balmerの手柄である（1885年）。彼は可視部にある4本のスペクトル線 H $_{α}$ 、H $_{β}$ 、H $_{γ}$ 及び H $_{δ}$ の波長がそれぞれ $\frac{9}{5} h 、 \frac{16}{12} h 、 \frac{25}{21} h 、 \frac{36}{32} h 、（ h = 3645.6 \overset{˚}{A} ）$ で表されることをを発見し、これが $λ = \frac{n ^{2}}{n ^{2} - 4} 、（ n = 3, 4, 5, 6 ）$ なる公式にまとめられることに気がついた。これに続いて、スウェーデンのRydbergも他の原子のスペクトルも細かく調べて、より複雑な状況でもいくつかの系列に分類すると、その一つ一つの系列に属する線が以下の公式に従って配列していることを1890年に発見した。 $\frac{1}{λ} = \frac{1}{λ _{\infty}} - \frac{R}{( n + b ) ^{2}} 、 n = 正の整数$ ただし、 $λ_{\infty}$ とはその系列の線の間隔が段々せまくなって行って、ついに一か所に集約するところ（系列端）の波長であり、上記の図においての右端 H $_{\infty}$ で示されたところの波長になる。他には、以下の図のようにNaのスペクトルは水素のものと比べて相当複雑であるが、これらの線は主系列、第1副系列、第2副系列という風に分類することで、各系列が上記の公式に従って配置していることが分かる。

特に、公式中の $R$ は原子の種類によらない不変定数（Rydbergの定数）であり $R = 109677.691 cm^{- 1}$ であり、Balmerの $h$ との関係は $R = h /4$ である。これに対して、 $b$ は物質ごとに又は系列ごとに異なる値をとるが、一つの系列内では近似的に一定の値を持つ定数である。さらに、Rydbergの研究は進み、系列端においても規則正しく並んでいることから、上記の公式において $\frac{1}{λ _{\infty}} = \frac{R}{( m + a ) ^{2}} 、（ m : 正の整数）$ であることも見出した。そのため、公式はRydbergの公式として以下の形でまとめられる。 $\frac{1}{λ} = \frac{R}{( m + a ) ^{2}} - \frac{R}{( n + b ) ^{2}}$ この公式に従うと、各系列の波長は $⎩ ⎨ ⎧ ・ m = 1 （ L y man 系列、 1906 年）・ m = 2 （ B a l m er 系列、 1885 年）・ m = 3 （ P a sc h e n 系列、 1908 年）・ m = 4 （ B r a c ka tt 系列、 1922 年）$ というように対応しており、水素の場合で図にまとめると以下のようになる（ $m$ が4以上の波長は省略している）。

　このように、スペクトルの波長は連続しておらず離散的な値をとることが分かる。これについては光の粒子説で述べたことが関係していると考えられる。具体的には、空洞輻射での光のエネルギーが $E_{n} = nh ν$ となっていることから、原子においてもエネルギーが離散的な値 $W_{n}$ しかとらないものとする。また、光電効果において電子は光から $h ν$ のエネルギーを吸収して起こるわけだが、そのとき電子のエネルギーが増加したと考えると逆にエネルギーが減少して $h ν$ の光を放出することも予想される。この２つのことから、仮に原子のエネルギーが $W_{n}$ から $W_{m}$ に減少し、その分 $h ν$ の光が放出されたとすると $h ν = W_{n} - W_{m} \to ν = \frac{W _{n} - W _{m}}{h}$ という関係が成り立つと考えられる（Bohrの振動数関係）。すなわち、あるエネルギー $W_{n}$ にある間は安定した状態（定常状態）であり、ある特定の作用が起こることで $W_{n}$ から $W_{m}$ のような変化（量子的飛躍、状態の遷移）が起こるものとする。すると、水素原子の場合では $\frac{1}{λ} = \frac{R}{m ^{2}} - \frac{R}{n ^{2}} \leftrightarrow ν = \frac{R c}{m ^{2}} - \frac{R c}{n ^{2}}$ であるが、整数が逆数をとった形になっており、1項目が端の波長に対応したものであることを踏まえて $m < n$ であるから $\frac{1}{m ^{2}} > \frac{1}{n ^{2}}$ となるが、この形を $W_{n}$ とすると空洞輻射の場合と逆で整数が大きくなるにつれてエネルギーが減少する形になってしまうため、 $W_{n}$ を以下のように負符号をつけたものにしてみる。 $\frac{W _{n}}{h} = - \frac{R c}{n ^{2}} \to W_{n} = - \frac{R h c}{n ^{2}}$ すると、 $n$ が大きくなるにつれて0に近づいて連続的に分布するようになるが、ここで連続で分布するようになるというのは古典的な状態に相当していることから、振動数に関する式においてまず $m = n - τ$ とすると $ν = R c (\frac{n ^{2} - m ^{2}}{n ^{2} m ^{2}}) = R c \frac{( n + m ) ( n - m )}{n ^{2} m ^{2}} = R c \frac{τ}{n ^{3}} \frac{( 2 - τ / n )}{( 1 - τ / n ) ^{2}}$ であるから、 $n \to \infty (τ / n \to 0)$ とする近似的に以下の式が得られる。 $ν ≃ \frac{2 R c}{n ^{3}} τ = \frac{2}{R c h ^{3}} ∣ W_{n} ∣^{3/2} τ 、 (∣ W_{n} ∣ = - W_{n})$ 一方で、水素原子は当時電子が円運動をしているものと考えられており、このときのエネルギーは $W_{n} = \frac{1}{2} m v_{n}^{2} - \frac{1}{4 π ϵ _{0}} \frac{e ^{2}}{r _{n}}$ であり、円運動する条件として $\frac{m v _{n}^{2}}{r _{n}} = \frac{1}{4 π ϵ _{0}} \frac{e ^{2}}{r _{n}^{2}}$ を利用すると、速度と動径は以下の通りになる。 $v_{n} = \frac{2}{m} ∣ W_{n} ∣^{1/2} 、 r_{n} = \frac{e ^{2}}{8 π ϵ _{0}} ∣ W_{n} ∣^{- 1}$ さらに速度に関して $v_{n} = 2 π r_{n} ν_{n}$ とおけるため $\frac{2}{m} ∣ W_{n} ∣^{1/2} = 2 π \frac{e ^{2}}{8 π ϵ _{0}} ∣ W_{n} ∣^{- 1} ν_{n} \to ν_{n} = 2 \frac{8 ϵ _{0}^{2}}{m e ^{4}} ∣ W_{n} ∣^{3/2}$ そのため、ある範囲 $τ$ でエネルギーが変化したときに放出する光の振動数は次のようになる。 $ν = 2 \frac{8 ϵ _{0}^{2}}{m e ^{4}} ∣ W_{n} ∣^{3/2} τ$ 以上の考察から、次の関係式が成り立つことが予想される。 $\frac{1}{R c h ^{3}} = \frac{8 ϵ _{0}^{2}}{m e ^{4}} \to R = \frac{m e ^{4}}{8 ϵ _{0}^{2} h ^{3} c}$ 各定数に値を代入すると、 $R$ の値が既知のものと一致していることが分かる。また、円運動の際は角運動量 $L_{n}$ が保存しているので求めてみると、 $L_{n} = m v_{n} r_{n}$ であるから $r_{n}$ と $v_{n}$ が $v_{n} = \frac{2 R h c}{m} \frac{1}{n} = \frac{e ^{2}}{2 ϵ _{0} h} \frac{1}{n} 、 r_{n} = \frac{e ^{2}}{8 π ϵ _{0} R h c} n^{2} = \frac{ϵ _{0} h ^{2}}{πm e ^{2}} n^{2}$ となることから、以下のように簡単な形になることが分かる。 $L_{n} = m v_{n} r_{n} = \frac{h}{2 π} n$ というように、簡単な関係式が導かれる。この関係式はBohrの量子条件と呼ばれ、角運動量が前回のエネルギーと同じように飛び飛びになっていることを意味している。さらに、ここで $m v_{n} = p_{n}$ とすると $2 π r_{n} = λ_{n} n 、 (λ_{n} = \frac{c}{ν _{n}} = \frac{h}{p _{n}})$ と前回の運動量と振動数の関係からおいてみると、電子の軌道の円周がちょうど波長の整数倍になっている（端が閉じた波になっている）ことが考えられる。これは、少しでもずれると波同士で打ち消しあいなくなることから結果的に、この関係を満たす波しかないということになる。これは、de Brogileが提唱した物質波であり、電子も波としての性質を有することを暗に示したものとなっている。

波動力学

　前回まで電子が波としての性質を有することを述べたが、粒子とされていたものが波でないかという考えは過去に光学にもあった。光学においては光を粒子としてその経路により光の振る舞いを説明していたが、多くの実験も踏まえて光が波であることが実証された。このとき、光学でも幾何光学から波動光学というように移動したわけだが、今回も同じように力学でも粒子から波動に置きかえてみると新しい法則が見つかるのではないかということが推測される。

　幾何光学の法則が力学の法則と同じ形をしているということは最小作用の原理のような進んだ理論をまたずとも、すでにNewtonの時代に注目され、Newtonはこれによって光の粒子論を展開しようとしていたのである。最も簡単な具体的な例として光の屈折を考えてみよう。

上図に示すように、屈折率 $n_{1}$ の媒質と $n_{2}$ の媒質とが平面 $XX$ を境にして接続しているとする。1本の（単色）光線がAからBに進んでいくとき、その進路AOBは、Snellの法則 $\frac{sin θ _{1}}{sin θ _{2}} = \frac{n _{2}}{n _{1}}$ を満たす。ただし、 $θ_{1}$ は入射角、 $θ_{2}$ は屈折角である。これがこの場合の幾何光学の法則である。これをNewtonと同じように粒子的な立場に立つと次のように説明される。光の粒子が媒質Aから媒質Bに入射するとき、その境界面 $XX$ を通過する際、その面に直角な力の作用を受ける。その力は、面に直角であるので粒子の速度の法線成分を変化させる（図の場合では増加する）。同時に、接線成分は変化しないはずなので以下の関係が成り立つ。

$v_{1} sin θ_{1} = v_{2} sin θ_{2} \leftrightarrow \frac{sin θ _{1}}{sin θ _{2}} = \frac{v _{2}}{v _{1}}$

当時の考え方としては、位相速度 $u$ と速度 $v$ を同じものとして

$\frac{sin θ _{1}}{sin θ _{2}} = \frac{v _{2}}{v _{1}} = \frac{u _{2}}{u _{1}} = \frac{n _{1}}{n _{2}}$

が得られるため、Snellの法則を満たさず粒子説が否定されていた。ところが、特殊相対性理論により光を粒子として考えた場合は速度ではなく運動量でしか取り扱えないため

$p_{1} sin θ_{1} = p_{2} sin θ_{2} \leftrightarrow \frac{sin θ _{1}}{sin θ _{2}} = \frac{p _{2}}{p _{1}} = \frac{n _{2}}{n _{1}}$

を満たす条件を考えると、次の比例式（比例定数 $k$ ）が成り立つことが予想される。

$n = k p$

　一方で、波として扱った場合にも、Huygensの原理に従って以下の図のように考えることができる。

このとき屈折する前後の波長の関係は、媒質1の中での波の波長を $λ_{1}$ 、媒質2の中でのそれを $λ_{2}$ とすると、

$\frac{λ _{1}}{sin θ _{1}} = \frac{λ _{2}}{sin θ _{2}} \to \frac{sin θ _{1}}{sin θ _{2}} = \frac{λ _{1}}{λ _{2}}$

となる。そのため、屈折率とはその場所での波長に逆比例するものであるといってもよい。すなわち、比例定数 $k^{'}$ を用いて $n = \frac{k ^{'}}{λ}$ と表わされることになる。

　以上のように、粒子的な立場あるいは波動論的な立場で考えることで運動量と波長の逆数が関係していることがうかがえる。このことから、光以外のものについても一般的に二つの側面から考察が可能であると予想できる。すなわち、電流（電子の流れ）というようなものを波動論的立場によって取り扱うことで、波としてみた時の光と同じように干渉や回折を示すことはなかろうか。光の進行がこれら波動特有の現象によって幾何光学の法則から外れるように、電子の運動も干渉や回折によって以前までの物理法則から外れて奇妙な振る舞いをすることはないであろうか。こうして、幾何学に対して波動光学が成立したことに相当して、通常の力学に対して波動力学が成立すると考えなければならないことになる。そこで、Compton散乱の説明でも登場した光の関係式

$E = h ν 、 p = \frac{h ν}{c} = \frac{h}{λ}$

がその他の物質（電子など）でも成り立つものと仮定してみる。ここで、 $E, p$ は光子（光を粒子とみなしたときの名称）のエネルギーと運動量、 $ν, λ$ は光波（光を波とみなしたときの名称）の振動数と波長である。その上で、電子における屈折がどうなるかを考えてみると、エネルギー $E$ が保存されており、運動量 $p$ も電子の質量 $m$ と速度 $v$ により $p = m v$ と書けるため、位相速度 $u$ は以下の通りになる。

$u = ν λ = \frac{E}{p} = \frac{E}{m v} \propto \frac{1}{v}$

そのため、電子においても屈折に関してはSnellの法則を満たすことが分かる。もう少し簡単な例で自由粒子の場合だと、エネルギーが $E = m v^{2} /2$ と書けるため、このときの位相速度 $u_{p}$ は

$u_{p} = \frac{E}{p} = \frac{1}{2} v$

というようになる。これは、相対論的力学的に書いた場合にも

$u_{p} = \frac{E}{p} = \frac{c ^{2}}{v}$

という形になる。どちらにせよ電子波の位相速度は電子の速度と異なることになる。一方で、位相速度とは異なる群速度（波が重なったときの全体的な波の移動速度） $v_{g}$ は角振動数 $ω$ と波数 $k$ により

$v_{g} = \frac{\partial w}{\partial k}$

と書けるが、ここで各物理量が

$ω = 2 π ν = \frac{E}{ℏ} 、 k = \frac{2 π}{λ} = \frac{p}{ℏ} 、 (ℏ = \frac{h}{2 π})$

となるため、エネルギーが運動量で

$E = \frac{p ^{2}}{2 m} 、 E^{2} = (p c)^{2} + (m_{0} c^{2})^{2}$ と書けるため、群速度は以下の通りになる。

$v_{g} = \frac{\partial E}{\partial p} = \frac{p}{m} = v 、 v_{g} = \frac{\partial E}{\partial p} = c^{2} \frac{p}{E} = \frac{c ^{2}}{u _{p}} = v$

相対論的量子力学

$⟨ a ∣ b ⟩$

参考文献

金原　寿郎（1963）「基礎物理学」
藤原　邦男（1984）「物理学序論としての力学」
Richard P. Feynman（1965）「ファインマン物理学 Ⅰ力学」
押田　勇雄、藤城　敏幸（1998）「熱力学」
砂川　重信（1997）「電磁気学」
櫛田　孝司（1983）「光物理学」

ベクトルについて

大きさと向きを持つ量

これは主に点の位置関係を表すために用いられる。ただ、一点のことを示すために使うことができないことに気を付ける。

例）図のように、ABCの三人がいた場合

Aさんから見てBさんはこの方向にこれだけの距離のところにいる。
Bさんから見てCさんはこの方向にこれだけの距離のところにいる。
Aさんから見てCさんはこの方向にこれだけの距離のところにいる。

互いの位置関係は実際、2つのベクトルから判明する。仮にAさんがCさんの位置を知らず、1のBに2を聞くと、直ちに3も判明してしまう（1と2を組み合わせたものが3になる）。 $A B + BC = A C$

この式はベクトルの方向のみ関わるもので、長さは三角形を描くことから、余弦定理などより求める必要があることに注意！

また、Cの代わりにBにしても $A B = A C + CB$ が成り立ち、前の式に代入すると以下のようになる。 $A C + (CB + BC) = A C$ ここで括弧の部分は結局B自身の場所を自分自身が見ているものとなっており、 $CB + BC = 0 \to CB = - BC$ 符号が変わると、向きが反転することが分かる。

$x y$ 座標でベクトルを表記してみる。前回の例に出てきたA,B,Cが次のように配置されているとする。 $A = (0, 0) 、 B = (1, 1) 、 C = (2, 1)$

すると、ベクトルをどのようにするかというと、

x

方向と

y

方向それぞれのほうから見た長さが判明すればよい。

A B = (1, 1) 、 A C = (2, 1) 、 CB = (- 1, 0) CB = - BC とおけ、 BC = (1, 0) であるからこのようになる。

これらを先ほどのベクトルの式に代入すると

(2, 1) = (1, 1) + (1, 0) {・ 2 = 1 + 1 （ x 成分） ・ 1 = 1 + 0 （ y 成分）

(1, 1) = (2, 1) + (- 1, 0) {・ 1 = 2 + (- 1) （ x 成分） ・ 1 = 1 + 0 （ y 成分）

ここまではベクトルの方向だけ扱ってきたが、ここでは大きさ（長さ）の方を見ていくことにする。図のように、ベクトル $a = (a_{x}, a_{y})$ があるとき、 $a = ∣ a ∣$ （ベクトル $a$ の大きさ）としてピタゴラスの定理から $∣ a ∣^{2} = a_{x}^{2} + a_{y}^{2} \to ∣ a ∣∣ a ∣ = a_{x} a_{x} + a_{y} a_{y}$ となる。このときの各辺を見ると ${・左辺：ベクトルを長さをかけたもの・右辺：各成分同士をかけて足したもの$ となる。

そうすると、右辺において一つが別のベクト $b = (a_{x}, a_{y})$ になった場合は $∣ a ∣∣ b ∣ = a_{x} b_{x} + a_{y} b_{y}$ であると推測される。これを確かめるため、図を長さだけで書くと、次のようになる。 $⎩ ⎨ ⎧ ・ a = ∣ a ∣ = a_{x}^{2} + a_{y}^{2} ・ b = ∣ b ∣ = b_{x}^{2} + b_{y}^{2} ・ c = ∣ b - a ∣ = ∣ a - b ∣ = (b_{x} - a_{x})^{2} + (b_{y} - a_{y})^{2}$

一方で、余弦定理を用いると

c^{2} = a^{2} + b^{2} - 2∣ a ∣∣ b ∣ cos θ

であることから、次の関係式が成り立つ。

(b_{x} - a_{x})^{2} + (b_{y} - a_{y})^{2} = ∣ a ∣^{2} + ∣ b ∣^{2} - 2∣ a ∣∣ b ∣ cos θ

これを整理することで以下の式が導かれる。

a_{x} b_{x} + a_{y} b_{y} = ∣ a ∣∣ b ∣ cos θ

これは内積と呼ばれるもので

a \cdot b = ∣ a ∣∣ b ∣ cos θ

というように書くと、ベクトル同士の関係を角度により確認することができる。